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Abstract 



In this thesis, we explore three alternatives to the current paradigm of the stan- 
dard inflationary big bang scenario. The three alternative themes are spin torsion 
(or Einstein-Cartan-Kibble-Sciama) theories, extra dimensions (braneworld cos- 
mology) and changing global symmetry. In the spin torsion theories, we found new 
cosmological solutions with a cosmological constant as alternative to the standard 
scalar field driven inflationary scenario and we conclude that these toy models do 
not exhibit an inflationary phase. In the theme of extra dimensions, we discuss the 
dynamics of linearized scalar and tensor perturbations in an almost Friedmann- 
Robertson- Walker braneworld cosmology of Randall-Sundrum type II using the 
1+3 covariant approach. We derive a complete set of frame-independent equa- 
tions for the total matter variables, and a partial set of equations for the non-local 
variables, which arise from the projection of the Weyl tensor in the bulk. The latter 
equations are incomplete since there is no propagation equation for the non-local 
anisotropic stress. For the scalar perturbations, we found solutions that reveal 
the existence of new modes arising from the two additional non-local degrees of 
freedom (extra dimensions) under the assumption that the non-local anisotropic 
stress vanishes. Our solutions should prove useful in setting up initial conditions 
for numerical codes aimed at exploring the effect of braneworld corrections on the 
cosmic microwave background (CMB) power spectrum. For the tensor perturba- 
tions, we set out the framework of a program to compute the tensor anisotropies 
in the CMB that are generated in braneworld models. In the simplest approxi- 
mation, we show the braneworld imprint as a correction to the power spectra for 
standard temperature and polarization anisotropies and similarly show that the 
tensor anisotropies are also insensitive to the high energy effects. Finally in the 
theme of changing global symmetry, we constructed a bounded isothermal solution 
embedded in an expanding Einstein de Sitter universe and showed that there is a 
possible phase transition in the far future. 
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Chapter 1 
Prologue 



/ want to know how God created this world. I am not interested in this 
or that phenomenon, in the spectrum of this or that element. I want 
to know His thoughts; the rest are details. " 

- Albert Einstein 

"The future of the empirical world (or of the phenomenal world) is 
completely predetermined by its present state, down to its smallest de- 
tail. " 

- Karl Popper, "Conjectures and Refutations" 

"There is a coherent plan in the universe, though I don't know what 
it 's a plan for. " 

- Fred Hoyle 

§1.1 Introduction 

How did our Universe begin? Where did it come from and where is it heading 
to? How will our Universe end? These questions about the origin and the end 
of our observable Universe have fascinated many great thinkers from the ancient 
Greeks to the dawn of the 21st century. There is, perhaps, no question other 
than that of the origin of our Universe, which so transcends cultural and temporal 
divides, inspiring the ancient and modern philosophers, and the modern cosmolo- 
gists. At a fundamental level, we yearn to find an explanation of why there is a 
Universe, how it has come to assume the form which we observe it to take, and 
what are the fundamental laws and principles that govern its evolution. Although 
some of the above questions are philosophical and metaphysical, there are still 
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others which are within the boundaries of scientific enquiry. Of course, we are 
interested here in the latter rather than the former. 

How do we describe our physical universe? The global description of the ob- 
servable Universe rests on two important attributes, namely the the cosmological 
parameters from a physical theory and the irregularities or inhomogeneities 
which we observe. The cosmological parameters describe quantitatively the vari- 
ous features of the Universe, for example, the baryon density, fib and the age of the 
Universe, t - The irregularities or density perturbations describe the physics 
of the evolution of the matter components and how they are distributed irregu- 
larly throughout the Universe. Einstein's theory of general relativity is a theory 
of gravity which is used to study cosmology, and is thought to be an effective low 
energy theory from a more fundamental theory that we have no grasp of. The 
dynamics of the expanding Universe are characterized by the Hubble parameter 
that describes the expansion rate, Hq, and the spatial curvature, Qk- The cosmo- 
logical parameters are obtained from experiments, for example, the measurement 
of the cosmic microwave background (CMB) radiation, which can be characterized 
by a thermal distribution at a temperature of about 2.728 K. However, not all the 
cosmological parameters are measured accurately by one type of experiment. For 
example, there are experiments such as large galactic surveys and measurement 
of type I supernova that can give us other accurate cosmological parameters to 
compare to the CMB experiments. For the matter content of our Universe, we ob- 
serve the visible baryonic matter. However, there is circumstantial evidence which 
suggest that there exists a large and dominant amount of nonbaryonic matter (or 
dark matter) in our Universe. To explain the issue of how the large scale struc- 
ture of clusters and galaxies has emerged, we assume the existence of cold dark 
matter which comprises particles with negligible velocity. There could be other 
possibilities such as hot dark matter (particles that decouple when their velocities 
are relativistic) or something exotic. The final possibility with increasing support 
is that our Universe possesses a nonzero positive cosmological constant, A. 

With these features, we have a reasonable description of our Universe. It re- 
mains a question how these features come to be the way they are. The answer lies 
in the search for an underlying theory of gravity which works from the Planck scale 
to the present size of our observable Universe. At present, there exists no theory of 
gravity which could work at the Planck scale although there are many candidates 
for these theories, namely M/string-theory and loop quantum gravity. On 
the other hand, an effective description beyond the Planck scale is expressed ele- 
gantly with Einstein's general relativity. General relativity provides an reasonable 
idea of how the cosmological parameters and the density perturbations come about. 
The resulting cosmology from general relativity is the standard inflationary big 
bang cosmology which assumes that there is a beginning of our Universe. In the 
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next section, we will summarize the central tenets of the standard inflationary big 
bang cosmology. 



§1.2 The Standard Inflationary Big Bang Cos- 
mology 

In this section, we briefly sketch the theory of the standard hot big bang cos- 
mology which is presented in many standard texts [109, 127, 131, 153, 154, 155]. 
The central premises of the theory are: 

1. The Universe is isotropic and homogeneous on large scales based on the 
Copernican or cosmological principle which states that there are no 
priviledged observers in the Universe. In general relativity, the observable 
Universe is described by the Friedmann-Robertson- Walker (FRW) metric* : 



ds 2 = -dt 2 + a(tf 



<h '~ +r 2 (d9 2 + sin 2 Ode/) 2 ) 



1-Kr 2 



where K is the curvature constant which is +1 for a closed Universe, —1 for 
an open Universe and for a flat Universe. 

2. The Universe is expanding and obeys an empirical relation known as the 
Hubble law. An example of the Hubble law is the relationship between the 
"distance" to a galaxy 1 ', di and the observed redshift of a galaxy, z, which 
can be expressed by the following power series: 

z = H d L + ^(q - l)(H d L ) 2 + . . . , (1.2) 

where H is the present value of the Hubble constant (H = (a/a)) which 
characterizes the age of the Universe, and q is the deceleration parameter 
which measures the rate of slowing of the expansion*. 

3. The Universe began from a state of infinite (or near infinite) density and 
temperature and then cooled with the expansion of the Universe. 



*In this thesis, we adopt the choice of natural units and set H = c = G = 1. 
tAn example of a distance indicator is the absolute luminosity L, which is defined as the 
energy per unit time produced by the source in its rest frame. 

$qo is defined in terms of the scale factor of the FRW metric to be: 

_ a(t ) 

qa "ifojfl* (L3) 



where to denotes the present age of the Universe. 
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The big bang model makes accurate and scientifically testable hypotheses. The 
agreement with data from astronomical observations has given us considerable 
confidence in the model. The four key observational successes of the standard Hot 
Big Bang model are the following: 

• Expansion of the Universe: Our Universe began about fifteen billion 
years ago in violent "explosions" at every point of spacetime. In an early 
super-dense phase, every particle started rushing apart from every other 
particle. At lower redshifts z < 1, the linear relationship between and z 
is clear and convincing. One may determine Hq from this linear relationship 
by making use of galaxies at relatively low redshifts. This linear relation at 
low redshifts was first discovered by Hubble [94]. 

• Origin of the cosmic background radiation: About 3 x 10 5 years af- 
ter the Big Bang, the temperature of the Universe had dropped sufficiently 
for electrons and protons to combine into hydrogen atoms. This is known 
as the recombination era. After the recombination era, the cosmic mi- 
crowave background photons just stream toward the observer at the present. 
Hence the temperature differences on this surface of last scattering become 
the anisotropies in the CMB temperature we see today. Since the recombina- 
tion era, the radiation has propagated freely as it is unable to interact with 
the electrically neutral background matter (hydrogen atoms). The radiation 
constantly loses energy because its wavelength is stretched by the expansion 
of the Universe. The radiation temperature has fallen to approximately 2.728 
K in the present day. This background temperature was discovered by Pen- 
zias and Wilson [157] in 1965. Subsequently, the temperature anisotropies 
((AT/T) ~ 10 -5 ) were discovered by the COBE satellite in 1992 [177]. The 
large angle CMB anisotropies probe the fluctuations in the early Universe, 
possibly due to the result of quantum mechanical processes during an epoch 
of inflation. The large scale structures of our present Universe are generated 
by these primordial fluctuations due to gravitational instability. A detailed 
description of the physics of these CMB anisotropies is given in [92, 93]. 

• Nucleosynthesis of the light elements: After the big bang occurred, the 
matter which consisted of free neutrons and protons, was very dense and 
hot. The temperature of the matter fell and some of the nucleons combined 
to form the light elements, for example, deuterium^, helium-3, and helium-4. 
Theoretical calculations predicted that 25 per cent of the matter should be 
in the form of helium-4, a result which is in good agreement with current 
stellar observations. The heavier elements, of which we are partly made, 

§The formation of deuterium depends crucially on the baryon/photon ratio. 
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were created later in the interiors of stars and spread widely in supernova 
explosions. 

• Formation of galaxies and large-scale structure: The formation of 
galaxies and large-scale structure of our observable Universe, can be under- 
stood in the framework of the hot big bang model. At about 10,000 years 
after the big bang, the temperature had fallen to the point at which massive 
particles dominate the energy density of the Universe. The matter content 
of the Universe can affect structure formation by different processes which 
modify any primordial perturbations, for example, growth under self gravita- 
tion and effects of pressure and dissipative processes. As a result, we observe 
the large-scale structure of our Universe today. 

Although the Big Bang model is successful, it is also plagued by various prob- 
lems. Since there is no quantum theory of gravity, there exist the problem of initial 
conditions for the standard big bang cosmology. It does not explain why the Uni- 
verse we observe is nearly flat (the flatness problem). Another interesting question, 
is that the comoving distance over which causal interactions can occur before the 
release of the cosmic microwave background is considerably less than the comov- 
ing distance that the radiation traverses after decoupling. It would imply that the 
microwaves arriving from regions seperated by more than the horizon scale at last 
scattering which typically subtends about a degree, cannot have interacted be- 
fore decoupling. The big bang model does not offer a satisfactory explanation for 
the nearly homogeneous temperature of the cosmic microwave background in the 
Universe. It requires that the homogeneity must be part of the initial conditions 
and this is known as the horizon problem. Finally, for the big bang to begin at a 
very high temperature, there are possible unwanted relics forbidden by observation 
which would survive till today, for example, the gravitino in supergravity theories 
or the moduli fields in superstring theory. 

To deal with these inherent problems, an inflationary phase is incorporated 
into the standard big bang cosmology. The inflationary cosmology was proposed 
by Guth [78], Albrecht and Steinhardt [3] and Linde [129, 130]. It assumes that 
there is an epoch during which the scale factor of the Universe is accelerating. It is 
often described as a rapid superluminal expansion driven by a source, for example, 
a scalar field (or the inflaton field). A physical characterization of inflation is 
given by the following expression 

where 1 / (aH) is the comoving Hubble length which is the most important charac- 
teristic scale of the expanding Universe. The condition for inflation is that the co- 
moving Hubble length is decreasing with time. With inflation, the above mentioned 
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problems of the big bang cosmology are more or less cured (see [109, 127, 153, 155]). 
In the case of the flatness problem, inflation requires the condition that the crit- 
ical density Q is driven towards 1 rather than away from it. It would explain 
why the Universe is nearly flat. As for the horizon problem, we can allow our 
present observable Universe to emerge from a tiny region which was well inside 
the Hubble radius early on during inflation because of the large reduction in the co- 
moving Hubble length during the inflationary phase. Finally, the relic abundances 
can be reduced by the expansion during inflation if they are produced before the 
inflationary epoch. 

Another example of the many successes of inflation is that it provides the mech- 
anism for setting the scale-invariant spectrum of density perturbations present in 
the CMB anisotropies. The discovery of the accoustic peaks in CMB anisotropics 
came from the recent balloon-based experiments such as BOOMERANG [115] and 
MAXIMA [14]. These accoustic peaks probe the information of the early Universe, 
since they depend mainly on the spectrum of initial fluctuations and fundamen- 
tal cosmological parameters. Although inflation is a successful phenomenological 
theory that accounts for the large scale structure and various other features of 
the observable Universe, the inflaton field (or the scalar field) is not related to 
the Standard Model which has unified electromagnetism, strong and weak nuclear 
interactions in high energy theories. As a result, there are many different models 
of inflation (see [109, 127, 132]). The future experiments in the CMB and grav- 
itational waves will be able to constrain the possible type of inflationary models 
that give us the present description of our observable Universe. 

Of course, the above summary does not do enough justice to the development 
of cosmology for the last few decades. However, we can formulate an updated 
picture of the history of our Universe (see fig. 1.1). In the picture, we could see 
that considering the period after inflation, we can account for most of the features 
of our Universe with the present standard inflationary big bang cosmology. 



§1.3 Three Alternatives to Standard Cosmology 

The central aim of this thesis is to look at three alternative themes to the 
present cosmology. Although there exist recent alternative models of inflation 
[7, 13, 17, 72, 90, 104] which are constructed in a top-down approach from a plau- 
sible theory of quantum gravity, the focus of this thesis is taking the themes in 
alternative cosmologies from the bottom-up approach (or phenomenological theo- 
ries) after the Planck scale. We adopt a modest approach in the subject and would 
like examine not all but some of these alternative models and their implications. 
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Figure 1.1: A diagram (not drawn to scale) illustrating the history of the 
standard inflationary big bang cosmology. 



We hope to explore phenomenological possibilities in these alternative models if 
possible. There are reasons why these alternative themes would be an interest of 
study (see fig. 1.2 to see how they fit into the big picture). The first reason is 
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to explore alternative mechanisms to the standard inflationary big bang cosmol- 
ogy. One could ask whether we could substitute other possible inflationary driven 
fields that might give us an alternative to the standard scalar fields, for example, 
this theme will be explored in the spin-torsion or Einstein-Cartain-Kibble- 
Sciama (ECKS) theories. The second reason is that some of these alternative 
themes are required as basic premises to be integrated into theories which are 
plausible candidates for quantum gravity theories. For example, the premise of 
extra-dimensions is an essential feature in M/String theory. Finally, we explore 
the topic of changing global symmetry where we look at the possibility of our 
Universe undergoing a phase transition. We will discuss these topics in detail in 
the following sections. 

§1.3.1 Spin- Torsion Theory 

The concept of torsion is introduced as an extension to Einstein's theory of 
relativity for various reasons [70]. These theories of gravitation included with 
torsion are built upon Cartan's earlier suggestion in [28] that the torsion (the anti- 
symmetric part of the connection) should be identified as a possible physical field. 
The first reason, perhaps a historical one, is that it creates an analogy of general 
relativity with the theory of dislocations for the study of materials in continuum 
mechanics. Yet, there is no reason to believe that the structures are inherently 
similar for both general relativity and continuum mechanics. The next and perhaps 
more plausible reason is that torsion theories are good candidates to formulate 
gravity as a gauge theory, for e.g. [38, 56, 81, 97, 98, 105, 142, 185]. Proceeding 
from the premise that torsion theories are possible candidates of a gauge theory of 
gravity, one possible extension of general relativity, known as spin-torsion theory 
can be found by the gauge-theoretic approach in [82]. In spin-torsion theory, the 
connection between torsion and quantum spin was identified in [105, 173, 192] 
when it was realised later that the stress energy tensor for a massive fermion field 
was not symmetric [43, 193]. 

In the presence of torsion, the standard cosmology would be significantly al- 
tered. The first results were discussed in [83, 103, 194]. The singularity theorems 
in the presence of torsion may be suppressed in a wide class of models [82]. In ad- 
dition, Kerlick [102] has shown that if a Dirac field provides the source of matter, 
then the energy condition in the singularity theorems is weakened by the presence 
of torsion, leading to an enhanced singularity formation rate. Hence it is of con- 
siderable interest to study cosmological solutions and their implications in such 
scenarios. 

Recently, a gauge theory of gravity (GTG) has been developed by Lasenby, 
Doran and Gull [121]. The gauge theory of gravity is based on the reformulation 
of physics in Clifford algebra pioneered by Hestenes and his collaborators in [85, 
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Figure 1.2: A diagram (not drawn to scale) illustrating how the three alter- 
native themes to the standard inflationary big-bang cosmology is altered with 
respect to fig. 1.1. 



86, 87, 88]. In the GTG approach, the gravitational effects are described by a 
pair of gauge fields defined over a flat Minkowski background spacetime. These 
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gauge fields ensure the invariance of the theory under arbitrary local displacements 
and rotations in the background spacetime, and the physical predictions are all 
extracted in a gauge invariant manner, i.e. the background spacetime does not 
play any dynamic role in the physics. In [121], the torsion sector was set to 
zero; generally, the GTG then reproduces most of the standard results in general 
relativity, but differs in the interpretation of global issues such as the role of 
topology and horizons. 

Further work has been done in [35, 49] to explore the torsion sector of GTG. 
One of the movitations was to look at spinor-driven inflationary models and see 
the comparisons with possible experiments and observations. Unlike in general 
relativity where torsion is interpreted as a property of a non-Riemannian manifold, 
the torsion in GTG is a physical field which is derived from the gravitational 
gauge fields. The added feature of GTG is that it constrains physically the form 
which the torsion can take. This constraint is a consequence of the demanding 
minimally coupled equations for the matter fields when derived from a minimally- 
coupled, gauge-invariant action. The spin-torsion theory in GTG led to interesting 
developments in the problem of finding cosmological solutions for a Dirac field 
coupled to gravity in a self consistent model. A new exact solution to spin torsion 
theory (or Einstein-Cartan-Dirac equations) was found in [35] that describes a 
homogeneous Universe with spin induced anisotropy. 

In Chapter 2, we extend the work in [35] to find the cosmological solutions for 
a Dirac field coupled self-consistently to gravity with a cosmological constant. We 
also discuss the implications of these cosmological solutions. 

§1.3.2 Extra Dimensions 

The holy grail of modern physics is to unify the four fundamental forces namely, 
gravity, electromagnetism, the strong and weak nuclear forces. Unfortunately, the 
two great pillars of theoretical physics, namely general relativity which is the 
theory of gravitation and quantum theory which describes the sub-atomic regime 
of matter, are mutually incompatible. The search for the theory of quantum gravity 
is still ongoing. 

The early Universe provides a testing ground for theories of gravity. The stan- 
dard cosmological model, based on general relativity with an inflationary era, is 
very effective in accounting for a broad range of observed features of the Universe. 
However, the lack of a consistent theoretical framework for inflation, together with 
the ongoing puzzles on the nature of dark matter and dark energy, indicate that 
cosmology may be probing the limits of validity of general relativity. 

M-theory (or superstring theory) is considered to be a promising potential path 
to quantum gravity. It assumes the basic premise that matter is made up of strings 
instead of point particles at the Planck scale (i.e. the scale where all the four forces 
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are unified). However, in order for M-theory to be a theory of quantum gravity, 
there are two pre-requisites, namely supersymmetry and extra dimensions. Both 
pre-requisites have not been experimentally verified. The Large Hadron Collider 
(LHC) in CERN (which will run in 2006) would test the possibility of finding 
supersymmetric particles at high energies. The question is whether we can find 
experiments that could demostrate the existence of extra dimensions 

As such, it is an important candidate for cosmological testing. In the absence of 
a sufficiently general M-theoretic model of cosmology, we can use phenomenological 
models that share some of the key features of M-theory, including branes. In brane 
cosmology, the observable Universe is a 1+3-dimensional "brane" surface moving 
in a higher-dimensional "bulk" spacetime. Standard-model fields are confined to 
the brane, while gravity propagates in the bulk. The simplest, and yet sufficiently 
general, phenomenological braneworld models are those based on the Randall- 
Sundrum II scenario [160]. These models have the additional advantage that they 
provide a framework for investigating aspects of holography and the AdS/CFT 
correspondence. 

The idea of using extra dimensions to unify fundamental forces was first ex- 
plored by Kaluza [99] and Klein [106] in the 1920s. Kaluza proposed that one 
can start from a five dimensional theory of general relativty, and assume that one 
of the dimensions is being curled up into a circle, leading to a theory that uni- 
fies gravity and electromagnetism. Klein took a step further after rediscovering 
Kaluza's theory, he noted the quantisation of the electric charge and hoped that 
Kaluza's theory would underlie quantum mechanics. He came close to the modern 
point of view in his discussion of higher harmonics and the size of the small circle, 
and insisted that the 5th dimension should be treated seriously. He assumed that 
the coordinate of the 5th dimension is periodic. It is difficult to envisage a 5D 
spacetime with such a topology. However, one can use a simple analogy provided 
by a hosepipe: at large distances, the hosepipe looks like a line but close inspection 
reveals at every point on the line, there is a circle (see Fig. 1.3). The Kaluza Klein 
picture requires the dimensions to be compact and essentially homogeneous. In 
this picture, it is the compactness which ensures that the spacetime is effectively 
4-dimensional at distances exceeding the compactification scale, i.e. the size of 
the extra dimensions. As a result, the size of the extra-dimensions must be very 
small, and the size was of the order of the Planck scale. With the Planck length 
lp ~ 10~ 35 m and the corresponding Planck energy scale Mp ~ 10 19 GeV, it remains 
almost impossible to probe the extra dimensions with present technology. 

In recent years, much attention has been focused on brane world scenairos. 
They are based on the idea that the standard model of particle physics (which 



^ For a more technical review on the subject of extra dimensions in the high energy physics 
framework, see [163]. 
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Figure 1.3: A diagram depicting the Kaluza-Klein concept of compactification 
of extra dimensions: (a) Compactifying a 3-D Universe with two space dimen- 
sions and one time dimension which simplifies the picture of Kaluza and Klein, 
(b) The Lorentz symmetry of the large dimension is broken by the compactifi- 
cation and what remains is the 2-D space plus the U(l) symmetry represented 
by the arrow, (c) On large scales we see only a 2-D Universe (one space plus 
one time dimension) with the "internal" U(l) symmetry of electromagnetism. 
The diagram is reproduced from the article "The search for extra dimensions" 
by Abel and Russell in [1]. 

unifies the electromagnetic, strong and weak nuclear forces) is confined to a 3- 
brane (an object or membrane which is a solution in string theory that corresponds 
to our observable Universe) and only gravity lives in the higher dimensional bulk. 
In general, the spatial dimensions transverse to the branes must be very small in 
order to avoid large deviations from Newton's law of gravitation (see 1.4). One 
interesting suggestion was put forward in [6] that two of the extra dimensions 
are large (mm scale). Since the four dimensional effective gravitational constant 
is related to the fundamental higher dimensional gravitational constant through 
the size of the extra dimensions, the fundamental gravitational constant can be 
reduced to the TeV scale, provided that the large extra dimensions have a size of 
~ 1 mm. The connection to superstring theory was made in [5] . This new proposal 
allows us to deal with the mass hierarchy problem in elementary particle physics 
without invoking supersymmetry. 

By projecting the brane world models to an effective 4D theory, we can study 
how the extra dimensions may affect the standard observational cosmology. Such 
an approach is the first step to understanding these new models and seeing whether 
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one can find any new features in them that may show up in future experiments in 
observational cosmology. In Chapter 3, we will set up the formalism of the (1+3)- 
covariant approach for studying cosmological perturbations in these braneworld 
models. Subsequently, in Chapter 4, we will explore the effects of the extra di- 
mensions on the cosmic microwave background for both scalar and tensor pertur- 
bations. 

§1.3.3 Changing Global Symmetry 

In standard general relativity, the Friedmann- Robertson- Walker (FRW) metric 
is a well-defined metric describing an expanding, homogeneous, isotropic Universe. 
The symmetry of a FRW spacetime is both translational and rotational about every 
point and there is no preferred centre (Copernican). 

However, local inhomogeneities grow to form galaxies and clusters of galax- 
ies and hence generate the large scale structure which we observe today. In an 
Einstein-de Sitter (flat, dust dominated) Universe, these clustered structures would 
grow continuously, leading to the merger of galaxy clusters to form superclusters. 
Subsequently, this leads to the formation of an ever-changing hierarchical struc- 
ture. Such an astrophysical mechanism results in a changing equation of state 
from a pressure-free expansion to a transition state where gravitational clustering 
induces a cosmologically significant pressure which feeds back into the metric. At 
this stage, the symmetry of the metric alters, and depends on the exact form of the 
equation of state. Extrapolating this picture into the future, we can envisage the 
following possible development. By taking the clustering model to the extreme, 
one could imagine a time in the far future when the whole Universe is basically 
one large cluster. The internal constituents of this cluster exchange energy. The 
whole Universe virializes on the largest possible scale and we get an asymptotic 
state represented approximately by an isothermal Universe. For such a state, the 
metric is well-defined, as the final symmetry is rotational about one point only 
and translational nowhere and there is a definite centre (Anti-Copernican). 

The effectiveness of gravitational clustering is related to the global structure 
of the standard Einstein-Friedmann Universes. There exist three general cases 
mentioned earlier for the future evolution of the Universe, and they depend on 
the total energy density. In a closed geometry the Universe will eventually rec- 
ollapse into a singularity and all large-scale structure will be destroyed in the 
big crunch. In an open geometry the rate of universal expansion is too great to 
allow complete clustering. In a flat geometry (Einstein-de Sitter) the expansion 
timescale is comparable to the clustering timescale. Gravitational clustering grows 
continuously. Therefore, although gravitational clustering of galaxies occurs in all 
Einstein-Friedmann cosmological models, only in the Einstein-de Sitter model does 
it grow continuously. Only the Einstein-de Sitter model allows the formation of 
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the isothermal asymptotic state. However, the existence of a cosmological constant 
would change this evolution. 

Saslaw and his collaborators[169, 171, 172] have found a solution to Einstein's 
field equations, which corresponds to a class of isothermal inhomogeneous Uni- 
verses in which the nonzero pressure balances gravity and discussed whether the 
infinite and unbounded isothermal sphere could provide a possible state of our 
Universe. They suggested various arguments which demonstrate the possibility 
that such spherical and static models might represent the ultimate state of an 
Einstein-de Sitter Universe (with A = 0). If that is the case, the Einstein-de Sitter 
Universe will undergo a dynamical symmetry breaking into the isothermal Uni- 
verse, i.e. from Copernican to anti-Copernican, which may correspond to a possible 
phase transition in the far future. This is (at least to our knowledge) the first time 
the question has been asked about what happens when a local condensation with 
one symmetry grows to eventually fill a Universe with another symmetry. 

In Chapter 5, we examine the changing global symmetry of the Einstein-de 
Sitter Universe and see how it could lead to a possible phase transition in the 
future. 
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Figure 1.4: The matter fields (the electromagnetic, strong and weak nuclear 
forces) in the standard model are confined on the 3-brane (our observable Uni- 
verse). Only gravity is allowed to propagate into the extra dimensions (i.e. 
the bulk). A collision between a proton and an antiproton can produce, for 
example, a single jet of matter particles plus graviton emission into the bulk. 
Such collisions might be seen in high-energy physics experiments. The diagram 
is reproduced from the article "The search for extra dimensions" by Abel and 
Russell in [1]. 



Part I 
Spin- Torsion Theories 



Chapter 2 

Spin- Torsion Cosmology 



"There could be no fairer destiny for any . . . theory than that it should 
point the way to a more comprehensive theory which it lives on as a 
limiting case. " 

- Albert Einstein 

In this chapter, we extend the gauge theory of gravity (GTG) approach of finding 
cosmological solutions in spin-torsion theories to include a cosmological constant. 
We find a new anti-particle solution in addition to the two solutions found in [35]. 
We also discuss the cosmological implications of these solutions. 

§2.1 Introduction 

The gauge theory of gravity (GTG) provides an alternative formulation to the 
standard spin-torsion or Einstein-Cartan-Kibble-Sciama theories. In the GTG ap- 
proach developed in [49, 120, 121], one can identify the torsion with the spin of the 
matter field that follows naturally from the gauge-theoretic approach to gravity 
(see [83]). Fundamentally in GTG, the gravitational effects are described by a pair 
of gauge fields which are defined over a flat Minkowski background spacetime. The 
gauge fields ensure the invariance of the theory under arbitrary local displacements 
and rotations in the background spacetime. All the physical predictions are ex- 
tracted in a gauge-invariant manner, thus ensuring that the background spacetime 
plays no dynamical role in the physics. In [49], the equations describing the Dirac 
field coupled self-consistently to gravity were given, and the analagoues of the 
Einstein-Cartan-Dirac (ECD) equations were derived from a minimially-coupled, 
gauged-invariant action. It is known that the Einstein-Dirac equations, which de- 
scribe a Dirac field coupled to gravity through the (symmetrised) stress-energy 
tensor only, cannot be derived from a minimally coupled action. Hence this serves 
as a compelling reason to regard the ECD equations as being more fundamental. 
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One of the main applications of GTG is in finding cosmological solutions for 
a Dirac field coupled to gravity in a self-consistent model. This problem was 
considered by several authors [2, 37, 46, 47, 84, 96, 152, 158, 159, 174] in the past. 
One of the main motivations was to see whether the Dirac field could offer an 
alternative type of inflationary model over standard scalar field models. However, 
the solutions found by these authors either solve the Einstein-Dirac equations 
[2, 37, 84, 152, 159] which do not include the effect of torsion induced by the spin 
of the Dirac field, or they do include the spin-induced torsion (i.e. they solve the 
Einstein-Dirac-Cartan equations) but only for massless fields [84, 96, 174] or for 
ghost solutions [46, 47, 158]. (A ghost solution has a vanishing stress energy tensor 
for the Dirac field.) The first massive, non-ghost solution was recently found for a 
Dirac field coupled self-consistent ly to gravity in [35]. 

The origin of a cosmological constant, A, is a mystery in gravitation and modern 
particle physics. It was introduced originally by Einstein to allow static homoge- 
neous cosmological solutions to the field equations in the presence of matter. The 
discovery of the expanding Universe by Hubble, made it unnecessary, but it did 
not change its status as a legitimate addition to the gravitational field equations, 
or as a parameter constrained by observation. In particle physics, the cosmolog- 
ical constant turns out to be a measure of the energy density of the vacuum i.e. 
the state of the lowest energy. This identification becomes crucial in dealing with 
two present unsolved problems: firstly, why the vacuum energy is so small and 
secondly, to understand why it is comparable to the present day density of the 
Universe. It turns out that there are compelling reasons to believe why our Uni- 
verse should have a non-zero A in both theory and experiment.* Recent results 
from the cosmic microwave background, supernova data and large scale structure 
seem to favour a cosmological constant with Q A = 0.7. 

In this chapter, we will extend the work in [35] to find a massive non-ghost 
solution for a Dirac field coupled self-consistently to gravity with the inclusion of a 
cosmological constant. We will briefly summarize the geometric algebra formalism 
and the gauge theory of gravity in §2.2. Subsequently, we will introduce the 
minimally coupled Einstein-Cartan-Dirac (ECD) action in §2.3, and the GTG 
analogue of the field equations. We will derive the cosmological solutions in §2.4 
and discuss their implications in §2.4.1. 



*A good review of the cosmological constant is given by Carroll [27]. 
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§2.2 Geometric Algebra and Gauge Theory of 
Gravity 

To find the cosmological solutions for a Dirac field coupled to gravity in a self 
consistent model with a cosmological constant, we shall employ the tools of Geo- 
metric Algebra (GA), which was originally developed from the spacetime algebra 
(STA) approach by Hestenes (see [85, 86, 87, 88]). We briefly summarize the 
geometric algebra formalism in §2.2.1 and the gauge theory of gravity in §2.2.2. 
Note that only in this chapter, we will use the other metric signature convention 

(+, - - -)■ 



§2.2.1 Geometric Algebra 

We start with the definition of a geometric (or Clifford) product in [39] that 
combines both symmetric inner product and antisymmetric outer product whuch 
is defined as follows: 

ab = a-b + aAb (2.1) 

where both a and b are vectors. 

In general, an arbitrary real superposition of the basis elements (2.7) is called a 
multivector. A multivector can be decomposed into sumns of elements of different 
grades, for example, scalars are grade zero objects, vectors grade one, bivectors 
grade two and etc. These elements inherit the associative Clifford product of the 
{7^} generators. For a grade-r multivector A r and a grade-s multivector B s , we 
define the inner and out products via 

A r .B s ={AB) lr _ s{ , (2.2) 
A r A B s = (AB) r+s , (2.3) 

where (M) r denotes the grade-r part of M. The subscript will be left implicit 
when taking the scalar part of a multivector. We shall make use of the commutator 
product 

A x B = -(AB-BA) . (2.4) 

The operation of reversion, denoted by a tilde, is defined by 

(AB)~ = BA , (2.5) 

and the rule that vectors are unchanged under reversion. We adopt the conven- 
tion that in the absence of brackets, inner, outer and commutator products take 
precedence over Clifford products. 
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The geometric algebra of spacetime is familiar to physicists in the form of the 
algebra generated from the Dirac-7 matrices. The spacetime algebra (STA) is 
generated by four vectors, {7^}, /x = 0...3, equipped with an associative (Clifford) 
product, denoted by juxtaposition. The symmetrised and antisymmetrised prod- 
ucts define the inner and outer products between vectors, and are denoted by a 
dot and a wedge respectively: 

7m • lv = \{l»lv + 7^7/0 = r)»v = diag(+ ) , 

2 (2.6) 

T/i A 7^ = 2 (T/iTv Ivlfi) ■ 

The outer product of two vectors defines a bivector — a directed plane segement 
including the two vectors. 

The full basis for a 16-dimensional STA is provided by 

1 {7^} {a k ,la k } {^7/J 1 

1 scalar 4 vectors 6 bivectors 4 trivectors 1 pseudoscalar 

(2.7) 

where a k = 7fc7o, k = 1. . .3, and X = 70717273 = o"iO"20"3- The pseudoscalar X 
squares to —1 and anticommutes with all odd-grade elements. The {a k } generate 
the geometric algebra of Euclidean 3-space, and are isomorphic to Pauli-matrices. 

We denote vectors in lower case Latin, a, or Greek for a set of basis vectors. A 
(coordinate) frame of vectors {e M } is generated from a set of coordinates {x M (x)} 
via e M = d^x, where <9 M = The reciprocal frame, denoted by {e^}, satisfies 
e M • e v = 5 The vector derivative V(= d x ) is then defined as 

V = . (2.8) 

More generally, the vector derivative with respect to a is denoted d a . Further 
details concerning geometric algebra and the STA may be found in [86, 88] and a 
summary of useful identities is given in appendix B. 

§2.2.2 Gauge Theory of Gravity 

The language of geometric algebra is used by Lasenby, Doran and Gull [121] 
to develop an alternative theory of gravity based on gauge principles alone. The 
theory is known as Gauge Theory of Gravity (GTG). The key ideas in GTG 
are summarized as follows: 

I. The gravitational effects are described by a pair of gauge fields, defined over 
a flat Minkowski spacetime namely, the position gauge field, h(a) and the 
rotation gauge field, Q(a). 
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(a) The first of these, the position gauge field, h(a) describes a position 
dependent linear function mapping the vector a to vectors. The overbar 
serves to distinguish the linear function, h(a) from its adjoint, h(a), 
where 

h(a) = d b h(b) ■ a . (2.9) 

The gauge-theoretic purpose of h(a) is to ensure covariance of the equa- 
tions under arbitrary local displacements of the matter fields in the 
background spacetime. 

(b) The second gauge field, the rotation gauge field, Q(a), is a position de- 
pendent linear function mapping the vector a to bivectors. The purpose 
of this gauge field is to ensure covariance under local (Lorentz) rotations 
of fields at a point in the background spacetime [121]. A gauge field 
a; (a) is introduced to map vectors to bivectors and is defined as 

u(a) = Vth(a) . (2.10) 

2. Physical predictions are extracted in a gauge invariant manner which ensures 
that the background spacetime plays no dynamical role in the physics. 

3. The two gauge fields are minimally coupled to the Dirac field to enforce in- 
variance under local displacements and both spacetime and phase rotations ' . 
One can formulate the Dirac action in a way that it ensures that the internal 
phase rotations and spacetime rotations assume equivalent forms. 

4. In the process, one can construct a Lagrangian density for gravitational 
fields. The first surprise of adopting such an approach is that the demand of 
the minimally-coupled Dirac equation restricts us to a single action integral. 
There is no freedom except for the inclusion of a cosmo logical constant. 

5. This leads to a set of field equations which are completely independent of 
how we choose to label the position of fields by a vector x. 

The resulting theory based on the above key ideas is a first order theory which 
is consistent with quantum mechanics at the first quantized level. One can de- 
velop additional and alternative insight to the general theory of relativity from the 
GTG. The ideas are applied in the context of black holes, collapsing matter and 
cosmology in the second part of [121]. 



^This requires a vector gauge field for electromagnetism. 
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The covariant derivative T> is assembled from the (flat space) vector derivative 
V and the gravitational gauge fields. The action of D on a general multivector M 
is given by 

VM = d a a ■ VM 

_ (2 11) 

= d a [a ■ h(y)M + u(a) x M] , v ' ; 

where we have introduced the rotational gauge field us (a) defined in (2.10). 

The covariant derivative contains a grade-raising and lowering component, so 
that we may write 

VM = V-M + VAM , (2.12) 

where 

V ■ M = d a -(a-VM) , V A M = d a A(a-VM) . (2.13) 

The single-sided transformation law of spinors under rotations requires us to in- 
troduce a seperate spinor covariant derivative Dip = d a a • Dip, where 

a ■ Dip = a ■ h(V)ip + \-uu(a)ip . (2.14) 

The field strength corresponding to the fl(a) gauge field is defined by 

R(a A b) = a ■ VQ(b) - b ■ Vfi(a) + Q(a) x Q(b) . (2.15) 

From this, we define the covariant Riemann tensor 

U(aAb)=Rh(aAb). (2.16) 

The Ricci tensor, Ricci scalar and Einstein tensor are then given by 

H(a) = d b ■ TZ(b A a) , (Ricci tensor) (2.17) 
TZ = d a - 11(a) , (Ricci scalar) (2.18) 

Q{a) = K(a) - \-oR . (Einstein tensor) (2.19) 

§2.3 The Self-Consistent Dirac Field 
and the Field equations 

The equations of motion for the Dirac field and the gauge fields are derived 
from a form of the minimally coupled action from [49, 121]: 



Sscv = J \d 4 x\ det(h) 1 



1Z ~ 

— — k 2 (Dipi^ip — mipip) — A 



^2.20) 
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Varying Eq. (2.20) with respect to fl(a), h(a) and ip respectively (refer to the 
procedure in [121]), the following equations describe a Dirac field of mass m coupled 
self-consistently to gravity [49] with the inclusion of the cosmological constant in 
this paper. 



K 2 

V A h(a) = k 2 S ■ h(a) + — [db ■ 5(6)] A h(a) , 
Q{a) -Aa = n 2 T(a) , 



covariant spin tensor, and is defined by 

S(a) = Sh^ 1 ^) , 
and S is the spin trivector and is defined by, 

5 = ^73^ , 
and the matter stress-energy tensor is defined to be 

T(a) = (a ■ D^i lz f) . (2.26) 



('wedge') 


(2.21) 


(Einstein) 


(2.22) 


(Dirac) 


(2.23) 


The term S(a) 


is the 




(2.24) 




(2.25) 



The above equations are the GTG analogue of the Einstein-Cartan-Dirac equa- 
tions. The second term on the RHS in Eq. (2.21) is set to zero because the 
minimally-coupled equations for the matter fields are obtained from a minimally- 
coupled Lagrangian only if the contraction of the spin tensor vanishes, i.e. d a • 
S(a) = 0. One could include cases for non- vanishing spin tensor since the spin 
tensor vanishes for scalar fields and Yang-Mills gauge fields. For spin-1/2 fields, 
we can use the form S(a) = S ■ a, where S is the spin trivector. Further results 
for a non- vanishing tensor are given in [49]. 

The wedge equation in the spin-torsion theory is algebraic in a; (a), and may 
be solved for a vanishing spin tensor to give 

uj(a) = -h\y A h-\a)} + -a ■ {d b A h[V A h~\b)}} + \n 2 a ■ S . (2.27) 



§2.4 Massive Non-Ghost Solutions to ECKS the- 
ory 

In this section, we derive a self-consistent solution for a massive Dirac field ip, 
which is both homogeneous and isotropic at the level of classical fields, with the 
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inclusion of a cosmo logical constant. We note that the classical fields do not feel 
the anisotropy of the u;-function which arises due to the spin of the Dirac field, as 
the classical fields couple to gravity via the fa-function only. We first summarize 
the procedure in [35] which we adopt and thn extend the solutions to the case with 
a cosmological constant. 

We begin by introducing a set of polar coordinates, 

x ' 7 

t = x • 70 , cos 9 



r = yj{x A 7o) 2 , tan = 



' 2 (2.28) 



x ■ 7 



x ■ 7 1 



We also make use of the following vectors which are members of the polar coordi- 
nate frame: 

_ x A 7070 fo on x 

e t = 7o , e r = . (2.29) 

r 

We choose the gauge based on the assumed symmetry for the /i-function for a 
vector b, and giving [121]: 

/ kr 2 \ 

h(a) = b ■ e t e t + a(t) 1 + — )b A e t e t , (2.30) 



where a(t) is the scale factor of the Universe, and the curvature constant k = 
— 1,0,-1-1 defined for open, flat and closed Universes respectively. Note that h(a) 
is globally defined for all k. 

The isotropic line element is generated by the /i-function 

/ kr 2 \~ 2 

ds 2 = dt 2 -\a{t)Y 2 (l + ^J [dr 2 + r 2 {d9 2 + sm 2 9d(t) 2 )} . (2.31) 

In this gauge, the fundamental observers have covariant velocity e t where t is the 
cosmic time, and their surfaces of homogeneity are t =constant surfaces. 
From Eq. (2.27), we obtain uj(b) in the following form : 

u(b) = H(t)b A e t - )-kra(t)e r {b A e t e t ) + ]-n 2 b ■ S , (2.32) 

where H(t) is the Hubble parameter in cosmic time, i.e. H(t) = a(t)/a(t) and a 
is the derivative of the scale factor with respect to cosmic time, t. 

We can cast the Dirac equation from Eq. (2.23) in the following form 



kr \ 3 1 3 

e t d t + a(t) ( 1 + — J e t e t A V + -He t + -kra(t)e r + -k 2 S 



ipia 3 = m^7o • 

(2.33) 
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Following Isham and Nelson [96], if the gauge invariant observables (for exam- 
ple, the the projection of the Dirac current onto the velocity of the fundamental 
observers constructed from ip) are to be homogeneous, we must have tp = ip(t). 
Then the only term with any dependence on the spatial coordinates would be 
(k/2)a(t)e r ipi<73 in Eq. (2.33). We thus arrive at the conclusion that the Universe 
must be spatially flat (i.e. k — 0), if we require the observables associated with the 
Dirac field to be homogeneous * [96]. The motivation for the above gauge choice 
is attributed to the fact that h must be globally defined, for all choices to k. 

From now onwards, we specialize to the flat Universe (k = 0). It is then 
convenient to work in a new gauge reached by the displacement, 

x' = f(x) = x ■ e t e t + a(t)x A e t e t , (2.34) 

which turns /i-function to the simple form 

h(b) = b + rH(t)b-e r e t . (2.35) 

This generates the line-element 

ds 2 = (1 - r 2 H 2 )dt 2 + 2Hrdtdr - [dr 2 + r 2 [dQ 2 + sin 2 6d(j) 2 )} . (2.36) 

The above gauge choice is known as the Newtonian gauge, and has an analogue 
in GR which employs the description of a set of geodesic clocks in a radial freefall, 
comoving with the fluid. 

The uj{a) transforms to give 

u(b) = H(t)bAe t + ^-(b-S), (2.37) 

which leaves the homogeneous Dirac field, ip unchanged. The surfaces of homo- 
geneity still have t =constant. The covariant velocity of the fundamental observers 
remains e t , but their radial coordinates are proportional to the scale factor a(t). 

We now consider the spin-torsion Universe with a non-zero cosmological con- 
stant A. The Riemann tensor evaluates to be 

11(B) = —HB ■ etet - H 2 B + \^B ■ SS - \k 2 (B ■ V) ■ S , (2.38) 

for any arbitrary bivector B, and the Einstein tensor Q{b) is found to be 

Q(b) = 2Hb A etet + 3H 2 b - ^n 2 b ■ (V ■ S) + l -n\b ■ S)S - ^S 2 b , (2.39) 

*One can arrive at the same conclusion if one attempts to solve the Dirac equation non-self- 
consistcntly on a homogeneous gravitational background. The reason is due to the fact that the 
only change to the equation (2.33) would be that S is now the torsion trivector of the background 
rather than the spin of the Dirac field. 
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The stress energy tensor of the self-consistent Dirac field evaluates to 

T(a) = (b-e t ij>in$ + HbAe t S + ^(b-S)S^ , (2.40) 

where we assumed that tp is a function of time alone. The Dirac equation reduces 
to 

e t d t + ^He t + ^k 2 s) ijia 3 = m^ 7o , (2.41) 
from which we deduce 

V '° ' (2.42) 

V-S = -2He t ■ S . 
Substituting Eq. (2.42) into Eq. (2.22), we obtain the following equation : 

2Hb A e t e t + 3H 2 b + ^K, 4 S 2 e t be t - mK 2 b ■ e t e t (ip^ - Ab = , (2.43) 

from which we extract the following set of scalar equations 

2H + 3H 2 - ^k a S 2 - A = , (2.44) 

3H 2 + ^k 4 S 2 - mn 2 (^fj - A = . (2.45) 

If we exclude torsion from the model, from Eq. (2.44), the equation becomes 

3H 2 + 2H - A = . (2.46) 

Since S 2 < 0, the effect of torsion would make H more negative for any given value 
of the Hubble parameter H . Consequently, torsion enhances singularity formation 
for a Hubble parameter at a particular instant. It follows that the Universe will 
be younger if we look at the singularity theorems in the presence of torsion [102]. 
Torsion effects become significant when the Compton wavelength of the field is 
larger than the Hubble radius. At earlier times, where p is singular, the torsion 
effects are dominant. However, the inclusion of a positive cosmological constant 
would provide a repulsive force which would dominate the late time behaviour 
of the Universe. The competing effects between the torsion and the cosmological 
constant lead to interesting consequences in the models which we are about to 
examine. 

Suppose one sets m — 0, then Eqs. (2.44) and (2.45) would be solved immedi- 
ately to give 

tf(t) = y|coth^. (2.47) 
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The cosmological constant alters the large t behaviour dramatically, as H would 
tend to a finite constant, i.e. H — > y^f- 

We consider the case where ipip ^ 0, and parametrise ip by the following relation 
[86]: 

i> = v/pexp (Ijpj R , (2.48) 

where p > and (3 are scalar functions of t. The rotor R is even-grade (and time 
dependent) and satisfies RR = 1. Substituting Eq. (2.48) into Eq. (2.41), and 
equating grades on either side, one finds the following equations 

p = Am sin /5J 7o AS- 3pH , (2.49) 
p/3 = A{m cos + 3vrp) J70 A S , (2.50) 
pRR = -2(me~ J/3 + 37rp) 7o • S . (2.51) 

where RR is a bivector. We simplify the calculation by taking sin (3 = 0. Since (3 
is then constant, a non-zero and constant p is forbidden by Eq. (2.49), since Eq. 
(2.50) gives 70 A S = 0. It follows that we must solve 

p = -3pH , (2.52) 

R = -{m cos (3 + 37rp)7 0j RJ73 , (2.53) 

subject to the constraint 

3H 2 - Ylix 2 p 2 - 8nmp cos (3 - A = . (2.54) 

In these equations, cos/5 = ±1. The solutions with cos/5 = 1 are regarded as 
"particle" (positive energy) solutions and those with cos/5 = — 1 as "anti-particle" 
(negative energy) solutions in the absence of gravity [48] . 

To solve the Eqs. (2.52) and (2.54), we make use of the substitution y — ^ 
to find 

y 2 = 1 ± Amy + 3Ay 2 

nk ( 2m V4:m 2 -3A\ ( , 2m VAm 2 -3A\ (2.55) 
= 3A y ± — — + — ) [y ± — 



3A 3A J \ 3A 3A 



One could obtain, by solving the above equation, the density parameter p and the 
Hubble parameter, H(t), for all t in both particle and anti-particle sectors. 
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§2.4.1 Cosmological implications 

In this section, we analyze how the inclusion of a cosmological constant would 
modify the solutions found in [35] for a spin-torsion Universe (see table 2.1 for the 
classification of the solutions). Before we proceed to the solutions, we examine the 
discriminant \/ Am 2 — 3A and the limits which these solutions can take. There are 
three possible cases: (i) Am 2 3> 3A, (ii) Am 2 = 3A and (iii) Am 2 <C 3A. In natural 
units, the cosmological constant, A ~ 1CT 122 . We compare this value to the known 
experimental values from the masses of the fermions in natural units. The mass of 
the electron, m e ~ 1CT 23 and the mass of the neutrino is estimated approximately, 
m v ~ lCT 7 m e (see the experimental bounds for the neutrino masses in [69]). It is 
found that the values of m 2 for the fermions are significantly much larger than the 
cosmological constant, hence, the preferred model would be case (i). For case (ii), 
m ~ 10~ 38 m e . However no such fermion is known to exist in nature, and hence 
case (ii) is not a plausible model. Since case (ii) is not a plausible model, case (iii) 
is effectively ruled out as well. In this section, we consider the solutions of interest 
for both case (i) and (ii). 

To solve Eq. (2.55), we demand the condition y ^ since p > 0. We note 
that these solutions would work for all values of t. We start with case (i) with 
4m 2 ^> 3A, and assume a positive cosmological constant i.e. a repulsive force. 
From Fig. 2.1, since p > 0, the requirement imposes an asymmetry between the 
particle and anti-particle solutions. For the particle case, we consider only one 
branch of the solution with y > [region (i)]. On the other hand, we observe 
that the anti-particle case has two branches of solutions. These two branches are 
labelled as regions (ii) and (iii). For A — > 0, the solutions found in [35] would 
correspond to regions (i) and regions (ii), with the appropriate limit taken. The 
presence of the cosmological constant has introduced an additional solution. 

For regions (i) and (ii), the Hubble parameter and density parameter for both 
particle and anti-particle respectively are: 

a 2 

^ ^ 67r sinh crf(a cosh erf ± m sinh erf) ' ^ ^ 

a(a cosh 2erf ± m sinh 2at) 

Hit) = 1 7 : r 2 -^ , 2.57 

3 sinh at(a cosh at ±m smh at) 

where a = y -^-- The density is singular at t — and at t — a" 1 coth~ 1 [±(m/a)] 
for both particle and anti-particle cases respectively. 

Starting from the particle case for t > 0, the solution in region (i) corresponds 
to an expanding Universe from an initial singularity at t — 0, and H(t) tends to the 
value for a flat Universe with a cosmological constant at later times. During this 
epoch, the scale factor of the Universe asymptotically approaches exp(t^A/3). If 
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Figure 2.1: A sketch of the equation (2.55) with y 2 vs y for both particle 
and anti-particle respectively, i.e. ±m. This is for the case, 4m 2 > 3A. The 
three regions of interest : (i) an expanding Universe with particle from an initial 
singularity (for t > 0) or a Universe contracting to a singularity in some finite 
time (for t < 0), (ii) an anti-particle Universe which expands to a finite value 
and contracts to a singularity (for t > 0), and (iii) an expanding anti-particle 
Universe (for t > 0) and a contracting Universe which bounces at t = and 
expands indefinitely (for t < 0). 



we move from £ — > — £, we found p(t) > 0, but H(t) < 0. For £ < 0, this solution 
describes a Universe contracting to a singularity in some finite time. 

Next, we examine the anti-particle solution in region (ii). For £ > 0, the 
Universe would expand from an initial singularity, turn around at some finite time 
and collapse to a singularity. A particle horizon is also present in the model which 
continues to exist right up to the singularity at the endpoint of the collapse. For 
Eqs. (2.56) and (2.57), taking the limit A — > 0, the solutions would be those found 
in [35]. For £ — * — £, the solution describes the same cosmology as before. 

From Fig. 2.1, the quadratic nature of Eq. (2.55) introduces a new anti-particle 
solution [region (iii)]. The Hubble parameter and density parameter for region (iii) 
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are found to be: 

a 2 1 



*<«> = T 



37T m + M cosh 2at ' 
sinh 2at 



M + cosh 2at 



(2.58) 
(2.59) 



where M = [mj \Jm 2 — a 2 ). This solution satisfies Eq. (2.55) for the anti-particle 
case and is geodesically complete. For t > 0, the solution describes an expanding 
Universe starting from some finite y. With the inclusion of a positive cosmological 
constant, the Universe generated by the anti-particle sector appears to be an eter- 
nal Universe with no singularities. The t < is an interesting case, since p(—t) > 0. 
This solution describes a contracting Universe from some finite —t, bounces at the 
instant at y = and expands indefinitely like the case for t > 0. The solution 
has the property that it smoothly matches a contracting phase (t < 0) onto an 
expanding phase (t > 0) which is reminiscent of the pre-Big-Bang scenario [72]. 

In principle, we can incorporate a negative cosmological constant by switching 
A — > —A in the above solutions. For both particle and anti-particle sectors, the 
solutions are similar to a closed FRW model, which corresponds to an expanding 
Universe starting from an initial singularity which turns around at a finite time, 
and collapses back into a singularity. 

Another solution of interest is the extremal case where 4m 2 = 3A. From Fig. 
2.2, we note that there are three regions of interest, one solution for the particle 
case and two solutions for the anti-particle case. The solution is degenerate at 
p = mj (37r). 

By setting the initial conditions when t = 0, and the initial density, p 
we obtain the solutions of (2.55) for 4m 2 = 3A. Starting from the particle case, 
for region (i), we obtain the density and Hubble parameter to be: 

Tfl 

P = 37r(e 2m* _ i) (P article ); (2-60) 

H(t) = 3(1 - e-^t) (P article )- (2-61) 

From the above solutions, as t — > oo, the particle Universe, H{t) would tend to a 
finite constant \J (A/3), which is typical of a Friedmann dust model with a positive 
cosmological constant. 

For region (ii), the solutions are: 

P = 37r(1 ™ e -2 mt) (anti-particle), (2.62) 
H ® = 3 ( e 2mt _ I) (anti-particle). (2.63) 
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Figure 2.2: A sketch of the equation (2.55) with y 2 vs y for both particle and 
anti-particle respectively, i.e. ±m. This is for the extremal case: 4m 2 = 3A. 
Like 2.1, we have three regions of interest : (i) an expanding Universe with 
particle from an initial singularity (for t > 0) or a Universe contracting to a 
singularity in some finite time (for t < 0), (ii) an anti-particle Universe which 
expands to a finite value and contracts to a singularity (for t > 0), and (iii) 
an expanding anti-particle Universe (for t > 0) and a static Universe which 
smoothly matches from —t to t (for t < 0). 



For the anti-particle Universe, H{i) would tend to zero with the same limit, which 
suggests the scale factor of the anti-particle Universe tends asymptotically to a 
constant (see Fig. 2.3). In this scenario, the Universe is a Minkowski space with 
torsion at late times. The sole parameter of the extremal Universe is dependent 
on the mass of the fermion, m. Both models represent an expanding Universe 
with an initial singular density which decays as it evolves in time (except that 
the anti-particle decays to m/3ii). In both cases, as we switch t — > —t, we get 
collapsing particle and anti-particle solutions which contract to a singularity. 



In addition, there is an additional anti-particle solution [region (iii)] for p > 
m/(3n) in the extremal case. The density parameter and the Hubble parameter 
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are found to be: 

Tfl 

H(t) = 30TF=) < 2 ' 65) 

This solution is similar to the region (iii) for the case Am 2 3> 3A. For t > 0, it is 
an expanding anti-particle Universe. The t < solution appears to be a Universe 
which smoothly matches a contracting phase onto an expanding phase (t > 0) at 
a finite density, with a bounce at p = m/(67r). However, for the t < solution, 
since the anti-particle solution starts at t = — oo, indicating that if(t) — > 0. This 
means that the Universe would not be expanding, and this scenario is similar to 
the Einstein static Universe, but in this case, the spin-torsion has balanced against 
the repulsive force of the cosmological constant. 
The stress-energy tensor T(a) is found to be 

T(e t ) = p(2vrp ± m)e t , 
7X73) = , 

r( 7 i) = +\ p h 12 - 7T P 2 7l , ( 2 - 66 ) 

7(72) = --pi?7i - 7rp 2 7 2 . 

As expected, e< is a timelike eigenvector of T(a). The only other real eigenvector 
is 73 which is dual to the spin trivector S. (Taking the dual is performed by multi- 
plying the pseudoscalar X.) The stress-energy tensor singles out the two directions 
in space to be algebraically special, and reflects the anisotropy of the solution at 
the level of the gravitational gauge fields. The energy density as measured by the 
fundamental observers is evaluated to be 

T(e t ) -e t = p{2irp±m) . (2.67) 

From the energy density, we evaluate the density parameter which is defined 
by the following relation 

n= s?lM^ t (268) 

which in regions (i) and (ii) evaluates to 

4 (a 2 + 3m 2 sinh 2 at ± 3ma sinh at cosh at) (2 69) 

3 (acosh2a:t ± msinh2at) 2 



and in region (iii) 



3 



3m 2 + 2a 2 + 3mb cosh 2at 
b 2 sinh 2at 



[2.70) 
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Figure 2.3: Graph of the scale factor a(r) vs conformal time (in units of Hqt) 
for the particle and anti-particle solutions for the extremal Universes at t > 
i.e. Am 2 = 3A for S7a = 0.7. The scale factor for flat Universe (with and 
without the cosmological constant) are plotted for comparison. Note that the 
scale factor for an extremal particle solution is similar to a flat universe with 
positive A, whereas the scale factor for an extremal anti-particle asymptotes 
towards a constant. 



where b = \/m 2 — a 2 . Setting A — > 0, we obtain the result in [35]. Hence at early 
times, Q ^ 1 since 4m 2 ^> 3A. 
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§2.4.2 Massive Non-Ghost Solution 

Finally, we arrive at the final form of the non-ghost solution. To solve the rotor 
equation (2.53), we would note the following result from [35]: 

7oA5 = =>• 7s • (R-foR) = . (2.71) 

The rotor R can be decomposed into a product of rotors R = $L, where $ and 
L are uniquely determined by the requirements $7 = 7o$ and L70 = 7o-£/. The 
condition in Eq. (2.71) restricts L to the form 

L = e ^( sin ^ <Tl+cos, ? <T2 ) (2 72) 

where £ and i] are scalars. However from the rotor Eq. (2.53), we find that the 
vector R'-foR is time independent. Since L is determined uniquely by this vector 
and 70, therefore L must be constant. Now we substiture R = $L into Eq. (2.53), 
and note that L commutes with Z73, to find that 

$ = -x$l<J3 , (2.73) 
where the scalar x(t) is defined by 

X = m cos f3 + 3np . (2.74) 
Hence Eq. (2.73) has the solution 

$(t) = $ oe - XCT «W ? (2.75) 

where $0 is a constant rotor which commutes with 70. The rotor may be eliminated 
by combining the global position-gauge transformation x' = $0^0, with the global 
rotation-gauge transformation defined by the rotor $ . The solution to Eq. (2.74) 
is 

X(t) = ±mt + 3tt J p(t)dt , (2.76) 

where we could substitute the various forms of p(t) found in Eqs. (2.56), (2.58), 
(2.60) or (2.62). The integration constant is chosen so that x is defined for p > 0. 
For regions (i) and (ii), we found that x is evaluated to be 



X(t) = ±mt + ^ In 



a 



1 +tanh 2 f ) ±2mtanhf 
tanhf 



(2.77) 



Similarly, for region (iii), which only applies to the anti-particle case, 



(2.78) 
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For the extremal Universes, for both particle and anti-particle cases we find x 
to be 

X(t) = ±mt + ]- m(l - e- 2mt ) [particle, region (i)] , (2.79) 
X 

X(t) = ±mt + - ln(l - e 2mt ) [anti-particle, region (ii)] , (2.80) 

X(t) = ±mt - - ln(l + e' 2rnt ) [anti-particle, region (iii)] . (2.81) 
2 

Hence we obtain the general self-consistent massive non-ghost solution for a 
Dirac field with the inclusion of a cosmological constant, ip with sin/5 = 0, with 
(3 = or 7r to be 



V>(*) = ex P (^Y ~ I( ?3X(t)^ exp Q(sin?7(7i + cos^)^) • (2. 



82) 



i2.5 Discussion 



Finally, we examine whether the solutions for the spin-torsion Universes are 
inflationary or not. From (1.4) in Chapter 1, the condition for inflation is 

d>0 ^H + H 2 >0. (2.83) 

Using (2.54), we obtain 

H + H 2 = ^ + g(p,0), (2.84) 



where 



g(p,(3) = -8n 2 p 2 T ^, (2.85) 

for the cases of particle and anti-particle respectively. For the particle case, 
g(p,P) < 0, and hence there is no inflation at earlier times since g(p, (3) ^> A. 
However, at late times, there will be de-Sitter inflation induced by the cosmolog- 
ical constant. On the other hand, it is observed for the anti-particle case that 
g(p, f3) > if p < {itl/Q-k). For p < (m/67r), the spin-torsion Universe would have 
an inflationary phase driven by the anti-particle. From Eq. (2.94) in Table 2.1, the 
lower bound for the density of a spin-torsion Universe is p = (to/37t). Therefore, 
the spin-torsion solutions do not incorporate an inflationary phase at earlier times. 
The only inflationary phase is provided by the cosmological constant at late times. 
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A Table of the cosmological solutions in Spin Tor- 
sion Theories 



Solution Type 


Hubble Parameter, H(t) 


Flat: 


(i) (K = 0,A = 0) 

H = | . (2.86) 

(ii) with A > 

2a 

H = —cothat. (2.87) 

(iii) with A < 

H=^-cotat, (2.88) 

where 

3A 

a = — . (2.89) 


Spin Torsion Universe: 




(a) Particle 


General solution [region (i)]: 

a(a cosh 2at + m sinh 2at) 

Hit) = . (2.90) 

3 sinh at (a cosh at + m smh at) 

(i) Setting 4m 2 > 3A, we obtain Eq. (2.87) for A > 
and Eq. (2.88) for A < 0. If we set A -> 0, we get the 
solution in [35]: 

H(t)= \ +2m \. (2.91) 
v ; 3t(l + mt) K ' 

(ii) For Am 2 = 3A, we get 

2m 

2.92 

3[1 - exp(-2mf)] V ; 
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(b) Anti-particle 



General solution [region (ii)] 

a(a cosh 2at — m sinh 2at) 



H(t) 



3 sinh at(a cosh at — m sinh at) 



(2.93) 



(i) Setting 4m 2 > 3A, we obtain (2.87) for A > and 
the inverse of (2.88) for A < 0. If we set A — > 0, we get 
the other solution in [35]: 



H(t) 



1 - 2mt 
3t(l - mt) 



(ii) For 4m 2 = 3A, we get 



2m 



3[exp(2mt) 



General solution II (Region (hi)) 
(i) Setting 4m 2 ^> 3A, we obtain 



*<«> = f 



sinh 2at 
M + cosh 2at 



where 



M 



m 



V: 



a" 



(ii) For the extremal case, we have 

Hit) = 

v ; 3[exp(2mt) + 1] 

There is no analogue for this case as A — > 0. 



(2.94) 



(2.95) 



(2.96) 
(2.97) 

(2.98) 



Table 2.1: A summary of the solutions for the flat Uni- 
verses and the spin-torsion Universes 



Part II 
Braneworld Cosmology 



Chapter 3 



Braneworld 
Formalism 



Cosmology I: 



"I do not wish to hide the fact that I can only look with repugnance 
.. upon the puffed up pretentiousness of all these volumes filled with 
wisdom, such as are fashionable nowadays. For I am fully satisfied 
that .. the accepted methods must endlessly increase these follies and 
blunders, and that even the complete annihilation of all these fanciful 
achievements could not possibly be as harmful as this fictitious science 
with its accursed fertility. " - Immanuel Kant 

In this chapter, we discuss the dynamics of linearized scalar and tensor per- 
turbations in an almost Friedmann-Robertson-Walker braneworld cosmology of 
Randall- Sundrum type II using the 1+3 covariant approach based on [123] and 
[124]- We derive a complete set of frame-independent equations for the total mat- 
ter variables, and a partial set of equations for the non-local variables which arise 
from the projection of the Weyl tensor in the bulk. The latter equations are incom- 
plete since there is no propagation equation for the non-local anisotropic stress. 
As a first step, we derive the covariant form of the line of sight solution for the 
CMB temperature anisotropics in braneworld cosmologies. We discuss possible 
mechanisms by which braneworld effects may remain in the low- energy universe in 
Chapter 4- 

§3.1 Introduction to Braneworld Cosmology 

It is understood that Einstein's theory of general relativity is an effective the- 
ory in the low-energy limit of a more general theory. Recent developments in 
theoretical physics, particularly in string theory or M-theory, have led to the idea 
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that gravity is a higher- dimensional theory which would become effectively four- 
dimensional at lower energies. 

Braneworlds, which were inspired by string and M-theory, provide simple, yet 
plausible, models of how the extra dimensions might affect the four-dimensional 
world we inhabit. There is the exciting possibility that these extra dimensions 
might reveal themselves through specific cosmological signatures that survive the 
transition to the low-energy universe. It has been suggested that in the context 
of braneworld models the fields that govern the basic interactions in the standard 
model of particle physics are confined to a 3-brane, while the gravitational field 
can propagate in 3 + d dimensions (the bulk). It is not necessarily true that 
the extra dimensions are required to be small or even compact. It was shown 
recently by Randall and Sundrum [160] for the case d — 1, that gravity could 
be localised to a single 3-brane even when the fifth dimension was infinite. As a 
result, the Newtonian potential is recovered on large scales, but with a leading- 
order correction on small scales: 

where the 5-dimensional cosmological constant A oc — /~ 2 . As a result, general 
relativity is recovered in 4 dimensions in the static weak-field limit, with a first- 
order correction which is believed to be constrained by sub-millimeter experiments 
at the TeV level [141, 160]. 

The aim of this chapter is to set up the evolution and constraint equations for 
perturbations in a cold dark matter (CDM) brane cosmology, presenting them in 
such a way that they can be readily compared with the standard four-dimensional 
results, and to provide approximate solutions in the high and low-energy universe 
under certain restrictions on how the bulk reacts on the brane. Our equations 
are clearly incomplete since they lack a propagation equation for the non-local 
anisotropic stress that arises from projecting the bulk Weyl tensor onto the brane, 
and our solutions are only valid under the neglect of this stress. However, our pre- 
sentation is such that we can easily include effective four-dimensional propagation 
equations for the non-local stress should such equations arise from a study of the 
full bulk perturbations. The lack of a four-dimensional propagation equation for 
the non-local stress means that it is currently not possible to obtain general results 
for the anisotropy of the CMB in braneworld models. Such a calculation would re- 
quire solving the full five-dimensional perturbation equations which is non-trivial 
since the equations can only be reduced to two-dimensional partial differential 
equations on Fourier transforming the 3-space dependence. 

The theory of cosmological perturbations using the metric-based approach, or 
the Bardeen approach has been extensively studied in [12, 20, 44, 45, 50, 100, 
101, 107, 108, 110, 111, 116, 147, 148, 149, 165, 167, 178, 186, 188] (with Z 2 
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symmetry) and in [162] (without Z 2 symmetry). On the other hand, it has been 
also explored by various authors within the 1+3 covariant methods in [21, 77, 123, 
124, 135, 136]. We will further summarise some of the qualitative implications for 
observational cosmology of the braneworld scenarios in the literature, and explore 
the possible implications for the CMB anisotropies in braneworld cosmology in the 
next chapter. 

We give a brief summary of the (l+3)-covariant approach in §3.2. In §3.3, we 
look at the field equations of the braneworld cosmology. Following that, we will 
look at the linearized perturbation equations in §3.4, where the local and nonlocal 
conservation equations are given in §3.4.1, and the propagation, divergence and 
constraint equations are given in §3.4.2. Following that, we discuss the cosmolog- 
ical perturbations on the brane in §3.5, distinguishing between scalar perturba- 
tions (§3.5.1), vector perturbations (§3.5.2) and tensor perturbations (§3.5.3). To 
conclude this chapter, we will derive a covariant expression for the temperature 
anisotropy in §3.6. 

§3.2 (1+3)- Covariant Approach in Cosmological 
Perturbation Theory 

The (l+3)-covariant approach, pioneered by Ehlers [57], Ellis [59] and Hawk- 
ing [79] provided an alternative gauge-invariant treatment of cosmological per- 
tubations [61, 65] other than the standard metric-based approach formulated by 
Bardeen[8]*. In the (1+3) treatment, based on the covariant formulation of hydro- 
dynamics and gravito-dynamics of Ehlers and Ellis (for an extensive review of the 
subject, see [11, 22, 57, 59, 60, 62]), the cosmological perturbations are described 
by gauge-invariant variables which have simple physical interpretations in terms 
of inhomogeneous measures of anisotropy of the Universe. 

The covariant approach has been applied in many areas of cosmology, for 
example, gravitational waves (see e.g. [54]), inflation [23, 53], magnetic fields 
[181, 182, 183] and radiative transport theory [63, 64]. This develop a strong case 
for its application in the study of cosmic microwave background physics. The line 
of sight calculation of CMB anisotropies under the instantaneous recombination 
based on the covariant approach was obtained by Dunsby [52] and Challinor and 
Lasenby [33] based on earlier work in [51, 61, 63, 64, 65, 180]. Subsequently, 
the covariant approach was extended to studying CMB anisotropies in the CDM 

*The Bardeen approach is inherently linear. It assumes a nonlocal decomposition of the 
perturbations into scalar, vector and tensor modes at the outset, with each of them treated 
seperately. Even though the Bardeen variables are gauge- invariant, they are not physically 
transparent, i.e. they do not characterise perturbations in a way that is amenable to simple 
physical interpretation. 
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model in [34], tensor anisotropies induced by gravitational waves in [31], polar- 
ization in [29, 32], the mode and multipole representations and their application 
to temperature anisotropies in [73, 74] and obtaining model-independent limits on 
the inhomogeneity and anisotropy of the CMB anisotropy on large scales in [139]. 

The covariant approach has the advantage that it does not employ the non- 
local decomposition into scalar, vector and tensor modes. If required, that can 
be adopted at later stages of the calculation in aid to solve the equations. Even 
if one denies that working with gauge- invariant quantities is a significant merit, 
a key advantage of the covariant approach is that one is able to work exclusively 
with physically relevant quantities, satisfying equations which make manifest their 
physical consequences. Another advantage of the covariant approach is that the 
full extension to non-linear methods [140] is straightforward. 

We start by choosing a 4- velocity u a . This must be physically defined in such 
a way that if the universe is exactly FRW, the velocity reduces to that of the 
fundamental observers to ensure gauge-invariance of the approach. From the 4- 
velocity u a , we construct a projection tensor h ab into the space perpendicular to 
u a (the instantaneous rest space of observers whose 4- velocity is u a ): 

Kb = 9ab + u a u b , (3.2) 

where g a b is the metric of spacetime. The operation of projecting a tensor fully 
with h ab , symmetrizing, and removing the trace on every index (to return the 
pro jected-symmetric-tr ace-free or PSTF part) is denoted by angle brackets, i.e. 

T(ab...c)- 

The symmetric tensor h ab is used to define a projected (spatial) covariant 
derviative D a which when acting on a tensor T b '" c d... e returns a tensor that is 
orthogonal to u a on every index, 

D a T b - c d ... e = h%h\ . . . h c s h l d . . . h u e V v T r -\.. u , (3.3) 

where V a denotes the usual covariant derivative. 

The covariant derivative of the 4-velocity can be decomposed as follows: 

V a u b = u ab + a ab + ^Qh ab - u a A b , (3.4) 

where u ab is the vorticity tensor which satisfies u a uj ab = 0, o ab = a^ ab ) is the shear 
which is PSTF, ee V a w a = 3H measures the volume expansion rate (where H 
is the local Hubble parameter), and A a = u b V b u a is the acceleration. In addition, 
we define the vorticity vector, oj a = — (l/2)curlu . 

Gauge-invariant quantities can be constructed from scalar variables by taking 
their projected gradients. Such quantities vanish in the FRW limit by construction. 
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The comoving fractional projected gradient of the density field p^ of a species i 
(for example, photons) is one important example of this construction: 

AW ^ jsDj* , (3.5) 

where a is a locally defined scale factor satisfying 

a = u b V b a = Ha , (3.6) 

which is included to remove the effect of the expansion on the projected gradi- 
ents. Another important vector variable is the comoving projected gradient of the 
expansion, 

Z a = aD a Q , (3.7) 

which provides a measure of the inhomogeneity of the expansion. 

The matter stress-energy tensor T ab can be decomposed irreducibly with respect 
to u a as follows: 

Tab = pu a u b + 2u( a q h) + Ph ab + TT ab , (3.8) 

where p = T ab u a u b is the energy density measured by an observer moving with 
4-velocity u a , q a = —h b a T bc u c is the energy flux or momentum density (orthogonal 
to u a ), P = h ab T ab /3 is the isotropic pressure, and the PSTF tensor rr ab = is 
the anisotropic stress. 

The remaining gauge-invariant variables are formed from the Weyl tensor C abc d 
which vanishes in an exact FRW universe because these models are conformally 
flat. The ten degrees of freedom in the 4-dimensional Weyl tensor can be encoded 
in two PSTF tensors: the electric and magnetic parts defined respectively as 

E ab = C abcd u b u d , (3.9) 
H ab = ^C acst u c rj st bd u d , (3.10) 

where r) abc d is the 4-dimensional covariant permutation tensor. 

In the covariant approach, we defer from making the frame choice until all 
the relevant equations, valid for u a are derived. There are many different frame 
choices and we state the two main examples''' which we employ in this thesis: 
(i)the pressure-free CDM frame, where the rest frame of the CDM defines a 
geodesic frame i.e. the peculiar velocity of the CDM is zero, and it defines a 
zero-acceleration frame, A a = 0; and (ii) the energy frame, which is defined by 
vanishing total energy flux q a = 0. 



t Other frame choices are discussed in a review paper by Gebbie, Dunsby and Ellis [74]. 
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3.3 Field Equations of Braneworld Cosmology 



In a recent paper, Maartens [135] introduced a formalism for describing the 
non-linear, intrinsic dynamics of the brane in Randall-Sundrum type II braneworld 
models in the form of bulk corrections to the 1+3 covariant propagation and 
constraint equations of general relativity. This approach is well suited to iden- 
tifying the geometric and physical properties which determine homogeneity and 
anisotropy on the brane, and serves as a basis for developing a gauge-invariant 
description of cosmological perturbations in these models. 

An important distinction between braneworlds and general relativity is that the 
set of 1+3 dynamical equations does not close on the brane. This is because there is 
no propagation equation for the non-local effective anistropic stress that arises from 
projecting the bulk Weyl tensor onto the brane. The physical implication is that 
the initial value problem cannot be solved by brane-bound observers. The non-local 
Weyl variables enter crucially into the dynamics (for example, the Raychaudhuri 
equation) of the intrinsic geometry of the brane. Consequently, the existence 
of these non-local effects leads to the violation of several important results in 
theoretical cosmology, such as the connection between isotropy of the CMB and 
the Robertson- Walker geometry. 

The field equations induced on the brane are derived by Shiromizu et al [176] 
using the Gauss-Codazzi equations, together with the Israel junction conditions 
and Z 2 symmetry. We will discuss this approach in detail in the later part of this 
section. In their elegant geometrical approach which they employed to analyze 
the dynamics of the RSII model, the bulk is a 1+4-dimensional spacetime with a 
non-compact extra spatial dimension. What prevents gravity from 'leaking' into 
the extra dimension at low energies is the negative bulk cosmological constant 
A 5 = — Q/£ 2 , where £ is a curvature scale of the bulk. 

In the weak-field static limit, null results in tests for deviations from Newton's 
law impose the limit i < 1 mm. The negative A 5 is offset by the positive brane 
tension A, which defines the energy scale dividing low from high energies. The 
limit I < 1 mm implies A > (100 GeV) 4 , and the energy scales i.e. the effective 
cosmological constant on the brane are related by the following relations: 



where k 2 = 87r/Mf , with M 5 being the fundamental 5-dimensional Planck mass, 
A 5 is the cosmological constant in the bulk and A is the tension of the brane. 
The constant k 2 = 8nG = 8vr/M 4 2 , and M 4 ~ 10 19 GeV is the effective Planck 
scale on the brane. A further intriguing feature of the braneworld scenario is that, 



A 




(3-11) 



A= ^(A 5 + /t 2 A) , 



(3.12) 
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because of the large extra dimensions, the fundamental energy scale of gravity 
can be dramatically lower than the effective Planck scale on the brane - as low 
as ~ TeV in some scenarios [141]. In generalised RSII models, the lower bound 
arising from current tests for deviations from Newton's law is 

M 5 > 10 5 TeV, A > (lOOGeV) 4 (3.13) 

and is related to M 4 via Mf = M\ji. 

A less stringent constraint given by nucleosynthesis [18, 141] implies that A > 
(lMev) 4 which gives 

/l MeV\ s 

M 5 >^^J , M 5 >10TeV. (3.14) 

At energies well above the brane tension A, gravity becomes 5-dimensional and 
significant corrections to general relativity occur. There are also corrections that 
can operate at low energies, mediated by bulk graviton or Kaluza-Klein (KK) 
modes. 

From the above discussion, we proceed to summarise the approach in [176] to 
derive the effective Einstein equations induced on the brane. The induced metric 
on the brane on all the hypersurfaces orthogonal to the unit normal n A (A spacelike 
vector), is defined to be 

9ab = 9ab - n A n B (3.15) 

where we use tildes to denote the 5-dimensional generalization of standard general 
relativity quantities. We adopt the choice of coordinates to be x A = (x a , y), where 
x a = (t, x l ) are spacetime coordinates on the brane and ha = 5 v a- 
Starting from the 5-dimensional Einstein equations, 

G AB = £ 2 [-A 5 g AB + S(y)(-\g AB + T AB )\ , (3.16) 

where the fields are confined to the brane that make up the brane energy momen- 
tum tensor, with Tab^ 3 = 0. 

Applying the Gauss-Codazzi relations to (3.16), 

Rabcd = Refgh9a E 9b F 9c G 9d H + ^Ka[cKd\b , (3-17) 
R B c9A B n c = V B K B A - V A K B B , (3.18) 

we obtain the following equation, 

Gab = ~y A s9ab + K C C K AB - K A C K CB + ^[K CD K CD - (K c c ) 2 ]g A B ~ Sab , 

(3.19) 
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where 

Sab = C ACBD n c n D , (3.20) 

is the projection of the bulk Weyl tensor orthogonal to n A , with £\ab] = = £a A - 
We evaluate (3.19) n the brane as y — > ±0 to get the standard Einstein equation 
on the brane. To do that, Kab has to be determined at the brane. The Israel 
junction conditions across the brane imply that qab is continuous, while Kab 
undergoes a jump due to the energy- momentum on the brane: 



-Kab K AB — k 2 



Tab + \9ab{\-T c c ) 



(3.21) 



We impose Z 2 symmetry* which implies that 

Kab = -K\b , (3-22) 



and obtain 



K 2 



Kab = ~ — 



Tab + \g AB {\ - T c c ) 



(3.23) 



where we dropped the (+) and evaluate the quantities on the brane by taking the 
limit y — > +0. 

Finally, we arrive at the result in [176]. The standard Einstein equation which 
is modified with new terms carrying the bulk effects on the brane: 

G ab = -Ag a b + K 2 T ab + k 4 S ab - £ ab , (3.24) 

where k 2 = 8ir/M%. The bulk corrections to the Einstein equations on the brane 
are made up of two parts: (i) the matter fields which contribute local quadratic 
energy- momentum corrections via the symmetric tensor S ab ; and (ii) the non-local 
effects from the free gravitational field in the bulk transmitted by the (symmetric) 
projection £ ab of the bulk Weyl tensor. The matter corrections are given by 

S ab = ^T c c T ab - -T ac T\ + ^ 9ab [3T cd T cd - (T c c ) 2 ] . (3.25) 

We note that the local part of the bulk gravitational field is the five dimensional 
Einstein tensor Gab, which is determined by the bulk field equations. Conse- 
quently, £ ab transmits non-local gravitational degrees of freedom from the bulk to 
the brane that includes both tidal (or Coulomb), gravito- magnetic, and transverse 
traceless (gravitational wave) effects. 

•'•The choice of Z% symmetry is motivated by a result in string theory by Horava and Witten 
in [91]. They showed that the 10-dimensional Eg x Eg, heterotic string is related to an 11 
dimensional theory on the orbifold R 10 x S 1 jZ^- In that model, the standard model is confined 
on the 4-dimensional spacetime while the gravitons propagate in the full spacetime. 
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The bulk corrections can all be consolidated into an effective total energy den- 
sity, pressure, anisotropic stress and energy flux. The modified Einstein equations 
take the standard Einstein form with a re-defined energy-momentum tensor: 



G ab = -Ag ab + K 2 T%\ 



(3.26) 



where 



rptot 

1 ab 



k 4 1 

Tab H — 2^ ab 2^ ab ■ (3.27) 



K K 

Decomposing £ ab irreducibly with respect to u a by analogy with Eq. (3.8), in 
[77, 135, 136] 



ah 



Uu a u b + 2u( a Q b ) + ^h ab + Vab ) ■ 



(3.28) 



(the prefactor is included to make e.g. U have dimensions of energy density), it 
follows that the total density, pressure, energy flux and anisotropic pressure are 
given as follows: 



tot 



p 

ptot 



7T. 



tot 
ab 



P + 



k 4 
K 
k 4 



K 

24 



(2p 2 -3n ab n ab )+U 



^ (2p 2 + APp + n ab n ab - Aq a q a ) + 



qT = Qa + 



k 4 



24 



(Apq a - 67c ab q b ) + Q a 



— [-(p + W)ii ab - 37r c(a 7r 6) c + 3q {a q b) ] + V ab 



(3.29) 
(3.30) 
(3.31) 
(3.32) 



Making use of eq. (3.11), the above equations can be recast into the following form 
which is more convenient for our purposes: 



p tot = p+^(2p 2 -3n ab n ab ) + -^U , 

P tot = P + jr(2p 2 + ApP + W ab - Aq a q a ) + \u , 
4 A K A 

1 6 

qT = q a + 4^( 4 p^ - w 6 ) + , 

1 6 

TTaf = n ab + ^t~(P + 3P)n ab + 7l c{a n b) C + q {a q b) ] + ~^Pab 



(3.33) 
(3.34) 
(3.35) 
(3.36) 



It is immediately obvious from the above equations that we can regain the 4- 
dimensional general relativity results, when A -1 — > 0. 
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For the braneworld case, it is useful to introduce an additional dimensionless 
gradient which describes inhomogeneity in the non-local energy density U: 

T a = -DM . (3.37) 
P 

The Gauss-Codazzi scalar equation for the 3-curvature defined by TZ is given 

by 

K = k 2 P (2+P-)+^-U- §6 2 + 2A , (3.38) 

where 

TZ = (3) i? = h ab (3) R ab (3.39) 

with ^R a b the intrinsic curvature of the surfaces orthogonal to w a§ . In braneworld 
models, the Gauss-Codazzi constraint reduces to the modified Friedmann equation 

where the 3-curvature scalar is TZ = 6K/a 2 . In non-flat models (K ^ 0) TZ is not 
gauge-invariant since it does not vanish in the FRW limit. However, the comoving 
projected gradient 

Vb = \VbTZ (3.41) 

is a gauge-invariant measure of inhomogeneity in the intrinsic three curvature of 
the hypersurfaces orthogonal to u a . 



§3.4 Linearised perturbation equations for the 
total matter variables 

In this section we derive the linearized perturbation equations 1 for the study 
of cosmological perturbations (scalar, vector and tensor) on the brane. These 
linearized equations are essential for the study of CMB anisotropies and large 
scale structure, which provide an indirect probe of early universe cosmology. 

§If the vorticity is non-vanishing flow-orthogonal hypersurfaces will not exist, and TZ cannot 
be interpreted as the spatial curvature scalar. 

^The exact, non-linear local and nonlocal conservation equations, propagation and constraint 
equations are found in [135, 136] and these equations can be also used for both cosmological and 
astrophysical modelling, including strong gravity effects. 
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§3.4.1 Local and non-local conservation equations 

Based on the form of the bulk energy-momentum tensor and Z 2 symmetry, the 
brane energy-momentum tensor is still covariantly conserved: 

V h T ab = . (3.42) 

The contracted Bianchi identities on the brane ensure conservation of the total 
energy- momentum tensor, which combined with conservation of the matter tensor 
gives 

V a S ah = k A V a S ab = ^V a S ab . (3.43) 

The longitudinal part of S ab is sourced by quadratic energy-momentum terms in- 
cluding spatial gradients and time derivatives. As a result any evolution and inho- 
mogeneity in the matter fields would generate non-local Coulomb-like gravitational 
effects in the bulk which back react on the brane. The conservation equation (3.42) 
implies evolution equations for the energy and momentum densities, and these are 
unchanged from their general relativistic form. To linear order in an almost-FRW 
brane cosmology we have 

p + e(p + P) + D a q a = 0, (3.44) 

and 

q a + \®q a + (p + P)A a + D a P + D b u ab = . (3.45) 
The linearised propagation equations for U and Q follow from Eq. (3.43) in [135]: 

U + \®U + D a Q a = , (3.46) 

and 

Q a + l<3>Q a + \D a U + D b V a b + \uA a = ^(p + P) (-2D a p + 3D b ir ab + 20q a ) . 

(3.47) 

Taking the projected derivative of Eq. (3.44) we obtain the propagation equa- 
tion for A a at linear order: 

P A a + (p + P) (Z a + aQA a ) + aD a D b q b + aQD a P - QPA a = 0. (3.48) 

From equation (3.46), we obtain the evolution equation of the spatial gradient of 
the non-local energy density: 

*«=(-- I) 0T * " + aQA «) ~ -D a D b Q b . (3.49) 

VP 3/ 3p p 

From the propagation equations for U and Q it can be seen that the energy 

of the projected Weyl fluid is conserved while the momentum is not conserved; 

rather it is driven by the matter source terms on the right of Eq. (3.43). Note that 

no propagation equation for V ab is implied so the set of equations will not close. 
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§3.4.2 Propagation and constraint equations 

In this section we give the general linearised gravito-electric, gravito-magnetic, 
shear and vorticity propagation and constraint equations on the brane, which 
follow from the Bianchi identities, and the equations for the kinematic variables 
17^, and 6 and its gradient Z a which follow from the Ricci identity. 

1. Gravito-electric propagation: 

E ab + QE ab + 7^ 2 (P + P)(?ab + \^K 2 D {a q b) + ^K 2 Q-K ab + ^K 2 TX ab - C\ll\H ab 

-[K 4 {-6p(p + P)a ab + 3{p + 3P)n ab + 3{p + 3P)ir ab 



12/t 2 A' 

- 6pD {a q b) + 0[p + 3P]n ab } - 48U(T ab - 36V ab - 36D {a Q b) - 12QV ab ] ; 

(3.50) 

2. Gravito-electric divergence: 

D h E ab + X -K 2 D\ ab - l -K 2 D a p + ^K 2 Bq a 



1 



8k 2 \ 



^ [ -^P®q a + 2{p + 3P)D\ ab + % P D aP ) + \QD a U - 160Q a - 2AD b 

(3.51) 



3. Gravito-magnetic propagation: 

K 2 3 K 2 

H ab + QH ab + cur\E ab - — curlvr afe = — curlP a6 - ^curl [(p + 3P)-K ab \ ; 

(3.52) 

4. Gravito-magnetic divergence: 

2 2 1 

D b H ab - K 2 (p + P)u a + ^-curlg a = ^cu a + -^-(8Wcu a - 3curlQ a ) ; (3.53) 

Z A K A 

5. Gravito-magnetic constraint 

D( a u b ) + curios - H ab = ; (3.54) 

6. Shear propagation: 

2 1 1 

<?ab + ^CTab + E ab - ~K 2 TX ab - D {a A b) = ^^{« [~(P + ^P)^ab] + l2V ab } ; 

(3.55) 
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7. Shear constraint: 

D b a ab - curler - ^D a Q + n 2 q a = -^-( K A pq a + 6Q a ) ; (3.56) 

3 K z \ 

8. Vorticity Propagation: 

2 1 

u {a) + -Qu a + -curL4 a = ; (3.57) 

9. Vorticity Constraint: 

D a uj a - A a cu a = ; (3.58) 

10. Modified Raychaudhuri equation: 

e = -^0 2 - i/t 2 (p + 3P) + A - ^ [k a P {2 P + 3P) + 12W] + L>M a ; (3.59) 

11. Propagation equation for the comoving expansion gradient >2 a which follows 
from Eq. (3.59): 

2 k 2 1 

Z a +-0Z a -a0A a +yaP a (p+3P)-a J D a P% = -^-^{K 4 aD a [p(2p+3P)]+12aD a U} ■ 

(3.60) 

The above propagation and constraint equations reduce to general relativity 
when A^ 1 — > 0. These linearized equations together with Eqs. (3.42)-(3.47), govern 
the dynamics of the matter and gravitational fields on the brane, and incorporate 
both the quadratic energy-momentum (local) and the projected effects from the 
bulk (nonlocal). The local terms are proportional to p/\, and they are dominant 
at high energies. The nonlocal terms would introduce imperfect fluid effects onto 
the brane, even if the matter has a perfect fluid form. 

It was shown in [135, 136] that the bulk effects give rise to new driving and 
source terms from Eqs. (3.50)-(3.60). However, the vorticity propagation and 
constraint equations, together with the gravito-magnetic constraint have no direct 
bulk effects. The local and nonlocal energy densities act as driving terms in the 
expansion propagation. The spatial gradients of the local and nonlocal energy 
densities provide sources for the gravito-electric field. The nonlocal energy flux 
provides a source for the shear and gravito-magnetic field, and finally, the nonlocal 
anistropic stress acts a driving term in the shear propagation and the gravito- 
electric and gravito-magnetic fields. 

The spatial gradient of the 3-curvature scalar is an auxiliary variable. It can 
be related to the other gauge-invariant variables using Eqs. (3.38) and (3.41): 

Va = K 2 pA a (l + + J^pT a - 2 -QZ a . (3.61) 
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Taking the time derivative of Eq. (3.61), commuting the spatial and temporal 
derivatives, and then making use of Eqs. (3.59) and (3.60), we obtain the evolution 
of the spatial gradient of the 3-curvature scalar: 

r ja +l®r ] a+lK(Z a +aeA a )+leaD a D b A b = -k 2 (l + aD a D b q b --^-aD a D b Q b . 

(3.62) 

In general relativity, propagating i] a is a useful device to avoid numerical instability 
problems when solving for isocurvature modes in a zero acceleration frame (such 
as the rest-frame of the CDM) [126]. 

§3.5 Cosmological Perturbations on the Brane 

In this section, we study the covariant, gauge- invariant splitting of linear cos- 
mological perturbations into scalar, vector and tensor modes. In particular, we 
concentrate on the scalar and tensor equations on the brane, which are used for 
the analysis in the following chapter. The vector perturbations are placed here for 
completeness. 

The limiting case of the background FRW brane is characterized by homogene- 
ity and isotropy, i.e.: 

DJ = V a = W ab = 0, (3.63) 

where the quantities / = p, P, 6, U, V a = A a , u a , Q a , and W ab = a ab , E ab , H ab , V ab . 

The covariant and gauge-invariant splitting into scalar, vector and tensor modes 
is given by: 

V a = D a V + V a , (3.64) 

and 

W ab = D {a D b) W + D {a W b) + W ab , (3.65) 

where a hat denotes a transverse (divergence free) quantity (and W ab is assumed 
trace-free). We note that V a , W ab and the other derived quantities, e.g. / and V a 
are first order. 

§3.5.1 Scalar Perturbations 

The scalar perturbations are covariantly characterized by 

V a = W a = Wab = , (3.66) 

and consequently, we obtain 

curlV; = = curlWU D b W ab = ^D 2 (D a W) , (3.67) 
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The vorticity constraint equation and gravitomagnetic constraint equation show 
that 

cu a = = H ab . (3.68) 

For scalar perturbations, the magnetic part of the Weyl tensor H ab and the 
vorticity tensor uj ah vanish identically. The electric part of the Weyl tensor E ab and 
the shear a ab need not vanish. The non-vanishing variables satisfy the propagation 
and constraint equations on the brane. 



Scalar Harmonics 



The tensor- valued, partial differential equations presented in the earlier sections 
can be reduced to scalar-valued, ordinary differential equations by expanding in 
an appropriate complete set of eigentensors. For scalar perturbations all gauge- 
invariant tensors can be constructed from derivatives of scalar functions . Thus it 
is natural to expand in STF tensors derived from the scalar eigenf unctions of 
the projected Laplacian: 

D 2 QW = -^Q( fe ) , (3.69) 

satisfying Q^> = 0(1)". We adopt the following harmonic expansions of the 
gauge-invariant variables: 

k k 



k k 

E «t> = E K**Q* ' *<* = E ' 
k a k a 



(3.70) 



v. 



a 
k 



Here v$ is the 3- velocity of species % relative to u a ; for the CDM model considered 
here we shall make use of v$ for baryons b and CDM c. For photons 7 and neutri- 
nos v we continue to work with the momentum densities which are related to the 
peculiar velocity of the energy frame for that species by e.g. gi 7) = (4/3)pW^ 7) in 
linear theory. The scalar expansion coefficients, such as are first-order gauge- 
invariant variables satisfying e.g. D a A® = 0(2). 



"The notation 0(n) is short for 0(e n ) where e is some dimensionless quantity characterising 
the departure from FRW symmetry. A list of identities for the scalar, vector and tensor harmonics 
are given in Appendix C. 
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We expand the non-local perturbation variables in scalar harmonics in the 
following manner: 

T a = > 

k 

Q a = Y.P^ ] i (3.71) 

k 

k 

In addition, we can expand the projected gradient of the 3-curvature term: 

* = I>©(l (3-") 

The form of this expansion is chosen so that if we adopt the energy frame (where 
q a = 0) the variable rjk coincides with the curvature perturbation usually employed 
in gauge-invariant calculations. 



Scalar equations on the brane 

It is now straightforward to expand the 1+3 covariant propagation and con- 
straint equations in scalar harmonics (see Appendix C). We shall consider the 
CDM model where the particle species are baryons (including electrons), which 
we model as an ideal fluid with pressure and peculiar velocity Va \ cold dark 
matter, which has vanishing pressure and peculiar velocity vjf\ and photons and 
(massless) neutrinos which require a kinetic theory description. We neglect pho- 
ton polarization, although this can easily be included in the 1+3 covariant frame- 
work [30]. Also, we assume that the entropy perturbations are negligible for the 
baryons, so that D a P^> = c 2 s D a p^ where c 2 s is the adiabatic sound speed. A com- 
plete set of 1+3 perturbation equations for the general relativistic model can be 
found in [34]. We extend these equations to braneworld models here. 

In the following, perturbations in the total matter variables are related to those 
in the individual components by 

pA fc = £>«A« pq k = Y,P®$, pn k = J2p ( H\ (3-73) 

iik 

where = (1 + P^ / p^)v^\ = v£\ and tt^ vanishes for baryons and 
CDM. Similarly, the total density and pressure are obtained by summing over 
components, e.g. P = ^^P^- It is also convenient to write 



P = ( 7 -l)p, 



(3.74) 
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but 7 should not be assumed constant (in space or time). 
We begin with the equation for the gravito-electric field: 



^/t 2 p(70"fc - qk) + jU 2 p©(l - 37)7T fc + ^K 2 p-k k 



1 



— K 



12k 2 \ v 
- (3 7 - 2)p 2 Qn k 



6 ( - ) p (ja k - q k ) - 3(p + 3P)pn k - 3(37 - 2)p(pn k + pn k ) 



- 12 ( - ) (4Ua k - 3pQ k ) - 3Q(pV k + P V k ) - 12pQV k 



(3.75) 



We have written this equation in such a form that every term is manifestly frame- 
independent. The shear propagation equation is 



k 



k 



K 



-[cr k + -Qa k + - ($ fe + A k ) - —pir k = -^v[-(3 7 - 2)^p l ^ k + 12pV k ] . 



Ak 2 \ [ 



The shear constraint is given by 



2 (k 



2 r 



Z k - 1 - 



3K 



k PVk - - - 
3 V a 



The gravito-electric divergence is 

3 



K 2 \ 



(3.76) 



K A p 2 q k + 6pQ k ). (3.77) 



1 3X \ i 



8k 2 X 



--p 2 eqk + -(3 1 -2)p 2 (l- 



, , 3X 

3K\ k 



tv 1, + ir(r)p<i k 

"TTfe + o~P A fc 

; 3 a 



(3.78) 



+ 16-pr fe - 160pQ fe - 16p 



3^ , „ 



a 



The propagation equation for the comoving expansion gradient Z a is given by 

2(pW A« + pMA?) + (1 + 3c2)p<» A< 6) + p^A? 
K 4 [(2p + 3P)pA fe + p(3pW A™ + 3p^A^ + (2 + 3c 2 )p^ A? + 2p( c )A< c) )] 



(3.79) 



1 a A; 

Z fc + -6Z fc - -04 fc + -A k + y ^ 

1 a 



2n 2 \k 
+ 12pT fe l . 
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The non-local evolution equations for T k and Q k are 



T fc = o ( 3 7 - 4)©T fc - -<d-A k - Z k + -Q k , 

S S p Spa a 



(3.80) 



and 



& -1(37- 4)00, + ^ 



T fe + 2(l-^]n 



AkU A 

+ 7T A k 

Sap 



K 4 ( k 

—IP < Qqk + ~ 
6 ! a 



3K , 



(3.81) 



The spatial gradient of the 3-curvature scalar is 



3;n « 2 p A /. , p 



' P T k - ~GZ k , (3.82) 
3 a 



and it evolves according to 



^~ 1 - 2 -" -^(^Afc + BpQfc). (3.83) 



The evolution equations for the scalar harmonic components of the comoving, 
fractional density gradients for photons, neutrinos, baryons and cold dark matter 
(CDM) are 



» 



.(*>) 



k /4 c A 4 

-- f -Z fc - q^j - -GA k (photons) , 

~ (^ z k ~ ~ \® A k (neutrinos) , 



(3.84) 
(3.85) 



1 + 



p(b) 



--{Z k - v 
a 



0>h 
k > 



+ 



p(b) 



c^GAf (baryons), 
(3.86) 



--(Z k -v { k c) )-QA k (CDM). 

(X 



(3.87) 
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The evolution equations for the momentum densities and peculiar velocities are 



•(7) 
Qk 



•O) 
Qk 



lk 
3 a 

lk 

3 a 



A< 7) + 4A k + 2 ( 1 - 



+ AA k + 2(1 



3K 
~k? 

3K 



71 



(7) 



(b) 



,(7) 



+ n e a T \-v k - q k 



7T 



k 



(1 - 3^)6^ + - Ak 



(b) 



,(7) / 4 



- n e a T - 



.(*) 



3 w fc - <?fc 



(7) 



3 U% a fe ' 



(3.88) 
(3.89) 

(3.90) 
(3.91) 



where the Thomson scattering terms involving the electron density n e and Thom- 
son cross section a? arise from the interaction between photons and the tightly- 
coupled baryon/electron fluid. The remaining equations are the propagation equa- 
tions for the anisotropic stresses of photons and neutrinos, and the higher moments 
of their distribution functions. These equations can be found in [34] , and with po- 
larization included in [29, 30, 32], since they are unchanged from general relativity. 
However, we shall not require these additional equations at the level of approxi- 
mation we make in our subsequent calculations. 



§3.5.2 Vector Perturbations 

The vector perturbations are covariantly characterized by 

V a = V a , W ab = D {a W bh curlAJ = -2fu a , 
and it follows that 

D a W ab = ^D 2 W b , cm\W ab = ^D {a cuT\W b) . 



(3.92) 



(3.93) 



Vector equations on the brane 

There are no bulk effects in the linearized vorticity propagation equation (3.57), 
and hence the vorticity decays away with the expansion like the general relativity, 
and this can be interpreted as angular momentum conservation equation on the 
brane. However the gravito-magnetic divergence equation becomes 

K 2 1 

D b H ab = k\ P + P)u a - — curlg a + — {k 2 p(p + P) + SU}u a - 3curlQ a . (3.94) 
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The bulk terms provide additional sources for the gravito-magnetic field. The 
existence of the bulk terms made it possible to source vector perturbations even 
when the vorticity vanishes, since curlQ a may be non-zero. 



§3.5.3 Tensor Perturbations 

The (l+3)-covariant description of gravitational waves in a cosmological con- 
text has been considered in detail by Challinor [31] and Dunsby et al [54] for the 
non-flat and flat cases respectively. In this section, we will extend to the case of 
braneworld cosmology. 

Tensor perturbations are covariantly characterized by 



DJ = 0, A a = u a = Q a = 0, D a W ab = 0, W ab = W ab . (3.95) 



In this description, the linearized gravitational waves are described by the 
tranverse degrees of freedom in the electric and magnetic parts of the Weyl tensor. 
D a W ab indicates that the shear and anisotropic stress (both local and nonlocal) are 
transverse. Furthermore, the vorticity and all projected vectors vanish at linear 
order. The individual matter components all possess the same four-velocity which 
defines the fundamental velocity u a in a pure tensor mode. 

The constraint equation (Eq. (3.54)) becomes 



which determines the magnetic part of the Weyl tensor from the shear. 

For convenience, we define a rescaled nonlocal anisotropic stress n* ab by the 
following relation: 



H ab = cm\a ab , 



(3.96) 




6 



V ab . 



(3.97) 



The propagation equations for the tensor perturbations from Eqs (3.50)-(3.60) 
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are reduced to the following set: 

p = -( p + P)e, 



(3.98) 

6 = -ie 2 - 1 -k\ P + 3P) + A - 2^[«V(2p + 3P) + 12W] , (3.99) 
2 k 2 k 2 

C^b = -gQo-afc - #ab + y (^ab + < b ) ~ ^(P + 3P)7T afe , 

ij a6 = -0# a6 - curlP ab + ycurl(7r a6 + 7r* 5 ) - ^curl [(p + 3P)lX ab } , (3.101) 

-2 r j j 

(p + P)(T afe + ?r ab + -07r a6 + n* ab + -©7T* 6 



(3.100) 



#ab = -©^ab + CUrlH ab - 



(3.102) 



+ 



3/€ 4 {(p + 3P)7T afe + (p + 3P)7T a6 } + K\p + 3P)07T afe 



12k 2 A 

- 6/t 4 p(p + ^>ab - 48Wa a6 

The inhomogeneous wave equations for the shear and the magnetic part of the 
Weyl tensor follow from differentiating Eqs. (3.100) and (3.101) along the flow lines 
of u a , making use of Eq.(3.74) (assuming 7 7^ 0) and the identities (C.10)-(C.13) 
from Appendix C, we obtain: 



<jab + D 2 a ah + -9& ab + 



i(4-37)*V-f 



&ab 



(37 - 2)^- ^37T a6 - (37 - 2)07T a6 + ^ ab - 3pO"a6 



(3.103) 



and 



2 ? ■ 

-f^ab + -D i/ a b + ~®Hab + 2 



, 2 3iT 



ab 



curl7r ab + -6curl7r a6 



+ « 2 



curl7r* 5 + -0curl7r* b 



K 



u(S7-2)p 



Tensor Harmonics 



curl7r afc - 7 0curl7r ab + curl7r a5 + -0pcurl7r afe 
07 — z / o 



(3.104) 



Analogous to the scalar case, the electric and magnetic parts of the Weyl tensor, 
the shear, and the local and non-local anistropic stress can be expanded directly 
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in the electric and magnetic parity tensor harmonics defined in [31]: 

k ^ ' 

= E (I) M£> + 9 >Q$) • 

k 



(3.105) 
(3.106) 

(3.107) 
(3.108) 
(3.109) 



where J2 k denotes a sum over the harmonic modes. Note that the electric and 
magnetic parity tensor harmonics are related by Eqs (C.14) and (C.15) in Appendix 
C. 



Tensor Equations on the Brane 

Using the gravitomagnetic constraint in Eq. (3.54) and together with tensor 
harmonics above, we expand the first order equations (eqs. (3.100) and (3.102)), 
and obtain ** , 



( &k + 1^) + (£) Ek - Y p7Tk = ^x p[ - K ^ - 2)p1Tk] + Y pnl ' 



and 



fcV 3K k 2 



(3.111) 



K 2 K 2 

°k + yP 7 ^ ~ y ( 3 7 - l)p@7Tfe 



12k 2 \ 



K 



6 ( - ) 7P cr k - 3(f) + 3P)p-K k - 3(37 - ^)p(p^k + P^k) - (37 - 2)p Qn k 



48 I ~) Ua k \ ~ y (K + PK + \f&*l 



'Note that from eq. (3.96), we find that H k is algebraic in 



(3.112) 



(3.110) 



These curl terms lead to a coupling of different polarization states. 
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From Eq.(3.103), we get 

1 

-(37 - 2)QiT k - 7r fe 

9(2 7 2 - A ?)K 4 p 2 ir k - 3(37 - 2)K 4 p 2 7r k - 2(2 - 97)/t 4 p 2 07r fe . 

(3.113) 

This equation is equivalent to the form found by Maartens in [135, 136] for the 
tranverse traceless modes on the brane. From the above equation, the nonlocal 
bulk effects would provide driving terms which would act like the anisotropic stress 
in general relativity. For the local anisotropic stress, the evolution is determined 
by the Boltzmann equation. However, the evolution equation of the non-local 
anisotropic stress is not known. We will make an local approximation to the 
evolution equation of the non-local anisotropic stress in the subsequent chapter to 
see its effects on the CMB tensor power spectra. 



a k + 0cr fc + 



a 2 a 2 3 



a k + K P [ £ 



+ 



1 

6/t 2 A \k 



§3.6 A covariant expression for the temperature 
anisotropy 

In this final section, we discuss the line of sight solution to the Boltzmann 
equation for the scalar contribution to the gauge-invariant temperature anisotropy 
<5r(e) of the CMB in braneworld models. We employ the 1+3 covariant approach, 
and show that our result is equivalent to that given recently by Langlois et al [118] 
using the Bardeen formalism. 

Over the epoch of interest the individual matter constituents of the universe 
interact with each other under gravity only, except for the photons and baryons 
(including the electrons), which are also coupled through Thomson scattering. 
The variation of the gauge-invariant temperature perturbation &r(e), where e a is 
the (projected) photon propagation direction, along the line of sight is given by 
the (linearized) covariant Boltzmann equation (valid for scalar, vector, and tensor 
modes) [33]: 



5 T (e)' + a T n e 5 T (e) = -a ab e a e b - A a e a - - 



(3.114) 
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where the prime denotes the derivative with respect to a parameter A defined along 
the line of sight by d\ = —u a dx a . 

Following the steps in [33], we expand the right-hand side of Eq. (3.114) in 
scalar harmonics and integrate along the line of sight from the early universe to 
the observation point R. Neglecting effects due to the finite thickness of the last 
scattering surface, on integrating by parts we find that the temperature anisotropy 
involves the quantity 

GHO' + \\ {ak ~ Zk) + A ' k ~ HAk = ~ 2 ^ k + © 2 1 (3 - 115) 

integrated along the line of sight (after multiplying with Q^). In simplifying 
Eq. (3.115) we have made use of the derivative of the shear propagation equation 
(3.76), substituted for q k and Z k from Eqs. (3.75) and (3.77), and finally used 
Eqs. (3.40) and (3.59). The quantity I is the total sum of all the braneworld 
corrections: 



/ = 



a\ 2 
k) 



/ 1 + Ie/ 1 + 7 2 + ^ 



A: 



-cr k I, 



(3.116) 



where 
h = 
h = 



1 



An 2 \ 
1 



[-(37 - 2)K 4 p\ k + 12pV k ] , 
k 



6 



12k 2 X ^ 
- (37 - 2)p 2 07r fc 
1 



7P a k - 3(p + 3P)p7r k - 3(37 - 2)p(pii k + pn k ) - 6 



48 



Ua k - 36(pV k + pV k ) + 36 - ) pQ k - 12pBV k 



2k 2 X 



[(3 7 - l)« 4 p + 12U] , 



n 2 p 2 4 7/ 

h = ~^- a k H — 7rrUa k - 



3A 



K 2 \ 



Ak 2 \ 



(4K*p 2 q k + 2ApQ k ) 



(3.117) 



A lengthy calculation making use of the propagation and constraint equations 
shows that 1 = 0. The final result for the temperature anisotropies is then 



[*r(e)]* = -£ 



A™ + \a k + A k ) Q( fc ) 



16 



(3.118) 



A.i 



where the event A is the intersection of the null geodesic with the last scattering 
surface. 
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In retrospect, one could re-derive the result for the temperature anisotropy 
in braneworld models much more simply by retaining the effective stress-energy 
variables p tot , P tot 5 g* ot and ir^f in the propagation and constraint equations used 
in the manipulation of the left-hand side of Eq. (3.115), rather than isolating the 
braneworld contributions. 

If we adopt the longitudinal gauge, defined by a a b = 0, we find that the electric 
part of the Weyl tensor and the acceleration are related by = — if the total 
anisotropic stress 7r*^ vanishes. It follows that in this zero shear frame we recover 
the result found by Langlois et al [118]. 



Chapter 4 



Braneworld Cosmology II: 
Initial Conditions and CMB 
Anisotropics 

"Order and simplification are the first steps toward the mastery of a 
subject - the actual enemy is the unknown. " 

- Thomas Mann 

In this chapter, we discuss the possible implications of braneworld cosmology for 
the CMB, based on work described in [123] and [124]- For the scalar perturbations, 
we supplement the equations for the total matter variables with equations for the 
independent matter components in a cold dark matter cosmology, and provide solu- 
tions in the high and low-energy radiation- dominated phase under the assumption 
that the non-local anisotropic stress vanishes. These solutions reveal the existence 
of new modes arising from the two additional non-local degrees of freedom. Our 
solutions should prove useful in setting up initial conditions for numerical codes 
aimed at exploring the effect of braneworld corrections on the cosmic microwave 
background ( CMB) power spectrum. We show that the 3-curvature is only constant 
in the high energy limit for the modes which correspond similarly to the growing 
and decaying mode in general relativity. The CMB temperature anisotropics are 
insensitive to the brane tension A if the dark energy contribution is ignored. For 
the tensor perturbations, we set out the framework of a program to compute the 
tensor anisotropics in the CMB that are generated in braneworld models. In the 
simplest approximation, we show the braneworld imprint as a correction to the 
power spectra for standard temperature and polarization anisotropics and similarly 
show that the tensor anisotropics are also insensitive to high energy effects. 
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§4.1 From Perturbation Theory to CMB Anisotropics 
in Braneworlds 

The cosmic microwave background (CMB) currently occupies a central role 
in modern cosmology. It is the cleanest cosmological observable, providing us 
with a unique record of conditions along our past light cone back to the epoch of 
decoupling when the mean free path to Thomson scattering rose suddenly due to 
hydrogen recombination. Present (e.g. BOOMERANG [115], MAXIMA [14], and 
VSA [164]) and future (e.g. MAP and PLANCK) data on the CMB anisotropies, 
combined with large-scale structure data, provide extensive information on the 
amplitude and evolution of cosmological perturbations. Potentially this allows us 
to infer the spectrum of initial perturbations in the early universe and to determine 
the standard cosmological parameters to high accuracy. An obvious question to ask 
is whether there are any signatures of extra dimensions which could be imprinted 
on the cosmic microwave sky. 

There has been an explosion of interest in the theory of cosmological perturba- 
tions in braneworlds and their implications for observational cosmology. We will 
briefly summarize * most of the work which has explored these issues. 

The relation between braneworld cosmology, the AdS/CFT correspondence and 
quantum cosmology has been discussed in [4]. Binetruy et al [12] have studied the 
background cosmological dynamics of a Friedmann brane in a Schwarzschild-Anti- 
de Sitter bulk. The modifications to inflation are also explored by various authors 
in [41, 71, 75, 128], in [95] (for quintesscence) and in [76, 116, 166] (for gravitational 
waves). The large scale perturbations generated from quantum fluctuations during 
de Sitter inflation on the brane have been computed and studied in [20, 71, 80, 
117, 141]. High-energy inflation on the brane generates a zero-mode (4D graviton 
mode) of tensor perturbations, and stretches it to super-Hubble scales. This zero- 
mode has the same qualitative features as in general relativity, remaining frozen at 
constant amplitude while beyond the Hubble horizon, but the overall amplitude is 
higher [119]. The massive KK modes (5D graviton modes) remain in the vacuum 
state during slow-roll inflation. The evolution of the super-Hubble zero mode is 
the same as in general relativity, so that high-energy braneworld effects in the 
early universe serve only to re-scale the amplitude. However, when the zero mode 
re-enters the Hubble horizon, massive KK modes can be excited. Qualitative 
arguments [76, 80] indicate that this is a very small effect, but it remains to be 
properly quantified, so that the signature on the CMB may be calculated, and 
constraints may be imposed on the braneworld parameters. 

Other authors [19, 50, 100, 101, 107, 108, 110, 111, 165, 166, 178, 187] have 
considered the cosmological consequences of interactions between the brane and 

*We apologize in advance if we have missed some of the literature in this thesis. 
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the bulk. These authors constructed approximations from the bulk equations and 
explored the implications for observational cosmology on the brane. Their methods 
are based mainly on the Bardeen approach. 

For the study of cosmological perturbations and its relation to the cosmic 
microwave background (CMB) and large scale structure (LSS), the covariant for- 
malism has been well described in the previous chapter. We will summarize some 
of the main results for the scalar and tensor perturbations based on earlier work 
by various authors. 

For scalar perturbations, the bulk effects introduce a non-adiabatic mode on 
large scales. The additional non-adiabatic has been found in [21, 77, 118, 123, 135]. 
In additional, another growing mode was found in [123] in ACDM cosmology. The 
density perturbations on large scales can be solved on the brane without solving 
for the bulk perturbations [77], but the Sachs- Wolfe effect cannot be found on 
the brane because of the non-local anisotropic stress (see also [118]), which is 
underdetermined. There are possible changes on the scalar modes for both large 
and small scales. 

For tensor perturbations, we recall earlier that the bulk effects generate a mass- 
less mode during inflation and a continuum of KK modes [71, 117]. The massive 
modes stay in the vacuum state, and on large scales, there is a constant mode with 
enhanced amplitude in [117]. In this case, there is no qualitative change on large 
scales in the low energy regime (4D general relativity), but significant change in 
small scales. 

We have learned from the previous chapter that there are two types of cor- 
rections which will arise from braneworld cosmology. The first type of correc- 
tion occurs at energies well above the brane tension A, where gravity becomes 
5-dimensional. This will introduce significant corrections to general relativity. On 
the other hand, there are also corrections that can operate at low energies, medi- 
ated by bulk graviton or Kaluza-Klein (KK) modes. Both types of correction play 
an important role in both scalar and tensor perturbations. 

In §4.2, we work out the perturbative dynamics and initial conditions for the 
scalar perturbations in both the CDM (§4.2.1) and energy frames (§4.2.2). We 
first show that there are two additional non-local isocurvature modes present in 
braneworld cosmology and subsequently examine whether they will contribute sig- 
nificantly to the CMB anisotropy in §4.2.3. 

We develop the formalism to compute the tensor anisotropies in the CMB in 
§4.3, which incorporates the early-universe high-energy braneworld effects, and we 
carefully delineate what is known on the brane from what is required from bulk 
equations. Once the 5D solutions are provided, our formalism, with its modified 
CMB code (based on CAMB [125, 126]), is able to compute these anisotropies. We 
illustrate this by using a simple approximation in §4.3.1 to the 5D effects (refer 
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to the analysis of the braneworld scalar Sachs- Wolfe effect in [9]). We plot the 
braneworld correction to the CMB power spectra for temperature and polarization 
anisotropies and discuss their implications in §4.3.2. 



>4.2 Perturbation dynamics and Initial Condi- 
tions of Scalar Modes 

4.2.1 The CDM Frame 



In this section we specialize our equations to FRW backgrounds that are spa- 
tially flat" 1 " and we ignore the effects of the cosmological constant in the early 
radiation-dominated universe. To solve the equations it is essential to make a 
choice of frame u a . We adopted a frame comoving with the CDM similar to [34]. 
Since the CDM is pressure free, this u a is geodesic (A a = 0) which simplifies the 
equations considerably. We shall adopt this frame choice here also, though we 
note it may be preferable to use a frame more closely tied to the dominant matter 
component over the epoch of interest. This can be easily accomplished by adopt- 
ing the energy frame (q a = 0). For completeness, we give equations in the energy 
frame in the following section. 

We neglect baryon pressure (c 2 s — > and — > 0) and work to lowest order in 
the tight-coupling approximation (n e a T — > oo; see e.g. [133]). At this order the 
energy frame of the photons coincides with the rest frame of the baryons, so that 
Va^ = 3gi 7 ' ) /(4p^ 7 - ) ), and all moments of the photon distribution are vanishingly 
small beyond the dipole. 

With these approximations and frame choice we obtain the following equations 
for the density perturbations of each component: 




(photons) , 



(neutrinos) , 



(4-1) 





k 



(Z k -v ( k >) (baryons), 



(4.3) 



a 



k 



Z k (CDM). 



(4.4) 



a 



TMore generally, curvature effects can be ignored for modes with wavelength much shorter 
than the curvature scale, k 3> \/\K |, provided the curvature does not dominate the background 
dynamics. 
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The equations for the peculiar velocities and momentum densities are 



(4pW + 6p^)q i k l) 



•O) 
Qk 



4k 

3 a 
lk 



-—(AfN 27T. 
3 a 



P (6) e^ } 



» 



(4.5) 
(4.6) 



along with = and = 6q^ /4. The latter equation, together with Eqs. (4.1) 
and (4.3), implies that Aj^ = 3A^/4 so that any entropy perturbation between 
the photons and baryons is conserved while tight coupling holds. The effects 
of baryon inertia appear in Eq. (4.5) because of the tight coupling between the 
baryons and photons. 

The constraint equations are found to be: 



and 



K 2 pq k 



U-) (Z k -a k )- l -(-\ (K 4 p 2 q k + 6pQ k ) = 0, (4.7) 



6 V K 



21 * 

a 

1 



K 



K p [ ~ ) (A fc - 7T fc ) + K Qpq k 



16 V k 



K 



-\®p 2 q k + \~P 2 [(3 7 - 2)n k + 2A k 
6 6a 



(4.8) 



+ 16 ( - ) p(T k - V k ) - IQOpQk 



The propagation equation for the comoving expansion gradient in the CDM frame 
is 



k 2 a 



2(p^A^ + p^A^) + p^A^+p^A^ 



Zk+ 6 eZk+ i 2kl 

= ~ (?) 4 lUl(2p + 3P) P A fc + p(3pWAW + 3pMAM + (2 + 3c>^ + 2p( c )A( c) )] 



+ 12pT fe 



(4.9) 



The variables Q k and a k can be determined from the constraint equations so their 
propagation equations are not independent of the above set. The propagation 
equation for § k is unchanged from Eq. (3.75) since that equation was already 
written in frame-invariant form. The propagation equation for the shear in the 
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CDM frame is 

{b k + le<7 fc ) + (^j $ fe - yp7T fc = ^ Q [-(37 - 2)«V7T fc + 12pP fc ] • 

(4.10) 

Finally we have the non-local evolution equations for T k and Q k which in the 
CDM frame become 

t fc = ^(37 - 4)6T fe - \--Z k + -Q k , (4.11) 
3 Spaa 



and 



4 



Qfe - 3(37 - 4)0Q fe + --(T fc + 2V k ) = - 1P 



k 

<S>qk + -(7r fc - A fc ) 

a 



(4.12) 



Solutions for the CDM Frame in the radiation-dominated era 



We now use the above equations to extract the solutions of the scalar pertur- 
bation equations in the radiation-dominated era, 7 = 4/3. To simplify matters, 
as well as neglecting the contribution of the baryons and CDM to the background 
dynamics, we shall only consider those modes for which D a p^> and D a p^ make a 
negligible contribution to the total matter perturbation D a p. This approximation 
allows us to write the total matter perturbations in the form 

(pW + p^)A k = pMAW + pMAM (pW + pM)q k = p™q™ + p^q { k \ (4.13) 

and effectively removes the back-reaction of the baryon and CDM perturbations on 
the perturbations of the spacetime geometry. We note that in making this approx- 
imation we lose two modes corresponding to the baryon and CDM isocurvature 
(density) modes of general relativity, in which the sub-dominant matter compo- 
nents make significant contributions to the total fractional density perturbation 
(which vanishes as t — > 0). However, for our purposes the loss of generality is not 
that important, while the simplifications resulting from decoupling the baryon and 
photon perturbations are considerable. We also neglect moments of the neutrino 
distribution function above the dipole (so there is no matter anisotropic stress). 
This approximation is good for super-Hubble modes, but fails due to neutrino free 
streaming on sub-Hubble scales. 

We shall also assume that the non-local energy density U vanishes in the back- 
ground for all energy regimes [77]. Physically, vanishing U corresponds to the 
background bulk being conformally flat and strictly anti-de Sitter. Note that 
U = in the background need not imply that the fluctuations in the non-local 
energy density are zero, i.e. T a 7^ 0. 
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With the above conditions the following set of equations are obtained: 



- ) Q k + 3HV k - V k 



' ^ (Z k + HZ k ) + /t 2 p ( 1 + y 1 A/, 



(a k + Ha k ) + ) <!»/, 



gfe + 3a fc 



lib 



' + 3 a Afc 
AW , fc f 4 7 n (7) 



* - «r i = 



6 p 
A 



J3_p 

tfX fe ' 

o, 
o, 
o, 

o, 



(4.14) 
(4.15) 

(4.16) 

(4.17) 
(4.18) 
(4.19) 

(4.20) 



where recall H = 0/3. For the constraint equations we find 



3 K *(i + ^-2 0)\iW*) = -i^Q*, 



2 £ 



3#<?* 



-) 

a 



6 p 
^ A 



(4.21) 



-)(T fc -7> fc )-3ifQ fc 



(4.22) 



Finally the non-local evolution equations are found to be : 



Tfc — — Qk , 

a 

9Q k + 3[ ^ (T fc + 2V k ) = -2k a P f^A k - 3Hq k 



(4.23) 
(4.24) 



It is easy to show by propagating the constraint equations that the above set of 
equations are consistent. 
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By inspection, there is a solution of these equations with 



A 



<?k = , 



(7) 
Qk = 

<A = 



T fc 



a 
A A 

3 a ' 
AA 

3 a ' 

-«4 



0, 
0, 
0, 



(4.25) 
(4.26) 

(4.27) 
(4.28) 
(4.29) 
(4.30) 

(4.31) 

(4.32) 
(4.33) 
(4.34) 



where A is a constant. This solution describes a radiation-dominated universe 
that is exactly FRW except that the CDM has a peculiar velocity vjjp = (A/a)Q^ 
relative to the velocity of the FRW fundamental observers. [This form for 
clearly satisfies Eq. (3.91) with A a = 0.] Such a solution is possible since we 
have neglected the gravitational effect of the CDM (and baryon) perturbations 
in making the approximations in Eq. (4.13). The same solution arises in general 
relativity [34]. Including the back- reaction of the CDM perturbations, we would 
find additional small peculiar velocities in the dominant matter components which 
compensate the CDM flux. We shall not consider this irregular CDM isocurvature 
velocity mode any further here. 

Another pair of solutions are easily found by decoupling the photon/neutrino 
entropy perturbations. Introducing the photon/neutrino entropy perturbation (up 
to a constant) A 2 and relative flux q 2 : 



Al 7) - A fe 



(7) 

q 2 = q k v 



(") 
Qk 



the equations for A 2 and q 2 decouple to give 

k 

A 2 - -q 2 = , 

a 

1 k 

q 2 + --A 2 = . 
3 a 



(4.35) 

(4.36) 
(4.37) 
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Switching to conformal time (dr = dt/a) we can solve for A 2 and q 2 to find 

q 2 ( T ) = B cos (^] + C sin (^] , (4.38) 



A 2 (r) = B sin {^j - C cos {^j . (4.39) 

The constants B and C label the neutrino velocity and density isocurvature modes 
respectively [24, 34], in which the neutrinos and photons initially have mutually 
compensating peculiar velocities and density perturbations. The perfect decou- 
pling of these isocurvature modes is a consequence of our neglecting anisotropic 
stresses (and higher moments of the distribution functions) and baryon inertia. 

Having decoupled the entropy perturbations, we write the remaining equations 
in terms of the total variables A& and q k . The propagation equations for the non- 
local variables and Qk are redundant since these variables are determined by 
the constraint equations (4.21) and (4.22): 

« £ Tl = 2 g)* •> + 2H g) (* - at) - .V (l + £) A* + « (4.40) 

Substituting these expressions in the right-hand sides of Eqs. (4.14) and (4.15) we 
find 

i) 2 (* t + att ) + g) (l + £) . _ i g) 3 w _ „» . l£ (3 ,n 

(4.42) 

;) + "*) + ^ (t) a * + 2 (£> 2 *' + 2ff (;) <* - = "I* • 

(4.43) 

fc\ . / ' A;\ 2 3 p 



<7 fc + if<7 fc ) + - = — ^-Pfc , 

(4.44) 

fc /4 \ 
A fe + - ( -Z k -q k J = , 



1 /c . 

g fc + --A fc = . 
3 a 



(4.45) 
(4.46) 
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These equations describe the evolution of the intrinsic perturbations to the brane. 
The usual general relativistic constraint equations are now replaced by the con- 
straints (4.40) and (4.41) which determine two of the non-local variables. The lack 
of a propagation equation for V k reflects the incompleteness of the 1+3 dimensional 
description of braneworld dynamics. 

In the following it will prove convenient to adopt the dimensionless independent 
variable 

, = A, (4.47) 

which is (to within a factor of 2-k) the ratio of the Hubble length to the wave- 
length of the perturbations. Using the (modified) Friedmann equations for the 
background in radiation domination, and with U = 0, we find that 

dx = k ( 2 + 3p/X \ 

dt a V 2 + p/X J ' v ' ; 

The relative importance of the local (quadratic) braneworld corrections to the 
Einstein equation depends on the dimensionless ratio p/\. In the low-energy 
limit, p < A, the quadratic local corrections can be neglected although the non- 
local corrections S a t may still be important. In the opposite (high-energy) limit 
the quadratic corrections dominate over the terms that are linear in the energy- 
momentum tensor. We now consider these two limits separately. 

Low-energy regime 

In the low-energy regime we have dx/dt ~ k/a and x ~ kr. The total energy 
density p is proportional to x~ 4 . Denoting derivatives with respect to x with a 
prime, using p <C A, and assuming that we can neglect the term involving (p/X)A k 
in Eq. (4.43) compared to the other terms, we find 

27 

3x 2 $' fc + 3x<S> k + (6 + x 2 )a k - x 2 Z k = - T r{3V k - xV' k ) (4.49) 
x 2 Z' k + 3xZ k + 2x 2 $ fe - 2xa k = -^rV k (4.50) 

K A 

9 

x 2 a' k + xa k + x 2 $ k = —V k (4.51) 

K A 

A' fe + ^Z k - q k = (4.52) 

q' k + ^A k = (4.53) 
Combining these equations we find an inhomogeneous, second-order equation for 

3x^1 + 12$' fc + x$ k = F k (x) , (4.54) 
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where 



F k (x) 



27 



K 4 X 



V 'k - — + "3 - "2 + - ) V " 

x Xx 6 x A x 



(4.55) 



In general relativity the same second-order equation holds for $ fc but with F k (x) = 
0. 

The presence of terms involving the non-local anisotropic stress on the right- 
hand side of Eq. (4.54) ensure that cannot be evolved on the brane alone. The 
resolution of this problem will require careful analysis of the bulk dynamics in five 
dimensions. In this thesis, our aims are less ambitious; we shall solve Eq. (4.54) 
with Vk = 0. Although we certainly do not expect V a b = 0*, the solutions of the 
homogeneous equation may still prove a useful starting point for a more complete 
analysis. For example, they allow one to construct Green's functions for Eq. (4.54) 
which could be used to assess the impact of specific ansatze for V a b [9]. 

With Vk = we can solve Eqs. (4.49)-(4.53) analytically to find 



x 6 



V3 



3 sin f —= ) — xV3 cos ( — = 



V3 
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x \ I x \ Ac% 4 
A k = c 4 cos I —= + c 5 sin —= H + - 




x" 



X 



C\ COS 
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c 5 / x \ c 4 / x \ 4c 3 Ax 

?l= vi cos l7sJ"vi sm ivsJ + ^ + vs 



/ x 



C2 sin 



3 cos ( —= | + xy/3 sin ( —= . 

V3J \VsJ\ 



x z 



c 2 cos 



X 



v/3 

(X \ fx 
—= + ci sin —= 

x 

7s 



x \ 

• \ cos I — = J — c 2 sin 



+ 



x \ .fx 
c 2 cos I —= + ci sin — F 
WSJ Vv/3 



(4.56) 
(4.57) 

+ Ci sin 
(4.58) 

(4.59) 
(4.60) 



The mode labelled by C3 is the CDM velocity isocurvature mode discussed earlier. 
The modes labelled by C\ and c 2 are the same as in general relativity; they describe 
the adiabatic growing and decaying solutions respectively. However, in the low- 
energy limit we also find two additional isocurvature modes (c 4 and c 5 ) that are 



*We have not investigated the consistency of the condition V a b — with the five-dimensional 
bulk dynamics in the presence of a perturbed brane. 
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not present in general relativity. These arise from the two additional degrees of 
freedom and Qj, present in the braneworld model (with = 0). The mode 
c 4 initially has non-zero but compensating gradients in the total matter and non- 
local densities, and C5 initially has compensated energy fluxes. Formally these 
isocurvature solutions violate the assumption that the term involving (p/A)Afc be 
negligible compared to the other terms in Eq. (4.43) since all other terms vanish. 
In practice, there will be some gravitational back-reaction onto the other gauge- 
invariant variables controlled by the dimensionless quantity p/A, but the general 
character of these isocurvature modes will be preserved for p/A 1. 



High-energy regime 

We now turn to the high-energy regime, where the quadratic terms in the 
stress-energy tensor dominate the (local) linear terms. In this limit the scale 
factor a oc t 1 ^. The modification to the expansion rate leads to an increase 
in the amplitude of scalar and tensor fluctuations produced during high-energy 
inflation [9]. With U = in the background, and p> A, the Hubble parameter is 
approximately 

H2 * S* V = % < 4 ' 61 > 

and dx/dt rs 3k/ a. In terms of conformal time r, x ~ 3kr. The total energy 
density, p is proportional to re*" 4 / 3 ). 



Power series solutions for the high-energy regime 

It is convenient to rescale the non-local variables by the dimensionless quantity 
/t 4 p. Thus we define 
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The fractional total (effective) density perturbation and energy flux can be written 
in terms of the barred variables [e.g. Y a = Y a /(/t 4 p)p in the high-energy limit as 

"' ~'--2(X, ■ (iT„). (4.65) 



-.tot 



,tot 



2 if' 

(9a + 6e.). 



(4.66) 



Making the high-energy approximation p ^> A in Eqs. (4.42)-(4.46), we obtain 
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-36P fc , 




18- 
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o, 




. 





3a;cx fc + o k + a;$fc = 18 

. / 1 4^ 

fe ~ 3 qk + 9 * 
, 1 A 

The non-local quantities Y fe and Q k are determined by the constraints 

Tfc = To x2 ® k + VE x ( Z k ~ crk) ~ ^A fc + V k , 
lo lo O 

Qk = -^x 2 (Z k - o k ) - ^q k . 



(4.67) 
(4.68) 
(4.69) 
(4.70) 
(4.71) 

(4.72) 
(4.73) 



We can manipulate Eqs. (4.67)-(4.71) to obtain a fourth-order equation for the 
gravitational potential $ k : 



729x 



<9a; 4 



+ 3888a; — £ + (1782 + 54a; 2 )' 



+144a;^ + (90 + a; 2 )$ fe = F fe (a;) , 

(4.74) 



where 



54 



ar 



243a; 



dx 4 



- 810a; 3 — -± + 18a; 2 (135 + 2a; 2 
ox 



dx 2 



-30a;(162 + a; 2 )- 



dx 



+ [x 4 + 30(162 + x 2 )}P k . 



(4-75) 

§ A general and useful identity for these variables is employed in deriving the equations for 
the high energy regime: 

K 4 p 
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Since we do not have an evolution equation for Vk we adopt the strategy taken in 
the low-energy limit and look for solutions of the homogeneous equations (Vk = 
0). In principle one can use these solutions to construct formal solutions of the 
inhomogeneous equations with Green's method. 

To solve Eq. (4.74) with Vk = we construct a Frobenius (or power) series 
solution for 



For each value of m we substitute into Eq. (4.74) and solve the resulting recursion 
relations for the {a n }. We then obtain the other gauge- invariant variables by 
direct integration. The original set of equations (4.67)-(4.71) has five degrees of 
freedom, so we expect one additional solution with = 0. This solution is the 
CDM isocurvature solution discussed earlier, and has a finite series expansion: 



oo 




(4.76) 



n=0 



where a ^ 0. The indicial equation for m is 



m(m - l)(3m + 5) (3m - 4) = . 



(4.77) 



a k = Dx * , 

Z k = Dx-i(\2 + x 2 ) , 



(4.78) 
(4.79) 
(4.80) 
(4.81) 



(4.82) 



where D is a constant. The non-local variables vanish. 
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The first two terms of the mode with m = are 
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(4.83) 
(4.84) 
(4.85) 
(4.86) 
(4.87) 
(4.88) 
(4.89) 



where a\ is a constant. The form of this solution is similar to the adiabatic growing 
mode of general relativity. 

The mode corresponding to m — 1 is 



$k = a 2 

&k = «2 

Zk = a 2 
A fc = a 2 
qt = a 2 
T fc = a 2 
Qk = a 2 
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(4.90) 
(4.91) 
(4.92) 
(4.93) 
(4.94) 
(4.95) 
(4.96) 



with a 2 a constant. This mode has no analogue in general relativity, except that <3>fc 
is a growing mode. As t — > there are non-zero but compensating contributions 
to the effective peculiar velocity g* ot /p tot from the matter and the non-local en- 



84 



4. BRANEWORLD COSMOLOGY II: 
INITIAL CONDITIONS AND CMB ANISOTROPIES 



ergy fluxes. The contributions of these components to the fractional total density 
perturbation aD a p tot / p tot vanish as t — > 0. 

The mode corresponding to m = — | is singular as t — > (it is a decaying mode 
with a 3 a constant): 



= a 3 x~l (l - j^x 2 \ , (4.97) 

a k = a 3 x-i (\ + ^rr 2 ^) , (4.98) 

(U 4 1217 io \ 

* = ° 3 («" " 1590435^) ' (4 "» 

( 4 10 4 is \ _ , 

= 03 a; 3 x 3 4.101 

yfc 3 \4455 1301265 / ' V ; 

1 7 7 13 ^ 0-3 1 1 7 2 



T, = as ( - jjg*. + — *t j = _„ ^ + , (4.102) 

a = 03 1 -h xi + is; 1 *) = I* 1 (- 1 + s* 2 ) • (4103) 



A similar mode is found in general relativity but there the decay of is more 
rapid ($£ oc a; -3 ) on large scales. The solution (4.103) describes an isocurvature 
velocity mode 1 where the early time matter and non-local (Weyl) components 
have equal but opposite peculiar velocities in the CDM frame. The existence of 
such isocurvature modes was predicted in [77] and [118] for large-scale density 
perturbations. 



'We note that we have identified this mode incorrectly in [123] since the variable Q is rescaled 
byx 4 / 3 . 
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Finally, for m = |, with a 4 a constant, we have 

^k = a iX l (l-^La?^ , (4.104) 

-% X+ Am x ) ' (4105) 

Z k = a^ (^--^- X A , (4.106) 
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A k = a A x^ ( -4 + 4^ ) > ( 4 - 107 ) 



9* = «4^3 ( - ) , (4.108) 



4 



T(x) = a 4 *i ( ? + ±-x> ) = ( 1 + Is* ) , (4.109) 
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2 28 y 2 V 42 

29 3 \ 3a 4 4 29 



G(s) = -_«J=_x. ^-2106* J- (4-110) 

In this mode the universe asymptotes to an FRW (brane) model in the past as 
t — > 0. Note that this requires careful cancellation between aD a p tot / p tot and g* ot 
to avoid a singularity in the gravitational potential $^ (which would diverge as 
a; -2 / 3 without such cancellation). Like the velocity isocurvature mode (m = 1) 
discussed above, this mode has no analogue in general relativity. We identify the 
solution (4.109) to be the non-local isocurvature density mode. The pair of non- 
local velocity and density isocurvature modes (4.103) and (4.109) is analogous to 
the neutrino density and velocity isocurvature modes (see Eqs. (4.38) and (4.39)). 



§4.2.2 The Energy Frame 

In this section, we present a complete set of evolution equations for the total 
matter variables in the matter energy frame, q a = 0. Note that the four-velocity of 
the energy frame is not necessarily a timelike eigenvector of the Einstein tensor in 
the presence of the non-local braneworld corrections to the effective stress-energy 
tensor. We assume that the matter is radiation dominated, the non-local energy 
density vanishes in the background, and we ignore local anisotropic stresses. We 
also assume that the baryons and CDM make a negligible contribution to the 
fractional gradient in the total matter energy density and to the energy flux, thus 
excluding the CDM and baryon isocurvature modes. We also give the evolution 
equations for the non-local density gradient and energy flux in the matter energy 
frame. 
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Denoting the variables in the energy frame by an tilde on the quantity, for 
example, A a , the relevant equations for scalar perturbations are 

A a = leA a -|z , (4.111) 
z a = -§e*. - \d>K - %r a - (i + f ) , (4.ii2) 

f a = --D 2 Q a , (4.113) 
P 

ka = ~l<dQ a - £r a - ^A a - D b V ab . (4.114) 

Solutions of these equations are related to those in the CDM frame (§4.2.1) by 
linearising the frame transformations given in [140]. If the CDM projected velocity 
is v$ in the energy frame, the variables in the CDM frame are given by 

A a = A a -^aQv^ , (4.115) 

Z a = Z a + l -D a D%^ - (l + P-) , (4.116) 
T a = T a , (4.117) 

Qa = Qa, (4.118) 

Qa = ~pv a c) , (4.H9) 

where we have used U = in the background. Note that the quantities T a and 
Q a are frame-invariant, hence we can assume the same form for these quantities 
for the subsequent calculations in the energy frame like the CDM frame. 
The CDM peculiar velocity evolves in the energy frame according to 

W = 4©^ c) + ^Aa- (4.120) 

Energy frame solutions in the High-energy Regime 

In the high energy limit, p/X 1, using the variable x defined in Eq. (4.47), 
we obtain the following four equations: 



9xA k = 3A fe - 4xZ k , (4.121) 
3x 2 Z' k = -xZ k + (^-- 15] A k - 36T fc , (4.122) 



v 4 

3xf' k =Af k + xQ k , (4.123) 
3xQ k =4Q k -l (t k + \K k - 2V k ) . (4.124) 
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Setting Vk = 0, we decouple the above equations to find a 4th order differential 
equation, 



9 x ^K k + 3^9^ + /V _ 5^\ 9^ + f x + <9A* 



16 <9a; 4 4 <9x 3 \ 24 4 / & 2 V 36/ dx 

( 7 x 2 x \ ~ 
+ -2 + 36 + 1296 ) Ak = °- 



(4.125) 



Using eq. (4.76), we obtain the indical equation: 

(m - 2)(3m - 4) (3m - 7)(1 + m) = . (4.126) 

We proceed to substitute each value of m back in (4.125) and solve the recurrence 
relations for {b n }. The original set of equations (4.121)-(4.124) has four degrees 
of freedom, and we note that the peculiar CDM isocurvature mode vanishes in the 
energy frame. In addition, we compute the spatial gradient of the 3-curvature, r] k 
from eq. (3.82). 

The mode corresponding tom = 2 are 



* , ( 15 11 3 

Zu — Oi x H x 

k 1 V 4 200 

-Y" h f 4 3 6 

Tu — —x H x 

k 180 V 770 

Sfc= 45 ("* + 440* 
, ^27 7 2 



The mode corresponding to m = | are 



3 

~ , / 9 i 9 7 
Z k = h (--** + — *i 

- 7b 2 if 1 9 

T » = -2f" ( 1 + 42 X 

= 6 2 4 / 29 o 

Qu = is x x 

24 V 2106 



4 



, ■ / 1 19 2 

Vk = ^2^ 3 X 

' 5 4 14040 



4.127) 
4.128) 
4.129) 
4.130) 
4.131) 

4.132) 
4.133) 
4.134) 
4.135) 
4.136) 
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The mode corresponding to m — | are 

A k = /,,,-, (l - ^ , , 

~ , / 9 4 1 10 

^ = & 3 + ^ 3 

- 56 3 f 1 7 13 

Ti. = — a; 3 H x 3 

fc 24 V 270 

- 56 3 / 4 7 10 
Qt- = -is -| x 3 

k 8 V 90 

4 56 

The mode corresponding to m = — 1 are 



,i/45 9 2 

7] k = b 3 X'i — - —x 



K k = hx- 1 (l + ^x 2 ^j , 



3 1 

x 2 4 



Tfc = -18 ^ + 30 X 

k 18 V 6 
h ( 1 5 2 



4.137) 
4.138) 
4.139) 
4.140) 
4.141) 



4.142) 
4.143) 
4.144) 
4.145) 
4.146) 



with 6i,6 2 , b 3 and 64 being constants. 

From the solutions (Eqs. (4.127)-(4.145)) in the high energy frame, we observe 
that powers and coefficients of the solutions for T and Q match the ones in the 
CDM frame with constants of integration (up to some pre-factor). The solutions 
are checked by transforming them to the CDM using Eqs. (4.115)-(4.119). Eqs. 
(4.134) and (4.140) correspond to the pair of non-local isocurvature density and 
velocity modes found earlier in the CDM frame. From Eq. (4.131), we find that 
the spatial gradient of the 3-curvature for the adiabatic mode is a constant on 
large scales as x — > 0. We have arrived at the same result found in [21, 77]. With 
a more detailed analysis, we can characterize the m = 2 and m — — 1 modes as 
the growing and decaying modes that are similar to the case in general relativity. 
We found that the curvature is conserved for the GR-like modes for U = 0, but 
not for the other two additional braneworld modes. 
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§4.2.3 Discussion on the Scalar Modes 

Regarding the imprint of braneworld effects on the CMB, we note several pos- 
sible sources. Once the universe enters the low-energy regime the dynamics of 
the perturbations are essentially general relativistic in the absence of non-local 
anisotropic stress (see §4.2.1). If V a b really were zero, the only imprints of the 
braneworld on the CMB could arise from modifications to the power spectrum 
(and cross correlations) between the various low-energy modes. Since there are 
two additional isocurvature modes in the low-energy universe due to braneworld 
effects, it need not be the case that adiabatic fluctuations produced during high- 
energy (single-field) inflation give rise to a low-energy universe dominated by the 
growing, adiabatic, general- relativistic mode (See 4.1 for a history of the scalar 
modes from the high to low energy regime). Hence there exist the possibility of 
exciting the low-energy isocurvature (brane) modes from plausible fluctuations in 
the high-energy regime. 

However, if we consider linear scales at last scattering which project onto an- 
gular scales corresponding to I ~ 200 and h = 0.75 in a = and A = 
universe, assuming the standard nucleosynthesis constraint [refer to eq.(3.14) in 
the previous chapter] and ignoring the contributions of the dark radiation term, 
we find these modes enter the Hubble radius at approximately 10 4 . This large 
redshift would indicate that the modes are outside the Hubble radius. If this is 
the case, the matching from the high energy regime (brane-world) to the low en- 
ergy regime (general relativity) would render the modes to be super-Hubble. An 
immediate consequence is that the CMB anisotropics would be insensitive to the 
brane tension, A, i.e. the high energy effects. 

Similarly in practice, we do not expect V a b = 0. In this case the non-local 
anisotropic stress provides additional driving terms to the dynamics of the fluctu- 
ations, and we can expect a significant manifestation of five-dimensional Kaluza- 
Klein effects on the CMB scalar anisotropies. Barrow and Maartens have shown 
in [9] that the 5-dimensional KK graviton stresses can slow down the decay of 
the shear anisotropy on the brane to observable levels. They found that with a 
suitable approximation for V a b that the initial shear to the Hubble distortion of 
~ IVPVLqHI at 5D Planck time would allow the large-angle CMB signal to be a 
relic of purely KK effects. These possibilities may reveal new braneworld imprints 
that would be more realistic for observational testing. 



4.3 Braneworld Tensor Anisotropies in the CMB 



In this section, we investigate the dynamics of tensor perturbations in braneworlds. 
Recalling Eqs. (3.111) and (3.112) from Chapter 3, we rewrite these equations in 
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Figure 4.1: A figure which depict history of the different modes for the 
braneworlds. The KK mode, the CDM and the neutrino isocurvature modes 
will trail from the high energy to the low energy regime, after inflation. In the 
high energy regime, where p/X S> 1, there are 4 degenerate modes in the case 
of the braneworld. When the instantaneous transition from high to low energy 
regime (p/X <C 1) at p/X = 1, the modes will be split up from four modes 
to two modes which correspond to the case in general relativity and the other 
remaining two modes which are the non-local isocurvature modes. The two re- 
maining non-local isocurvature modes from the high energies are redshifted by 
approximately 10 4 if there is no dark energy term. Hence the only effect on the 
CMBR anisotropies is induced only by the KK mode. 



conformal time: 



k k 2 ^ 

- 2 K + na k ) + -E k - ^-pn k = k 2 {2 - 3 7 )^7r fe + ^-pvr* , (4.147) 
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k 2 , / k 2 3K 

- (E' k + HE k ) -k^- + — - Y 7PJ <r k + Y P< - Y {Sl ~ 1)Hp7Vk 

= -^{2k 1P 2 a k - (3 7 - 2)pV fc - [3 7 ' - (3 7 - 2)(6 7 - l)H]p 2 n k } (4-148) 

2 k 2 \ , I 

- -kn 2 p*a k - — y P n* k + (1 - 37)Wp7r^J , 

where the non-local energy density is defined as in [9]: 

P* = ^U. (4.149) 

A prime denotes d/dr, 7i = a' /a, and the (non-constant) parameter 7 is defined 
in Eq. (3.74). Equations (4.147) and (4.148), which have all the braneworld terms 
on the right-hand sides, determine the tensor anisotropies in the CMB, once 7i k 
and ix* k are given. 

§4.3.1 A Local Approximation to the Nonlocal Anisotropic 
Stress 

The general solution for the non-local anisotropic stress will be of the form 

tt£(t) oc J df g(r,f)E[n k ,a k }\ f , (4.150) 

where Q is a retarded Green's function evaluated on the brane. The functional F 
is known in the case of a Minkowski background [167], but not in the cosmological 
case. (An equivalent integro-differential formulation of the problem is given in 
[148]; see also [178].) Once Q and F are determined, Eq. (4.150) can in principle 
be incorporated into a modified version of Boltzmann codes such as CAMB [126] 
or CMBFAST [175]. It remains a major task of braneworld cosmological pertur- 
bation theory to find this solution, or its equivalent forms in other formalisms. 
In the meanwhile, in order to make progress towards understanding braneworld 
signatures on tensor CMB anisotropies, we can consider approximations to the 
solution. 

The non-local nature of 11 k1 as reflected in Eq. (4.150), is fundamental, but is 
also the source of the great complexity of the problem. The lowest level approx- 
imation to -rcl is local. Despite removing the key aspect of the KK anisotropic 
stress, we can get a feel for its influence on the CMB if we capture at least part of 
its qualitative properties. The key qualitative feature is that inhomogeneity and 
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anisotropy on the brane are a source for KK modes in the bulk which "backre- 
act" [135] or "feed back" [147], onto the brane. The transverse traceless part of 
inhomogeneity and anisotropy on the brane is given by the transverse traceless 
anisotropic stresses in the geometry, i.e. by the matter anisotropic stress n a b and 
the shear anisotropy a a b. The radiation and neutrino anisotropic stresses are in 
turn sourced by the shear to lowest order (neglecting the role of the octupole and 
higher Legendre moments). 

Thus the simplest local approximation which reflects the essential qualitative 
feature of the spin-2 KK modes is 

K 2 Tr* ab = -(Ha ab , C' = 0, (4.151) 

where ( is a dimensionless KK parameter, with ( = corresponding to no KK 
effects on the brane, and ( = = A^ 1 giving the general relativity limit. [Note that 
for tensor perturbations, where there is no freedom over the choice of frame (i.e. 
u a ), there is no gauge ambiguity in Eq. (4.151). However, for scalar or vector per- 
turbations, this relation could only hold in one frame, since 7r* b is frame-invariant 
in linear theory while a a b is not.] 

The approximation in Eq. (4.151) has the qualitative form of a shear viscosity, 
which suggests that KK effects lead to a damping of tensor anisotropies. This is 
indeed consistent with the conversion of part of the zero-mode at Hubble re-entry 
into massive KK modes [117, 76]. The conversion may be understood equivalently 
as the emission of KK gravitons into the bulk, and leads to a loss of energy in the 
4D graviton modes on the brane, i.e. to an effective damping. The approximation 
in Eq. (4.151) therefore also incorporates this key feature qualitatively. 

With this first approximation, we can close the system of equations on the 
brane by adding the equation 

K 2 pn* k = -CH^a k . (4.152) 

cr 

We will also assume K = = p* in the background. The parameter ( (together 
with the brane tension A) then controls braneworld effects on the tensor CMB 
anisotropies in this simplest approximation. 



Initial Conditions 



Ignoring the photon anisotropic stress (i.e. n ab = 0), we differentiate Eq. 
(4.147) and decouple it with Eq. (4.148) and the modified Raychaudhuri equation 
[see Eq. (3.59) of Chapter 3]. Using the variable Uk = a 1+ ^ 2 , the shear satisfies 
the following equation of motion 



u'! + 



k 2 + 2K 



-l-C/2 



u k = 0, 



(4.153) 
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where we have used Eq. (4.152). In flat models (K = 0) on large scales there is a 
decaying solution a k oc a~^ 2+c -\ Since Eq. (4.153) contains no first derivative term 
the Wronskian is conserved. On large scales we can use the solution a k oc a~( 2+ ® 
to write the conserved Wronskian as W = a' k + (2 + ()Ha k . (The Wronskian 
vanishes in the decaying mode.) Integrating gives the following two independent 
solutions on large scales in flat models: 



CTk 



B k a-V + V f T dfa(f) 2 +t, 



(4.154) 



where A k and B k are constants of integration. If we let ( — > 0, we recover the 
results in [125] for the general relativity case. 

The conserved Wronskian is proportional to the metric perturbation variable, 
H T , characterising the amplitude of 4D gravitational waves. In flat models H T is 
related to the covariant variables quite generally by 



Tk 



<^k 
k 



2E k 



(4.155) 



Ignoring photon anisotropic stress, we can eliminate the electric part of the Weyl 
tensor via the shear propagation equation (4.147) to find kHxk = ~ (J ' k — (C + 
2)7ia k = —W. The fact that Ht is conserved on large scales in flat models in the 
absence of photon anisotropic stress can also be seen directly from its propagation 
equation, 

H'r k + (2 + QHH' Tk + k 2 H Tk = . (4. 156) 

We can solve Eq. (4.153) on all scales in the high-energy (p> A and a oc r 1 / 3 ) 
and low-energy (/) € A and a oc r) radiation-dominated regimes, and during 
matter-domination (a oc r 2 ). The solutions are 

u k (r) = Vkr di Ji( 5+ ^(A;t) + d2Yi^ 5+ ^(kr) (high energy radiation), 



u k (r 



kr 



(4.157) 

^3^i(3+C)(^ T ) + ^4^i(3+c)^ r ) (l° w ener gy radiation), (4.158) 



u k {r) = Vkr d^Jb + ^{kr) + dgYs^Arr) 



(matter domination), (4.159) 



where di are integration constants. The solutions for the electric part of the Weyl 
tensor can be found from Eq. (4.147). For modes of cosmological interest the 
wavelength is well outside the Hubble radius at the transition from the high- 
energy regime to the low-energy. It follows that the regular solution (labelled by 
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di) in the high-energy regime will only excite the regular solution (d 3 ) in the low- 
energy radiation-dominated era. Performing a series expansion, we arrive at the 
appropriate initial conditions for large-scale modes in the low-energy radiation era: 

g " = 1 -2g^) + 8(3 + ( S + +OK * rn <4160) 
g > = -3TC + 2(3 + 0^5 + + ° [( * r) ' 1 - <4161) 

£t= 2(3T0"4(3 + 0(5 + +O[(lr)1 - <4 - 162) 

In the limit £ — * 0, we recover the general relativity results [31]. 

For modes that are super-Hubble at matter-radiation equality (i.e. kr cq <C 
1), the above solution joins smoothly onto the regular solution labelled by d 5 in 
Eq. (4.159). For kr cq 3> 1, the shear during matter domination takes the form 

a k = -2f +( T(§ + C)(kry^J lH (kr). (4.163) 

In the opposite limit, the wavelength is well inside the Hubble radius at matter- 
radiation equality. The asymptotic form of the shear in matter domination is 
then 



Gk „ F[ ^ C)] l+C/2 (krp 1+ ^ sin(kr - ttC/4). (4.164) 

We use the initial conditions, Eqs. (4.160)-(4.162), in a modified version of 
the CAMB code to obtain the tensor temperature and polarization power spec- 
tra. The temperature and electric polarization spectra are shown in Figs. 4.2 and 
4.3 for a scale-invariant initial power spectrum. The normalisation is set by the 
initial power in the gravity wave background. Figs 4.2 and 4.3, together with 
Eqs. (4.157)-(4.162), are the main result of this work, and we now discuss the 
physical conclusions following from these results. 



§4.3.2 Discussion on the Tensor Power Spectra 

As expected, we find that the power spectra are insensitive to high-energy 
effects, i.e. effectively independent of the brane tension A : the ( = curve in 
Fig. 1 is indistinguishable from that of the general relativity model (both power 
spectra are identical at the resolution of the plot). For the computations, we have 
used the lowest value of the brane tension A = (100 GeV) 4 , consistent with the 
tests of Newton's law. 
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There are three notable effects visible in Figs. 1 and 2, from our approximate 
model of the KK stress: (i) the power on large scales reduces with increasing 
KK parameter (; (ii) features in the spectrum shift to smaller angular scales with 
increasing (; and (iii) the power falls off more rapidly on small scales as ( in- 
creases. Neglecting scattering effects, the shear is the only source of linear tensor 
anisotropies (see e.g. [31]). For 1 Z < 60 the dominant modes to contribute to 
the temperature C/s are those whose wavelengths subtend an angle ~ 1/7 when 
the shear first peaks (around the time of Hubble crossing). The small suppression 
in the C/s on large scales with increasing ( arises from the reduction in the peak 
amplitude of the shear at Hubble entry [see Eq. (4.163)], qualitatively interpreted 
as the loss of energy in the 4D graviton modes to 5D KK modes. 

Increasing ( also has the effect of adding a small positive phase shift to the 
oscillations in the shear on sub-Hubble scales, as shown e.g. by Eq. (4.164). The 
delay in the time at which the shear first peaks leads to a small increase in the 
maximum I for which l(l+l)Ci is approximately constant, as is apparent in Fig. 4.2. 
The phase shift of the subsequent peaks in the shear has the effect of shifting the 
peaks in the tensor CjS to the right. For I > 60 the main contribution to the 
tensor anisotropies at a given scale is localized near last scattering and comes 
from modes with wavenumber k ~ Z/t , where r is the present conformal time. 
On these scales the gravity waves have already entered the Hubble radius at last 
scattering. Such modes are undergoing adiabatic damping by the expansion and 
this results in the sharp decrease in the anisotropies on small scales. Increasing the 
KK parameter ( effectively produces more adiabatic damping and hence a sharper 
fall off of power. The transition to a slower fall off in the C/s at I ~ 200 is due to the 
weaker dependence of the amplitude of the shear on wavenumber at last scattering 
for modes that have entered the Hubble radius during radiation domination [179]. 
[The asymptotic expansion of Eq. (4.163) gives the shear amplitude oc k~( 2+ ^ 
at fixed r, whereas for modes that were sub-Hubble at matter-radiation equality 
Eq. (4.164) gives the amplitude oc £r (1+c/2) .] 

Similar comments apply to the tensor electric polarization Cf , shown in Fig. 4.3. 
As with the temperature anisotropies, we see the same shifting of features to the 
right and increase in damping on small scales. Since polarization is only generated 
at last scattering (except for the feature at very low / that arises from scattering 
at reionization, with an assumed optical depth t c = 0.03), the large-scale polar- 
ization is suppressed, since the shear (and hence the temperature quadrupole at 
last scattering) is small for super-Hubble modes. In matter domination the large- 
scale shear is = —kr/(5 + 2£); the reduction in the magnitude of the shear 
with increasing KK parameter ( is clearly visible in the large-angle polarization. 
The braneworld modification to the tensor magnetic polarization Cf has the same 
qualitative features as in the electric case (see Fig. 4.4). 
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Figure 4.2: The temperature power spectrum for tensor perturbations in 
braneworld models using the approximation in Eq. (4451), with £ the dimen- 
sionless KK parameter. Models are shown with £ = 0.0, 0.1, 0.25, 0.5 and 
1.0. The initial tensor power spectrum is scale invariant and we have adopted 
an absolute normalisation to the power in the primordial gravity wave back- 
ground. The background cosmology is the spatially flat ACDM (concordance) 
model with density parameters Oj, = 0.035, f2 c = 0.315, £1\ = 0.65, no massive 
neutrinos, and the Hubble constant Hq = 65 km s Mpc -1 . 




Figure 4.3: The electric polarization power spectrum for tensor perturbations 
for the same braneworld models as in Fig. 4.2. 
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Figure 4.4: The magnetic polarization power spectrum for tensor perturbations 
for the same braneworld models as in Fig. 4.2. 



Part III 
Changing Global Symmetry 



Chapter 5 

Changing Global Symmetry 



Some say the world will end in fire 
some say in ice. 
From what I've tasted of desire 
I hold with those who favor fire. 

- Robert Frost 

This is the way the world ends, 
Not with a bang but a whimper. 

- T.S. Eliot 

In this chapter, we examine the changing global symmetry of the Einstein-de 
Sitter Universe and see how it could lead to a possible phase transition in the future 
as in reference [10]. The largest galaxy cluster in a flat Einstein-de Sitter universe 
may grow indefinitely to encompass most of space after an extremely long time. 
We derive a general relativistic metric for a very large, bounded nearly isothermal 
cluster of galaxies which is embedded in such a universe. The embedding is done 
by means of a "Schwarzs child membrane" . Pressure is important, unlike previous 
Tolman-Bondi models of inhomogeneities. The cluster's expansion, represented by 
a sequence of models, alters the average global symmetry of an increasing volume of 
space-time. Initially the metric is homogeneous and isotropic, having translational 
and rotational symmetry around every point. As the metric evolves, it eventually 
loses its translational symmetry throughout larger regions and retains rotational 
symmetry around only one point: the centre of the cluster becomes the centre of 
the universe. Our sequence of hybrid models transforms between the Einstein-de 
Sitter and the isothermal cluster limits and illustrates how a changing equation of 
state of matter can alter the global symmetry of space-time. 
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§5.1 Introduction 

It is well known in standard texts [40, 156, 189] that one can have a changing 
equation of state in a single metric. For example, one can study the behaviour 
of the FRW metric in the radiation-dominated and the matter-dominated era 
which correspond to different equations of state. In these situations, the equation 
of state does not influence the global symmetry of the metric, and there is a 
smooth transition from the radiation-dominated to the matter-dominated era. In 
other previous inhomogeneous models, the equation of state has zero pressure 
[15, 58, 112, 113, 122]. Therefore, it remains an interesting fundamental question 
to ask whether such a dynamical symmetry breaking process can be modelled in 
general relativity. To paraphrase this question in a more concrete form, can one 
find a metric to describe a changing global symmetry, with a changing equation of 
state that will alter the form of the metric? Here, the basic effect of the pressure 
on the metric arises not from its contribution to the energy density or to the 
expansion, but from its support of growing inhomogeneity. 

At first sight, it may not seem possible to solve the Einstein field equations 
directly for a changing equation of state which induces symmetry changes in the 
metric. Even the form of such an equation of state has not been examined previ- 
ously. Hence, to answer the above question, we need to adopt a different approach 
to find a more general solution. Our method is to find a continuous hybrid metric 
which is consistent with the known limit-solutions. The hybrid solution is parame- 
terised by two constants and applies only to the quasi-static phases of the universe 
when it is dominated by one gigantic isothermal cluster. As a second cluster of 
comparable size enters the horizon the hybrid description breaks down. The clus- 
ters merge and eventually form a new, larger isothermal sphere. After relaxation 
of this larger system the hybrid solution is valid again with changed values of its 
characteristic constants. 

This allows us to separate the problem into two parts. First, we find the 
parameterised hybrid solution which describes those periods in the evolution of 
the universe, which are dominated by a large isothermal core. We then describe 
the temporal evolution as a discrete sequence of such stages i.e. as a sequence of 
hybrid solutions with changed parameterisation. 

Due to the different equations of state, it is not possible to find direct matching 
conditions for the generalised isothermal metric and the FRW metric. The reason 
can be understood physically. Clustering induces a change from a universe with 
zero pressure to a universe with non-zero cosmologically significant pressure. There 
is a symmetry change from the spherical, homogeneous FRW models, which have 
no centre to the isothermal model, which singles out the point with highest density 
as its centre. 



§5.2. THE ISOTHERMAL UNIVERSE 
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§5.2 The Isothermal Universe 

§5.2.1 The Isothermal Metric 

We begin with a brief summary of the isothermal universe found by Saslaw et 
al[172]. The starting point is to consider an asymptotically static universe that 
satisfies the general static, spherically symmetric line element: 

d S 2 = _ e Hr) dt 2 + e X(r) dr 2 + ^^2 + ^2 M(j) 2y ^ 

The Einstein field equations are given by 

G ab = R ab - X -Rg ah = MT ab) (5.2) 
where the energy-momentum tensor of a perfect fluid is 

Tab = (p + P)u a u b + Pg ab . (5.3) 

To derive the solution for the isothermal universe, an equation of state is required 
to determine the metric. The isothermal equation of state is characterized by a 
pressure gradient which balances the mutual self-gravity of its constituent particles 
i.e. idealized point galaxies. The dispersion of the particles' peculiar velocities is 
independent of position. Therefore the equation of state is: 

P = ap, (5.4) 

where a is a constant for the whole spacetime and < a < 2/3. 

There are two assumptions in obtaining the solution of the isothermal sphere. 
Firstly, we assume e A<r - ) to be independent of time. This has the geometrical in- 
terpretation that the radial coordinate of the metric scales radially as a constant. 
Secondly, we assume the density distribution of a finite isothermal sphere (which 
realistically would be bounded by an external pressure or tidal disruption) is de- 
scribed by the relation p oc r~ 2 at large radii, as was first shown by Emden [67] for 
the Newtonian case and subsequently extended to the relativistic case by Chan- 
drasekhar [36]. To obtain the cosmological metric describing an infinite isothermal 
universe one has to make the simplifying approximation that p oc r~ 2 throughout 
the entire sphere. Solving the Einstein field equations (5.2) using the line element 
(5.1) together with Eqs. (5.3) and (5.4) and the density distribution given above, 
gives the following relations: 

e u = Ar^, (5.5) 
4a 
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The line element for the isothermal universe is therefore stated as follows: 



ds 2 = -Ar^dt 2 + ( 1 + . — , 9 ) dr 2 + r 2 (d6 2 + sin 2 6d(j) 2 ), (5.8) 
V (! + «)/ 

where A is an arbitrary constant. 

Here we find that the transformation f = r u casts the metric in Eq. (5.8) into 
its isotropic form 

rfs 2 = _ Af V dt 2 + f 2W [ df 2 + f 2 dQ 2} ? ( ^ g) 

where the functions U, V and W are defined as 

4a 

l ' = f/(TT^)' < 5 ' n) 

W = ^ - I- (5.12) 

For completeness, we also state the conformal form of Eq. (5.8), which relates to 
minimally curved spaces, as found by Dadhich [42]: 

ds 2 = [-dt 2 + k\dr 2 + r 2 (d6 2 + sin 2 Odcj) 2 )] . (5.13) 
where e A = k\ and e v = r~ 2m . Note that k\ and m are constants relating to a. 



4a 



§5.2.2 Properties and Implications of the Isothermal solu- 
tion with a Global Phase Transition of the Universe 

Before we proceed further, we note a few interesting features with the isother- 
mal solution, (5.8), found by Saslaw et al in [172]. First of all, the isothermal 
universe in their infinite isothermal sphere without a boundary. The 

solution is not a realistic one which we could use to model the ultimate phase 
transition of the universe. However, even with that limitation, from examining 
the physical interpretation of the unbounded isothermal sphere, we could infer 
why a global phase transition would eventually happen. 

From an extensive study in the dynamics of galaxy clustering (which was well 
discussed in the references [169, 170, 171]), there is a discussion on whether the 
many-body clustering of galaxies in our expanding Universe would be represented 
as a phase transition. The conclusion is reached in [170] that galaxy clustering 
lacks the near simultaneity of the transition over all scales, even though it has some 
basic features of such a phase transition. In Einstein-Friedmann cosmologies with 
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Qo < 1, and in related cosmologies with > 0, the universe would eventually 
expand so rapidly relative to the gravitational clustering timescales such that the 
pattern of galaxies would freeze out on large scales, and correlations cease to grow, 
and the attempt at a second order phase transition would die away before it could 
ever happen. 

However, in the case of the flat Qq = 1, Q\ = 0, Einstein-de Sitter cosmology, 
there is a possibility for this phase transition to occur. In present day observational 
cosmology, in the viewpoint of the data from the observed fluctuations in the 
cosmic microwave background anisotropics [14, 115], it seems very likely that our 
Universe is nearly flat. However, whether there is a cosmo logical constant, is still 
being debated, but if there is, the consequence of a phase transition would have 
been different. 

One would ponder how the phase transition would enter in the unbounded 
isothermal sphere scenario. The first thing is to recognise that there is a change of 
fundamental symmetry between the Einstein-de Sitter model and the isothermal 
model. We recall the form of the Friedmann- Robertson- Walker metric (Eq.(l.l) 
in Chapter 1). 

To match Eqs.(l.l) and (5.8), we need to show that the pressure must be 
discontinous at the junction across the hypersurface where r is a constant. Com- 
parison of the two metrics would show that for a > 0, there is no such hypersurface 
where this match occurs. Therefore this change must be continuous. Since there 
is a time dependence in the global expansion a(t), and the isothermal metric has 
no time dependence, it suggests that in the limit that t — > oo, the transition would 
take place everywhere. Such an event would correspond to a phase transition. Of 
course, a more complicated metric that allows growing inhomogeneities would be 
required to follow this transition. 

This has set the scene and motivation for this chapter on the theme of changing 
global symmetry in our thesis on alternative cosmologies. In [10], we have found 
and studied a hybrid model which corresponds to the global symmetry change 
mentioned above. In the next section, we will subsequently derive a generalized 
form of their metric to include a bounded isothermal sphere using an approach 
by Tolman [184]. We review and discuss the properties of the bounded isother- 
mal sphere. This is followed by a discussion of the approximations involved in 
the general solution (§5.3) and the thermodynamic consequences implied by the 
modified equation of state (§5.4). In §5.5, we first match the bounded isothermal 
sphere to a Schwarzschild background, and then embed the combined metric into 
the expanding Einstein-de Sitter universe, adapting earlier techniques (McVittie 
[145, 146] and Hogan [89]). In §5.6, we describe the evolution of the universe as 
a sequence of static stages. Each stage is described by a hybrid metric with char- 
acteristic parameterisation. Finally, in §5.7, we discuss the implications of such 
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a cosmological solution, linking the dynamical breaking of symmetry to a phase 
transition. In the same light, we mention an additional role for the cosmological 
constant which counteracts the formation of the isothermal sphere. 



§5.2.3 Generalization of the Isothermal Metric to Finite 
Spheres 

The basic problem in the search for a solution by direct matching of the isother- 
mal metric and the FRW metric, lies in the finite pressure-mismatch between the 
FRW models, which are based on the assumption of zero pressure, whereas our 
FRW fitting uses zero pressure, but the isothermal sphere could have a nonzero 
uniform pressure. To circumvent this difficulty, one has to generalise the infinite 
isothermal sphere in [172] to a finite sphere, where the pressure drops to zero at the 
boundary. Then it is possible to match this finite sphere solution to the expanding 
space surrounding it. 

To generalise the isothermal solution we employed the technique used by Tol- 
man [184]. By starting with the general static, spherically symmetric line element 
described by Eq. (5.1) one arrives at the following field equations for the perfect 
fluid case: 

dr ( r 2 ) dr ( 2r ) dr ( 2r ) ' ^ ^ 



87rp = e- A (- + V) -4 (5.15) 



ny ny 2 / ^"2 



8nP = e- x [--\)+\. (5.16) 



ty ty 2 / ,y*2 

We assume e~ A = constant in order to obtain a simple generalization of the isother- 
mal equation of state, so that a very large central cluster can have a finite boundary. 
This allows us to integrate Eq. (5.14), and leads to the corresponding solution: 

e u = (Cr 1 -" 1 - Dr 1+n ) 2 , (5.17) 
e A = 2-n 2 , (5.18) 

1 — n 2 1 

2 — n A r z 



o n 2 r 2 Q _ £) r 2n 



where n, C and D are arbitrary constants of integration. 
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To relate Tolman's solution to the isothermal metric, one makes the following 
identification 

\ 4a o . 

1 + (TTa) 5 = 2 ~ ™ ■ 

where 

n = ±^^. (5.22) 
1 + a v 7 

By taking the positive root* of the above relation the isothermal metric is gener- 
alised to the form 

e v = (Cr^ - Dr&y , (5.23) 

4a 1 . 

4a + (1 + a) 2 r 2 

2(l-a) 

4a 2 1 C-\Dr^+^ 

8ttP = -r — a m , . 5.26 

4a + (1 + a) 2 r 2 c _ Dr ^ V ; 

D is interpreted as a deviation term for the isothermal sphere. The original isother- 
mal metric without a boundary is included in this form in the limit D — > 0. 

The specific case for n — 1/2 was discussed by Tolman [184]. In the limit 
D/C — > 0, the ratio of pressure to density (a) approaches one-third throughout and 
the sphere grows without limit. Physically, it represents the blackbody radiation 
solution in the Oppenheimer-Volkoff analysis. In the limit of large p, the equation 
of state for the sphere goes over to the approximate form, p — 3P oc p 1 ^ 2 , as for a 
highly compressed Fermi gas. 

In the generalised form of the bounded isothermal metric, we still require the 
fixed density distribution p oc r -2 , but allow the pressure-density relation to devi- 
ate from the perfect isothermal equation of state. Hence, we deduce the following 
generalised isothermal equation of state: 

, 2(l-q) 

P C-\Dr^r 

- = « — - o ' ( 5 - 27 ) 

P C-Dr^r 

where a is defined as the constant of proportionality between pressure and den- 
sity at the centre of the sphere. Note that the above equation of state allows the 
pressure to drop to zero at a finite radial distance, rv Hence, our metric incorpo- 
rates the feature of describing a finite sphere, each shell of which is approximately 

""Choosing the negative root gives the same solutions with C and D interchanged. This reflects 
an implicit symmetry of the Tolman solution. 
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isothermal. It is the requirement that the isothermal solution has a boundary, 
which introduces deviations from the perfect isothermal equation of state. A re- 
lated class of spherically symmetric solutions of Einstein's equation which consist 
of a perfect fluid with uniform density but non-uniform pressure, was found by 
Wesson and Ponce de Leon [190, 191]. In their solution, the equation of state is 
similar to Eq. (5.27) except that their equation of state for the metric is both 
dependent on time, t and radial position, r. These solutions are recognised to be 
scale free and self-similar and were discussed in detail by Carr [25, 26]. The nature 
of equation (5.27) is that at any given radius and for any given value of a, the 
pressure P is proportional to p, but the constant of proportionality varies with 
position (see §5.4 for a discussion). We denote this as a generalised isothermal 
sphere. Formally, this is a result of the finite boundary. Physically, this could re- 
sult from incomplete dynamical relaxation in the outer regions of the cluster where 
the densities are lower and the relaxation timescales longer, but a detailed model 
of these dynamics is beyond the scope of this chapter. The pressure decreases to 
zero at the boundary and it remains zero for larger radii. This means that infall or 
expansion around the cluster is radial, moving mainly with the Hubble flow, and 
has negligible random velocity. 



The boundary of the generalised isothermal sphere is physically defined by 
P(r = r b ) = 0: 



The pressure parameter at the centre, a, is directly related to the scale of the 
generalised isothermal sphere. For vanishing a, Eq. (5.28) agrees with the fact 
that in this limit there is no isothermal sphere. Larger values of the constant a 
correspond to a sequence of static solutions for larger spheres which have grown 
through the merger of clusters at later stages. We find that C and D can be used to 
normalise the scale of the generalised isothermal sphere. We will also discover that 
C and D are not independent, since the matching to the Schwarzschild exterior 
provides a relation between them. 

The mass of the generalised isothermal sphere can be determined by integrating 
its density up to the boundary: 



Given the functional form of the isothermal boundary, r b , in terms of C, D and 
a we can rewrite the generalised isothermal metric in the following form which 
closely resembles the original isothermal metric, but with a radially dependent 




l+q 



(5.28) 




(5.29) 
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deviation term. 



e = 



4g 

Ar 1 + a 



1 + 



1 - 
4a 



a 



2(1-") 
1 + a 



8np 



1 + a) 2 ' 
4a 1 



4a + (1 + a) 2 r 2 ' 



P = 



ap 



1 _ 21 



2(1-") 
1 + a 



1 - a 2 



2(l-a) ' 
1+a 



(5.30) 

(5.31) 
(5.32) 

(5.33) 



where we have set C 2 = A. 

Hence we have found the interior metric for a generalised isothermal sphere of 
finite radius: 



ds = —Ar l + a 



for r < n,. 



I -a' 



2(1- a) -i 2 

y \ 1 + CK 

n 



dt 2 + l 



4a 



1 + a) 2 



rfr 2 + r 2 (^ 2 + sin 2 ^0 2 ) 
(5.34) 



§5.3 The Isothermal Approximation 

Our generalization of the isothermal universe involves two approximations 
which we now want to consider in turn. First we review the approximation of 
the generalised metric to the ideal isothermal counterpart. For the generalised 
isothermal metric to be a good approximation to Eq. (5.8), Eq. (5.30) requires 
the condition 

2(1-") 



j = a 2 « 1, (5.35) 

where < a < |. This condition on a incorporates the possibilities of relativistic 
pressure (a = 1/3) and of a perfect monatomic gas (a = 2/3). 

A contour plot of j (r/r b , a) is given in Fig. 5.1. It shows that j < 1 for all a and 
(r/r 6 )-ratios. The generalised metric therefore resembles the perfect isothermal 
metric very closely. Deviations from the perfect isothermal metric express the 
feedback of the finite boundary of the isothermal sphere on the metric. 
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Figure 5.1: Contour plot of j. Here x = (r/rj,). For j <C 1, which characterizes 
almost all combinations of a and x, the deviation of the generalised metric from 
the ideal isothermal form is neglegible. 



The generalised equation of state can be rewritten in the form P = ap [1 — A] 
where we define the deviation term 

2(l-q) 

(1 - a 2 ) (-) 1+a 
A = W 2(1 _ a) • (5-36) 



1 - a 2 - 



2(1-") 
1 + a 



Next we define the transformation x = (r/r^) to obtain the normalised radial 
gradient of the deviation: 

OA 2(l-a) 2 x 



8X (*2- n ^_\ 2 ' 

x x + a — a z x l + a 



(5.37) 



For a good approximation to the isothermal equation of state p = ap, we suggest 
two necessary limits: (a) A — > which scales the pressure- density dependence and 
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(b) \dA/dx\ <C 1, which expresses the fact that in an ideal isothermal sphere, the 
form of the equation of state is independent of position. 



. 6 
0.5 
. 4 

a 

0.3 
0.2 
0.1 
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X c 

Figure 5.2: Contour plot of the deviation term A. The plot gives the char- 
acteristic core region where the deviation from the perfect isothermal equation 
of state is less than a particular contour value. Note that although the ratio of 
the core radius to the boundary radius x c = (r c /rb) decreases with increasing a 
(and hence with step number), both r c and r b are monotonically increasing with 
each step. The decreasing x c merely reflects their relative rates of increase. 

Fig. 5.2 shows a contour plot of the deviation term A. We see that for a partic- 
ular value of the constant a the approximation gets worse with increasing distance 
from the center, as the radius- dependent perturbation becomes more significant. 
The plot gives a characteristic core region where the deviation from the perfect 
isothermal equation of state is less than a particular contour value. Note that 
although the ratio of the core radius to the boundary radius x c = (r c / r b ) decreases 
with increasing a, at later evolutionary steps, both r c and r b are monotonically 
increasing with each step. The decreasing x c merely reflects their relative rates of 
increase. 

Fig. 5.3 shows diagrams of the radial derivative of the deviation term dA/dx. 
This derivative becomes very small for both small a and small r/r b . For large a 
and r/rt, there is a significant variation of the equation of state with position. The 
derivative dA/dx shows the complex functional form expected for a finite sphere. 
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Figure 5.3: Contour plot of the derivative of the deviation term (dA/dx). 
The plots gives the core region (small a, r <C r^) where the derivative is less 
than a particular value and the ratio of pressure to density is therefore nearly 
independent of position. 



The contour plots of Figs. 5.1 - 5.3 indicate that there exists a core region 
within our generalised sphere for which the metric and the equation of state are 
close to the ideal isothermal limit. There are several alternative ways of denning 
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this core region. All of these arbitrary definitions, however, have in common that 
they describe the same qualitative picture. In that picture our metric describes a 
sphere which tends to the perfect isothermal limit at the centre, has a core region 
where it is close to the isothermal state and then a transition region approaching 
the boundary where the isothermal approximation breaks down. 

For small a, a good approximation to the ideal isothermal equation of state 
requires 



We may ask for the range of r < r c (relative to r&) that satisfies a fixed constraint, 
for example, k c = const = (1/3). This gives a core sphere, r c , within the isothermal 
sphere for which the approximation does not deviate by more than the set limit. 
This critical core radius is given by 



As r — > 0, Eq. (5.38) implies that A — > 0, and the approximation to the isothermal 
equation of state goes to its ideal limit. For larger r the approximation gets 
continuously worse at a rate described by Eq. (5.38). Both the outer radius and 
the core radius r c grow with each step as a and A change. Moreover, from equation 
(5.39), we see that r b expands faster than r c and the 'wobble' region between the 
isothermal core and the boundary radius grows. 



5.4 Thermodynamic Consequences implied by 
the Generalised Equation of State 



Although an isothermal equation of state is well-known to be a reasonable 
zeroth-order approximation for rich clusters of galaxies (like the Coma cluster) 
it is also known not to be realistic, especially in the outer regions, for clusters 
with finite boundaries. From a physical point of view, this is because the longer 
dynamical relaxation timescales at the lower densities in the outer regions do not 
allow sufficient energy exchange for thermalization. From a mathematical point of 
view, the departure from isothermality is forced by the outer boundary at a finite 
radius. 

The generalised equation of state can be written as 




(5.38) 



l + q 




(5.39) 




(5.40) 



where 



a e g(a,r) 



= a ■ 




2(1-") 



(5.41) 
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We determine the limits: (i) r = 0, a c g — > a, (ii) (r/r&) <C 1, a c g- ~ a, and (iii) 
(r/r b ) — > 1, a c fi — > 0. Note that we have implicitly introduced a new definition of 
a: 



a ^ ^4 (5.42) 
p(r = 0) y ' 



i.e. a characterizes the constant of proportionality between pressure and density 
at the centre of the sphere. 

Eq. (5.40) can be regarded as an isothermal equation of state with a tem- 
perature that depends on r according to the fractional term in (5.41). This type 
of temperature behavior is quite common in problems involving boundary condi- 
tions. In our case, the pressure is zero at the boundary since the temperature is 
zero there for a good physical reason, namely the galaxies there have not relaxed 
and randomised their peculiar velocities significantly since the density of galax- 
ies is low and there are fewer gravitational scatterings. The peculiar motions of 
galaxies near the boundary are mainly streaming motions, probably mostly infall 
into the approximately isothermal cluster. 

Note that for any constant r (or given temperature), Eq. (5.40) shows that 
dPj dp is positive so this equation of state has a positive isothermal compressibility, 
as needed for mechanical stability. Because Eq. (5.40) is close to an isothermal 
equation of state, and equivalent to it at any fixed radius, it should have reasonable 
thermodynamic properties as long as the number density of galaxies is large enough 
that one can define the thermodynamic properties in a volume smaller than that 
of the range of radii one is considering. This is similar (though not exactly the 
same) as defining the local temperature in a gas which has a temperature gradient. 

All that the equation of state is saying is that the temperature is a function of 
position. This is similar to the case of polytropes in stars, where it is sometimes 
generalised to include inhomogeneous polytropes [36] . In real stars the temperature 
dependence of the equation of state on position occurs because of the nuclear 
reactions and energy transport processes. In our largest cluster it occurs because 
the relaxation timescales are dependent on the local density. Calculating this 
is a separate problem. For the equation of state to be physically acceptable, it 
must satisfy the first and second laws of thermodynamics locally. An isothermal 
equation of state does this, although our case is complicated by the long range 
nature of gravity in the cluster. We can regard the isothermal case with variable 
temperature as an approximation, just as it is in stellar structure. 
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§5.5 A Solution of an Expanding, Centrally 
Isothermal Universe 

In this section, we derive a hybrid metric which encompasses the behavior 
of both an expanding universe and a static nearly isothermal core. Physically, 
we consider the evolution of such a hybrid metric to occur through three stages. 
First, the universe follows the usual globally homogeneous Einstein-de Sitter ex- 
pansion, which tends asymptotically to stationarity, i.e. the expansion timescale, 
a/ a, of the universe tends to infinity at late times. Second, the condensation 
stage (which is attributed to galaxy clustering) gradually produces a significant 
cosmological pressure which alters the initial global homogeneity and supports a 
radially inhomogeneous global distribution of matter. Consequently, the equation 
of state changes with time, resulting in a dynamical feedback on the metric. Third, 
the postulated asymptotic final state is the generalised isothermal sphere with a 
boundary characterized by the equation of state in Eq. (5.27). Our solution only 
applies to stages one and three. It only describes the relaxed states of the universe 
with a static approximately isothermal core. The temporal evolution will be a 
sequence of such relaxed, static stages. We do not have a continuous description 
in time, but only discrete steps between quasi-static periods with constant a. Our 
solution breaks down in between. We note that the sequence of hybrid metrics 
(with different values of the constants a and A) is constrained by various boundary 
conditions in both space and time, namely, (a) at early times (for t — 0) the metric 
is FRW for all r, and (b) at late times (for t — > oo) the universe would have a 
bounded approximately isothermal metric with a horizon expanding proportional 
to t, and be FRW on large scales (r — > oo) which increase at tW. 

A possible counter argument to this proposal would be the following: the 
horizon of the Einstein-de Sitter can expand rapidly enough that the matter inside 
is not consolidated into a single bound system, no matter how long one waits. 
The only loophole is if the primordial perturbation spectrum tilts towards long 
wavelengths. However it would seem at first sight that the scenario we propose is 
not valid. 

However, the key point against the skeptic's argument is that we are consid- 
ering the very non-linear case of structure formation after very long times when 
the clustering parameter b — > 1 and nearly all the matter is in a few dominant 
massive clusters. The metric still expands at the Einstein-de Sitter rate which is 
proportional to t( 2 / 3 ) since this is determined to a first approximation by an aver- 
age over all horizons, so the approximate expansion timescale, t expansion , is given 
by the usual 

2a 1 



116 



5. CHANGING GLOBAL SYMMETRY 



We wish to compare this expansion with the timescale, t merge , for the clusters to 
merge. We are interested in the non-linear merging of smaller galaxy clusters 
with a larger cluster in the expanding Einstein-de Sitter Universe. In the simplest 
approximation, as an example to estimate the basic merging timescale, t merge , for 
collapse into a single cluster, consider a massive cluster surrounded by a roughly 
spherical distribution of less massive clusters, where the whole configuration has 
an average large initial overdensity, ft, compared to the initial critical average 
density p ci in that volume. For simplicity, suppose the clusters have negligible 
initial comoving velocity. We know that this region which has gradually become 
overdense will subsequently expand more slowly than the universe, turn around 
and collapse leading to a merger of the cluster. 

To estimate the collapse (merger) timescale, we use the illustrative example in 
[168] (See Eqs. (37.10)-(37.19) in the reference). For p t > p ci (and thus E « 1 
where E = a max ), the time scale to maximum expansion is 



(5 - 44) 

where pi is the initial average density in a volume containing the clusters, and 

2 71-2 ^ 37T 

tmax = T~8^Gp~ = 32Gp~ ' ( ' 



Thus 



t2max ~* x (6ttGVw) = » 5.65^ (5.46) 



Expansion 32G P* ' ' Hi p 



and we obtain 



tmax ~ 2.36 p ^ ^expansion ■ (5-47) 

If the clusters are a factor of ~ 2 closer than the average separation, pi pa 8p cr u 
and t max ~ 0-8t exp ansion- The timescale for the clusters to merge is less than about 
2t max . If the clusters have initial velocities (which they will probably acquire while 
they form) that are mainly toward the most massive cluster, then this merging 
timescale is shorter. Therefore, any clusters that have less than about half of 
the average separation of clusters (i.e. where pi > 8p cr it) are likely to merge 
within one expansion timescale. Clusters with larger initial separations will merge 
after longer times. On these longer timescales, the horizon scale will double on 
a timescale equal to t expansion , which becomes infinite. This may or may not be 
the timescale to introduce new clusters into a particular horizon, depending on 
where they are located in the adjacent horizons. In the asymptotic limit, t — > oo, 
which we are interested for this chapter, an arbitrary large horizon can contain an 
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arbitrary large dominant cluster, separated from such clusters in other horizons 
by arbitrarily large distances. 

What is the effect of the matter that enters the horizons? At any given time 
within a horizon centred on a dominant supercluster, the galaxies not in the super- 
cluster will be accreted by it on a timescale less than the expansion timescale. If 
the accretion is nonlinear, e.g. if the supercluster dominates the local gravitational 
field, this accretion timescale is shorter than the average expansion timescale by 
a factor of about \/2. At any finite time, the infinite E-deS universe will contain 
an infinite number of these expanding horizons. The positions of their central su- 
perclusters throughout the whole universe may have a Poisson distribution since 
they have not had time to influence each other. However, they could, have long 
wavelength primordial perturbations which would speed up their clustering. Some 
of the horizons at any given time will be dominated by one supercluster, others 
may contain more, or may contain merging superclusters. Outside any of these 
horizons, the other clusters may be averaged over an increasingly large scale to 
give the effective average density of the E-deS universe. 

In the limit of infinite future time (which is what we consider), for a Pois- 
son distribution of superclusters initially outside each other's horizons, they come 
within each other's horizons and can merge into our isothermal sphere solution 
which grows to ultimately incorporate all the matter of the universe in the limit 
of infinite time. If long wavelength initial perturbations are important, this could 
happen sooner. In either case, the average separation of superclusters increases 
more slowly than the horizon expands, which is what our model requires. 

§5.5.1 The derivation of the hybrid metric 

Bounded isothermal universes with non-zero pressure (but with zero-pressure 
and non-zero density at the boundary at r = r h ) cannot be matched directly to 
the Einstein-de Sitter model across a specific r = constant hypersurface. How- 
ever, one can achieve a hybrid solution by first matching the isothermal metric to 
a Schwarzschild exterior. This is physically justified by Birkhoff 's theorem. A for- 
mal mathematical treatment is given by Fayos et al [68] for the general matching 
of two spherically symmetric metrics for various pressure and density conditions. 
One can perceive the Schwarzschild metric as a kind of membrane (equivalent to 
the requirement of a closed trapped surface discussed by Fayos et al [68] for the 
spherically symmetric solution, such as the Vaidya metric, to match onto the FRW 
metric) which facilitates the matching between the bounded static isothermal sphere 
and the expanding Einstein-de Sitter universe. Such a procedure would be analo- 
gous to the astrophysical case for an isothermal star whose external Schwarzschild 
metric can be embedded in its small region of the universe. 

The continuity condition of the metric coefficients, requires the following rela- 
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tions to be satisfied: 

( 2a 2 \ 2 Offl 

Cr i+° _ Dr i+° j = 1 - — , (5.48) 

^) = l + 7 4 \ =fl-— V 1 . (5.49) 

(l + «) 2 V 1 J 

where r fe is denned by Eq. (5.28). These provide the matching relation between C 
and D : 

D(a)=E(a)C«, (5.50) 

where 



E(a) 



l-a 

4a 

1 + 



2a l-a 

a 2 ri ~ 2 



[l-a 2 ] Q . (5.51) 



l + a) 2 

Substituting the matching condition into Eq. (5.28) we find 

r b {a) = F(a)A-^ (5.52) 

where A = C 2 and 

l + ct 

( a 2 \ 2 ( 1 -") 
F(a) S (^) (5.53) 

with £'(a) defined in Eq. (5.51). The matching eliminates the constant D and 
therefore shows that the boundary radius r b is only dependent on the parameters a 
and A. Hence, we also deduce immediately that the generalised isothermal metric 
(5.34) is parameterised by the two constants a and A. 

This procedure also determines a Schwarzschild mass from the matching con- 
dition: 

m(a) = ( Aa + \* +a)2 ) n = M lso (5.54) 

Hence, the Schwarzschild mass m(a) is equal to the mass of the bounded gener- 
alised isothermal sphere, M iso , as found independently by direct integration of the 
density. The equality of the results indicates consistency in our approach. The 
corresponding Schwarzschild radius is 

r s = 2m=(- — — — ) r b , (5.55) 



4a + (1 + af 

which indicates that r s < r b , i.e. the Schwarzschild radius is always inside the 
bounded isothermal sphere (as is the case for a star). 

To complete the picture, we adapted the method in McVittie [145] to enable 
us to embed the finite generalised isothermal sphere with given properties in the 
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Einstein-de Sitter model. A sequence of embeddings with different properties is 
later used to describe the time evolution of the universe. McVittie's solution 
applies to a mass particle in an expanding universe and has been applied by No- 
erdlinger and Petrosian in [150] to investigate the motion of galaxies in galaxy 
clusters under the influence of cosmic expansion. In [145], he found a way of em- 
bedding a Schwarzschild geometry in a FRW spacetime. Having matched the gen- 
eralised isothermal metric to a Schwarzschild exterior, we use McVittie's approach 
to embed this exterior spacetime in an expanding Einstein-de Sitter cosmological 
background. The most general form of the McVittie metric was summarised in a 
later paper by McVittie in [146] which incorporates a larger class of such solutions. 
Hence, we find 

(i _ m(a) \ 2 / / n. \ 4 

rft 2 + ^ 1 + ^Jj a 2 (t)[df 2 + f 2 (d6 2 + sin 2 9d<p 2 )] (5.56) 

for the exterior (r > r^) matched solution of a generalised isothermal sphere of 
mass m(a) in an expanding FRW universe. Note that the solution is stated in 
isotropic coordinates. Here a(t) is the scale factor of the standard cosmological 
solution and m(a) is defined by Eq. (5.54). 

If one makes the substitution f\ = a ■ f, the line element in Eq. (5.56), at 
any particular instant t±, reduces to the Schwarzschild line element for small f 
(isotropic coordinate). McVittie also showed that the transformation of the time- 
dependent metric in Eq. (5.56) into the static Schwarzschild field is independent 
of the moment in time at which the transformation is made in [161]. This proves 
that the embedding is valid for all times t. 

To relate the interior generalised isothermal metric and the exterior metric, we 
give the transformation between the spherical polar radial coordinate r used in 
Eq. (5.34) and the isotropic coordinate f\ in Eq. (5.56): 

\ 2 

m 



^(1 + — ) , (5.57) 



1 



y/r 2 — 2mr + (r — m) 



(5.58) 



2 

It can then be shown that the interior and exterior metrics match at r = r b , for 
a particular instant t\ and that the matching is valid for all t in [151]. For a 
given value of a = constant, the derivatives of the metric coefficients also match 
according to the junction conditions described by Bonnor and Vickers [16]. In 
the next section, we will discuss the matching of the metrics in a more rigourous 
manner. 

At the beginning of this section we discussed the limits in space and time that 
the sequence of hybrid solutions has to satisfy. We now compare these constraints 
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to the corresponding limits of Eqs. (5.34) and (5.56). For vanishing a, the mass 
and the boundary radius of the isothermal sphere tend to zero and the metric is 
FRW on all scales. For small a (i.e. early evolutionary steps) the applicable hybrid 
metric has the generalised isothermal form for small r and is FRW for large r. For 
relatively large a (i.e. late evolutionary steps) the metric tends to the isothermal 
form on much larger scales. We therefore see that all limits of our hybrid solution 
in both time (discrete steps between static periods) and space agree with the limits 
required by physical arguments. 



§5.5.2 The validity of the matching of the metric 

We may ask more generally whether using the generalised isothermal sphere pa- 
rameters in the Schwarzschild interior metric could lead to a problem in matching 
the FRW exterior. Mars [143] has shown that the only static, spatially compact, 
vacuum region in a Friedman- Lemaitre spacetime must be a 2-sphere comoving 
with the cosmological flow and with the Schwarzschild- exterior geometry. In our 
model, at the boundary between the Schwarzschild metric and the bounded isother- 
mal sphere, there are no problems that arise from the energy momentum tensor. 
This is because at that boundary between the bounded isothermal sphere and the 
Schwarzschild exterior, the pressure P is zero. Therefore, there is no problem em- 
bedding the matched system of a bounded isothermal sphere and the Schwarzschild 
metric into the FRW metric, if the solution of the combined system is expanding. 

Since the energy momentum tensor is not an issue in the matching, we can 
work out the intrinsic and extrinsic curvatures of the metric by matching on a 
hypersurface E, at a constant comoving radius, r. The smooth matching of the 2 
spacetimes across a hypersurface E is guaranteed if the Darmois junction condi- 
tions are satisfied, i.e. if the first fundamental (intrinsic metrics) forms and second 
fundamental (extrinsic metrics) forms are calculated in terms of the coordinates on 
the hypersurface, and they are identical on both sides of the hypersurface. If both 
the first and second fundamental forms are continuous on E, then the matching 
is not piecewise i.e. it will match for all times rather than just at single instants. 
The Darmois junction conditions allow us to use different coordinate systems on 
both sides of the hypersurface. We adopt the procedure in [55] and have calculated 
explicitly the interior and the exterior metrics and they match on a 3-dimensional 
timelike tube and exhibit no curvature discontinuity for all times. In the pro- 
cess, we applied the procedure in [55] to check the conditions for matching the 
Schwarzschild- Tolman system and the FRW metric on a hypersurface E, where 
r = To for the radial constant comoving coordinate in the FRW metric. 

We recall the Schwarzschild- Tolman metric (bounded isothermal sphere interior 
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with a Schwarzschild exterior), which is written in coordinates (T, p, 9, </)), 

d# = -(l- dT 2 2 _MM^ 1 d( ? + p 2 M 2 + p 2 sin 2 g d(j) 2 ^ (5>5g) 

where 

M (") = 4a + ( i + a) ^ • < 5 - 60 ' 

and Tb is a constant. The general FRW metric is given in Eq. (1.1). 

The first fundamental form is the metric which the hypersurface £ inherits 
from its spacetime in which it is imbedded: 

dx^ dx v 

^ = 9 ^frJ> (5 ' 61) 

where u l = (w 1 , w 2 , w 3 ) = (u,v,w) is the coordinate system defined on the hy- 
persurface. Note that i,j denotes only the coordinates of the hypersurface of 
3-dimensions while p, v denotes the coordinates of the 4-dimensional spacetime. 
The second fundamental form is defined by 

% = (r>,-^) — — , (5.62) 

where is the unit normal to the hypersurface, S. Finally if £ is given by the 
function f[x^{u % )] = 0, then n M can be calculated by: 

d f 

n„ = ; » J (5.63) 

V^djdfsf 

We will denote the FRW and Schwarzschild- Tolman (ST) frames by the subscripts 
or brackets with F and S in our quantities. 

We consider a spherical hypersurface £ given by the function Jf^f) = r ~ r o = 
0, where r is constant. 

For the FRW frame 

x F = t = u, x r F = r = r , x F = 6 = v , x F = <f) = w, (5.64) 
and similarly for the ST frame, 

x t s = T = T(u), x p s = p = p(u), x e s = 9 = v, x% = <f> = w. (5.65) 
The first matching condition ^(F) = 7^(5*) gives 

-^)(fy-o-^ro 2 -' <«•«> 

a(t) 2 r 2 = p 2 . (5.67) 
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Next we calculate the second fundamental forms. The (outward pointing) unit 
normal in the FRW frame follows from (5.63), and since /f(^ ) = r — r = 0, we 
get n^F)=5\n r (F). 

For the unit normal in the Schwarzschild case, we cannot derive the result 
directly from Eq. (5.63), but we use the result that n^(S) must satisfy the following 
conditions: 

n"(F)n„(F) = n^S)n^(S) = 1 , (5.68) 
n,(S)^ = 0, (5.69) 

where the second condition comes from the fact that partial differentiation of 
fs[x % s{u 1 )} = with respect to u\ With these conditions we get the following 
equations in component form: 

n e (S) = n^S) = , (5.70) 
2M(a)Y 1 2/CA / 2M(a)\ . 



1 ~—^ L ) r4(S) - [1 - ) ni(S) = -1 . (5.72) 
Comparing the last equation with Eq. (5.66), we derive n^(S) to be 

"*< s > = {-% "f • «• °) • <" 3 > 

We can compute f2 from Eq. (5.62) to be 

% = ---^=d r gij , (5.74) 

and this equation is valid for any coordinate hypersurface r = r =constant, in an 
orthogonal coordinate system and parametrised by x % = u l . 

From Eq. (5.69), we obtain the following equation by differentiating with 
respect to u\ 

^ } ra = -^» (5 - 75) 

and we get 

n (q\ dx u dx x d 2 x^ 

ihj lb = T,,^n u — — — — + n u - — . . 5.76 

We know that Qij(S) = Qij(F) = for % ^ j, and for the remaining diago- 
nal components, we have flij(S) = Qij(F) on E, and we get the following three 
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differential equations: 

n „ ^ fdT\ 2 ^ fdpV rp dTdp d 2 T d 2 

(5.77) 

n vv (S) = n vv (F) = {\a(t)\r = T p ee (S)n p , (5.78) 

tt ww (S) = n ww (F) = (\a(t)\r sin 2 6 = T P H (S)n p , (5.79) 

where in the Schwarzschild frame, the Christoffel connections are 

r - = ^ = ^SR) • < 5 ' 80 > 

P 

T^ = -(p-2M) , (5.82) 

r^ = Csin 2 ^ (5.83) 



and ( = \/l — &?"o- Note that this will also prove that the matching is continuous 
for the cases where k ^ 0, i.e. open and closed cases as well. In retrospect, if 
we follow the procedure of Hogan[89], we will also obtain the generalization of 
the Schwarzschild- Tolman metric with the open and closed FRW universes, but 
their forms will be complicated. By the same physical reasons mentioned in the 
discussion of the infinite isothermal sphere solution, the open and closed FRW 
universes would deter the phase transition from happening. 
Using Eqs. (5.78), (5.67) and (5.66), we get 

dT (p 



du p- 2M(a) 



(5.84) 



du) e-K-r^v ^ 



dp\ 2 (p-2M(a) 



Differentiating Eqs. (5.84) and (5.85) gives 

d 2 T 2M(a) dp 



(5.87) 



du 2 ~ s (p - 2M(a)) 2 du 
d 2 p M(a 



du 2 p 

Because M(a) =constant for values of a, we can substitute Eqs. (5.84) and (5.85) 
into (5.77) to show that Eqs. (5.77), (5.78) and (5.79) are consistent. From these 
solutions, one would realise that both the first and second fundamental forms are 
continuous on S, hence the matching is not piecewise, i.e. it will match for all 
times rather than just at single instants. 
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§5.6 Evolution of the Isothermal Universe - A 
Sequence of Hybrid Solutions 

Thus far, we have derived the hybrid metric to describe a generalised isothermal 
sphere embedded via a "Schwarzschild membrane" in the expanding Einstein-de 
Sitter universe at any particular cosmic time. Now we consider briefly how such a 
hybrid model might form and evolve. 

Greater than average densities in any volume of an Einstein-de Sitter universe 
produce local gravitational clustering timescales shorter than the current expansion 
timescale. This facilitates the growth and merging of a small number of large 
clusters within an apparent or particle horizon. Collective relaxation [168] could 
thermalize the velocities of the galaxies on a timescale ~ (Gp)~^ where p is the 
average density in a volume containing the merging clusters. Eventually, this 
would lead to an approximately isothermal virialized cluster. This process would 
continue indefinitely as other great clusters appear within the horizon and merge. 

During periods when one massive isothermal cluster dominates, our idealised 
hybrid metric would apply. More generally, a much more complicated metric would 
be needed to describe the periods when two or more great clusters dominate a 
horizon. However, we may envisage our hybrid metric applying at intervals, for 
particular values of a and A when one cluster dominates. "Snapshots" of the long 
term evolution would then be characterized by a discrete sequence of increasing 
values of a as virialization progresses, and consequently by decreasing values of A, 
as seen in Fig. 5.4. 

How long will it take an isothermal universe to form? We can estimate this 
timescale and compare it to timescales for several possible competing effects. The 
timescale for the formation of the isothermal sphere can be estimated from the evo- 
lution of iV-body galaxy clustering. For clustering in an adiabatically expanding 
universe this gives [169] 



where a is the time- dependent scale factor of the universe and b is defined as 
the average ratio of the gravitational correlation energy, W, to twice the kinetic 
energy, K, of peculiar velocities in an ensemble of galaxies. Inverting this gives 
the time-dependence of b(t) within the limits b(0) = (no interaction, perfect 
gas) and ft(oo) = 1 (complete virialization, isothermal). The transition to the 
isothermal state therefore corresponds to a b close to unity, where the universe is 
dominated by voids and a small number of very large clusters. We define r as the 
time after which the void probability of the universe becomes greater than 0.5. 
This describes the timescale required for an approximately isothermal cluster to 
grow to cosmologically significant size. From the galaxy distribution function in 



a(t) = do 



(1-6) 



7 ' 

8 



(5.88) 
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[169], we get the condition 

e -N(i-b) > (5>g9) 

where N is the total number of average mass galaxies in the observable universe. 
If we assume N pa 10 10 , this allows us to define the condition for transition to the 
isothermal state as b c = 1 — 1CT 10 . Hence, we find the scale factor a c (r) = 10 8 a , 
where a is the present Hubble radius. For an Einstein-de Sitter expansion with 
a(t) = ts this corresponds to 

t exp ^ 10 12 t = 3 x 10 29 seconds, (5.90) 

where we used an estimate for the present Hubble time to — 10 17 seconds. Hence, 
the universe would take approximately 10 29 seconds to form the isothermal sphere. 

The solution permits a sequence of generalised isothermal universes whose radii 
grow with each step according to 

fa 2 C\^^ 1+a 

r b = =F(a)A^, (5.91) 



D 



where 



1 + a 
v 2 \ 2(l-a) 



S ( ^ ] (5.92) 



and E(a) is defined by Eq. (5.51). The values of r& increase as a and A change 
their values with time (discrete steps: for each stage a and A are constant; see Fig. 
5.4). We explore next how this constrains the time- dependence through discrete 
steps of the sequence of a and A = C 2 values. 

We can normalize our metric by finding approximate limits on the sequence of 
A. These limits can be found by the following approximate argument: = at 
early times requires lim t ^ A t > 1 (see Fig. 5.4). We next consider the limit for 
late times. From the previous discussion we can set r b ~ a c = 10 8 a = O(10 36 )cm 
as an approximate scale of the isothermal universe at late times (t = t c ). 

This allows us to use Eq. (5.91) and find A tc at a = (2/3), 

A tc ~ 0(1(T 58 ) cm"!, (5.93) 

which gives an approximate limit for the time evolution of the sequence of the 
constants A t . The proposed time evolution of A t therefore goes as follows: A 
starts off at t = with A > 1. At t c (corresponding to a(t c ) = 2/3), A tc ~ 
O(10 -58 ) cm" 5. For t > t c , A t decreases further, so that lim^oo A t — > and 
lim^oo ^(t) — > oo, meaning that the isothermal sphere grows indefinitely. Fig. 
5.4 should therefore be interpreted as a map of the future scale of the isothermal 
universe for changes in a and A. Note, however, that our present state of the 
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universe is still very close to a — 0, r& = 0. As the universe evolves in the far 
future, changes of a and A correspond to an increasing scale of the isothermal 
part of the universe. This evolution as a sequence of quasi-static steps can be 
characterized by discrete points on the surface of Fig. 5.4. 



log(r&) 




The Present 

Figure 5.4: Surface plot of - the model's future laid out before us: It is 
interpreted as a map of the future scale of the generalised isothermal universe 
for particular values of a and A. The radius of the generalised isothermal sphere, 

— 4c* 

r;,, is in units of centimeters, A is measured in [cm 1 + a ] and a is dimensionless. 
rb and A are plotted on a logarithmic scale. The hypothetical present size of the 
isothermal sphere is taken to be ~ 10 25 cm (the size of the Coma cluster). 

There are competing effects, which could inhibit the formation of the gener- 
alised isothermal state, namely, (a) the conversion of all the matter into gravita- 
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tional radiation, (b) Hawking radiation and (c) proton decay and (d) a positive 
cosmological constant. 

Next we estimate the amount of gravitational radiation produced by the inter- 
action between a pair of galaxies. Using the standard quadrupole formula given by 
Landau and Lifschitz [114] and making reasonable estimates for masses and mean 
distances of galaxy binaries we find that it would require about 10 51 seconds for 
the universe to convert all its mass into gravitational radiation. This is 10 22 t exp . 
The mass which could be evaporated by Hawking radiation of small black holes 
in the time needed to reach the isothermal phase transition is calculated as 10 28 
grams which is equal to the mass of the Earth and can therefore safely be ne- 
glected. The isothermal state is therefore reached long before Hawking radiation 
could have influenced the mass content of our universe. Some grand unified the- 
ories suggest the instability of the proton. Experimentally, however, it was found 
that the lower limit for the half-life of the proton must be greater than 10 33 seconds 
(if not infinity). That leaves us with only one more possibility, a positive non-zero 
cosmological constant. If we assume = 0.7 in standard big bang cosmology, 
the cosmological constant would become dominant after the radius of the universe 
increases to about three times its current value. This is significantly faster than 
the timescale of the formation of the isothermal sphere. It suggests a possible 
role for A, as a buffer in preventing our universe from undergoing a possible phase 
transition. 



§5.7 Discussion 

We have presented a relativistic approach to describe the evolution of the uni- 
verse from an Einstein-de Sitter expansion with no centre to an isothermal sphere 
with a centre. We have shown how this changing symmetry of the universe can 
be incorporated into general relativity. This demonstrates that general relativity 
permits a changing symmetry on a global scale, a nontrivial result. 

Our solution required a Schwarzschild metric to act as an intermediate back- 
ground to match the static generalised isothermal metric to the FRW metric of the 
expanding universe. The necessity for this Schwarzschild intermediary can be un- 
derstood both intuitively and on a formal mathematical basis. The Schwarzschild 
metric acts as a kind of membrane linking a static space-time region to an expand- 
ing space-time. In [68, 143, 144], the general matching of two spherically symmet- 
ric spacetimes was investigated in detail. Their results provide a formal basis for 
understanding the role of the "Schwarzschild membrane". Applying the general 
matching conditions derived in their paper to the expanding, centrally isothermal 
universe proves that a direct matching without the "Schwarzschild membrane" 
between the generalised isothermal and the FRW spacetimes is not possible. An 
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extension of this idea of a "Schwarzschild membrane" to general spherically sym- 
metric spacetimes other than the isothermal metric might be possible and have 
interesting implications. 

In our hybrid model the field equations are solved for constant pressure/density 
ratio, a and, as stated, the solution therefore does not incorporate a time- dependent 
a. We therefore introduced a "sequence of quasi-static stages". We envisage the 
universe evolving through a sequence of models with increasing a as its central 
cluster grows in size and the isothermal approximation becomes better for larger 
regions. The main physical argument for this is that the local clustering timescale 
in a nonlinear positive density inhomogeneity in an Einstein-de Sitter model is 
shorter than the expansion timescale since the average local density is greater 
than the average global density. Our metric only describes those phases of the 
evolving universe when it is dominated by one relaxed, static isothermal cluster. 
As a second large cluster enters the horizon our description breaks down and only 
applies again after the system has relaxed to a larger static cluster with new values 
of a and the boundary radius of the cluster. This leads to a sequence of discrete 
changes which characterize the evolving model. Our solution therefore describes 
only discrete steps between static intervals and not continuous changes in time. 
Finding a continuous description remains a problem for the future. 



Part IV 
Epilogue 



Chapter 6 
Discussion 



"After a lifetime of crabwise thinking, I have gradually become aware 
of the towering intellectual structure of the world, One article of faith I 
have about it is that, whatever the end may be for each of us, it cannot 
be a bad one. " 

- Fred Hoyle 

"/ do not know what I may appear to the world. But to myself I seem to 
have been like a boy playing on the seashore, and enjoying myself in now 
and then finding a smoother peeble, or a prettier shell than ordinary, 
while the great ocean of truth lay all undiscovered before me. " 

- Sir Issac Newton 

From the Beginning to the End 

In the epilogue, we summarize the main results of this thesis and discuss pos- 
sible future directions. We began this thesis with a modest aim to study possible 
and viable alternatives to the standard inflationary big-bang cosmology, in the 
hope that we can understand these alternative cosmologies better and search for 
phenomenological effects if possible. As such, we have explored three interesting 
alternative themes, namely, spin-torsion theories, extra dimensions and changing 
global symmetry. 

Beginning from chapter 2 for spin-torsion theories, we found new massive, 
cosmological solutions for the Dirac field coupled self-consistently to gravity with 
a cosmological constant in a K = (flat) universe. In addition to the generalized 
particle and anti-particle solutions found in [35] with a cosmological constant, we 
discovered another anti-particle solution which has interesting implications. All 
these solutions have the presence of a particle horizon but they do not possess 
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an inflationary phase. Although these spin-torsion solutions are interesting, they 
are shown not to be viable alternatives to the standard scalar field inflation. It 
has been shown in [2, 152] that the vacuum polarization terms in the framework 
of Einstein-Dirac equations would provide an alternative to standard inflationary 
models. Hence the next step is to reconsider the model with vacuum polarisation 
terms. 

Next, we come to the theme of theories with large extra dimensions. We con- 
centrated on a sector of the braneworld cosmology in chapters 3 and 4, particularly 
on the study of linear perturbations using the (l+3)-covariant approach and the 
possible effects on the CMB from the extra dimensions. In chapter 3, we set up 
the (l+3)-covariant formalism for the study of cosmological perturbations, and 
provided a complete set of frame-independent equations for the total matter vari- 
ables, and a partial set of equations for the non-local variables which arise from 
the projection of the Weyl tensor in the bulk. We also derived the covariant form 
of the line of sight solution for the CMB temperature anisotropies in braneworld 
cosmologies. However, our setup still lacks the propagation equation for the non- 
local anisotropic stress, which appears to be an important ingredient in the study 
of CMB anisotropies. 

In the following chapter, we studied the scalar and tensor perturbations and 
their implications on the CMB anisotropies. For the scalar perturbations, we 
provide the equations for the total matter variables with equations for the inde- 
pendent constituents in a cold dark matter cosmology, and supplement solutions 
in the high and low-energy radiation-dominated phase under the assumption that 
both the dark energy and non-local anisotropic stress vanish. These solutions re- 
veal an additional pair of non-local isocurvature velocity and density modes which 
emerge from the two additional non-local degrees of freedom. These mode solutions 
would be useful in setting up initial conditions for numerical codes (for example, 
CAMB [126] and CMBFAST [175]) aimed at exploring the effect of braneworld 
corrections on the cosmic microwave background (CMB) power spectrum. For the 
tensor perturbations, we set out the framework of a program to compute the tensor 
anisotropies in the CMB that are generated in braneworld models. In the simplest 
approximation, we show the braneworld imprint as a correction to the power spec- 
tra for standard temperature and polarization anisotropies and similarly show that 
the tensor anisotropies are also insenitive to the high energy effects. 

In principle, future observations can constrain the KK parameter (, which 
controls the generation of 5D modes within our simplified local approximation, 
Eq. (4.151). The other braneworld parameter A, the brane tension, is not con- 
strained within our approximation. In practice, the tensor power spectra have not 
been measured, and the prospect of useful data is still some way off. What is more 
important is the theoretical task of improving on the simplified local approxima- 
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tion we have introduced. This approximation has allowed us to encode aspects 
of the qualitative features of braneworld tensor anisotropics, which we expect to 
survive in modified form within more realistic approximations. However, a proper 
understanding of braneworld effects must incorporate the nonlocal nature of the 
KK graviton modes, as reflected in the general form of Eq. (4.150). A possible 
direction is to investigate the tensor perturbations in the bulk and determine how 
it would actually affect the tensor power spectrum. 

It is also necessary to investigate the scalar anisotropies, which have a dominant 
contribution to the measured power spectra. There are a few possible directions 
for the scalar anisotropies. First, there is the possibility of exciting these non- 
local velocity and density isocurvature modes. Another possibility comes from 
the fact that the dark radiation, U is not necessarily zero. It has been shown by 
Langlois et al [119] that that the dark radiation, Uq is asymptotically a constant 
upon reaching the low energy regime, while it is time- dependent as U ~ a 4 at high 
energies. Hence it would be natural to extend the mode solutions we found for 
cases W ^ and work out the scalar anisotropies. All these possibilities for both 
scalar and tensor anisotropies may reveal new braneworld imprints that are more 
amenable to observational testing. 

Finally we come to the last theme of changing global symmetry. In chapter 
5, we constructed a possible relativistic description of how an Einstein-de Sitter 
expansion with no centre would evolve to an isothermal sphere with a centre. 
We have shown how this changing symmetry of the universe might possibly lead 
to a phase transition in the future. This construction demonstrates that general 
relativity permits a changing symmetry on a global scale, a nontrivial result. 

The largest galaxy cluster in a flat Einstein-de Sitter universe may grow in- 
definitely to encompass most of space after an extremely long time. We derive 
a general relativistic metric for a very large, bounded nearly isothermal cluster 
of galaxies which is embedded in such a universe. The embedding is done by 
means of a "Schwarzschild membrane". Pressure is important, unlike previous 
Tolman-Bondi models of inhomogeneities. The cluster's expansion, represented by 
a sequence of models, alters the average global symmetry of an increasing volume 
of space-time. Initially the metric is homogeneous and isotropic, having trans- 
lational and rotational symmetry around every point. As the metric evolves, it 
eventually loses its translational symmetry throughout larger regions and retains 
rotational symmetry around only one point: the centre of the cluster becomes the 
centre of the universe. Our sequence of hybrid models transforms between the 
Einstein-de Sitter and the isothermal cluster limits and illustrates how a changing 
equation of state of matter can alter the global symmetry of space-time. 

Gravitational clustering creates a cosmologically significant pressure, which can 
feed back into the spacetime structure of the universe. This induces a dynamical 
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symmetry breaking effected by the redistribution of matter. Is this transition be- 
tween distinct symmetry states of the universe an actual phase transition? The 
answer depends on the timescale and the definition of a phase transition. The 
experimental study of numerous condensed matter systems has demonstrated that 
transitions, which show a macroscopic discontinuity, might fundamentally be con- 
tinuous when examined on a microscopic scale (e.g. a ferromagnetic phase tran- 
sition). This may be analogous to the transition in the isothermal universe in 
which gravitational clustering provides the continuous, "microscopic" description 
of what is "macroscopically" a discontinuous phase transition when viewed on a 
long enough timescale. The different symmetries of the universe simply arise as 
different limits of the sequence of hybrid solutions. It has also suggested a new 
role for the cosmological constant, that is to deter the Universe from going into a 
phase transition in the future. 

This possible end of our Universe is conditional on the extremely long term 
stability of the matter. If the matter decays into radiation before the isothermal 
metric develops, the expansion can continue and possibly accelerate. Since matter 
can be converted into radiation in the isothermal universe, the equation of state 
would change and produce another possible phase transition. 

Of course, it is speculative to conceive the end of the simple Einstein-de Sit- 
ter Universe without a cosmological constant, just as it is difficult to conceive 
the beginning of our Universe that might be endowed with spin-torsion or extra- 
dimensions. It would be appropriate to conclude this thesis along the same thought 
as [171] that perhaps T. S. Eliot had some intuition of this idea in East Coker, 
where he said, "In my beginning is my end. " 
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Appendix A 



Conventions, Units, Sign 
Manifesto 

and Table of Symbols 

In this thesis, we would like to follow the spirit in [109] and [127] to state 
the conventions used and provide (i) the units for the conversion factors and fun- 
damental constants and (ii) cosmological parameters, the sign manifesto for the 
covariant approach (which is useful for the modifying the CAMB code) and the 
table of symbols adopted throughout this thesis. 



§A.l Conventions 

We would like to follow the conventions of the metric signature (— , +, +, +) 
for Chapters 1,3,4 and 5 and (+, — , — , — ) for Chapter 2, and we would adopt 
the choice of natural units, i.e. we set the fundamental constants h — c — ks — 1, 
and there is one fundamental dimension, energy, [E] such that 

[E] = [M] = [T E ] = [L] = [T]- 1 (A.l) 

where [M], [T E ], [L] and [T] are the dimensions of mass, temperature, length and 
time respectively. 

The unit of energy is taken to be a GeV= 10 3 MeV = 10 6 keV = 10 9 eV. For 
the case of braneworlds, the more popular unit for energy is in the TeV scale i.e. 
TeV= 10 3 GeV. 
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AND TABLE OF SYMBOLS 



jA.2 Units 

A. 2.1 Conversion Factors 

Temperature: 
Mass: 
Length: 
Time: 

Megaparsec: 
A. 2. 2 Fundamental 



Planck's constant: 


h = 


6.6261 x l(r 34 m 2 kg s" 1 


Speed of light: 


c = 


29972458 m s" 1 


Boltzmann's constant: 


k B ~- 


= 1.3807 x 10~ 23 J K- 1 


Newton's constant: 


G = 


: 6.672 x 10~ n kg m 3 s~ 2 = rripf 


Planck energy: 


mpi 


= y^M = 1.1211 x 10 19 GeV 


Planck mass: 


mpi 


= ^ = 2.1768 x 10~ 8 kg 


Planck time: 


tpi 


= x fi& = 5.3904 x lO" 44 s 


Planck length: 


I pi 


= ^§ = 1.6160 x 10~ 35 m 


Electron mass: 


m e 


= 0.5110 MeV 


Neutron mass: 


m n 


= 939.566 MeV 


Proton mass: 


m p 


= 938.272 MeV 


Thomas cross section: 


Gt 


= 6.652 x 10~ 29 m 2 


Solar mass 


M & 


= 2 x 10 30 kg 



jA.3 Physical Parameters 
A. 3.1 Cosmological parameters 



1 GeV = 1.1605 x 10 13 K 
1 GeV = 1.7827 x 10~ 27 kg 
1 GeV" 1 = 1.9733 x 10~ 16 m 
1 GeV" 1 = 6.5822 x 10~ 25 s 

1 Mpc = 10 6 pc = 3.0856 x 10 22 m = 1.5637 x 10 38 GeV" 1 
Constants 
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Hubble constant: H = 100 h km s" 1 Mpc" 1 

Present Hubble distance: cH 1 = 2998 h' 1 Mpc 

Present Hubble time: H Q l = 9.78 x lO 6 /^ 1 yr 

Present critical density: p c , = I.88/1 2 x 10~ 32 kg m~ 3 

§A.4 Sign Manifesto 

In order to incorporate the equations in the CAMB code, we need to re-write 
the equations in the convention (+, — , — , — ). To do this, we adopt the convention 
for the change of signs in [31] and expand the range of quantities in this thesis, for 
the following quantities: 

9ab —> —Qab, V a — > V a Rabcd ~^ Rabcd Rab ~^ Rab, 

R —> —R, h a b — > —hob, 1Z — > —TZ, D a — > —D a , 
u a —u a , q a — > — q a , A a — > —A a , a a b — > —c a b 
p^p, P^P, 0^0, U^U, 

A a — > — A a , Z a — > — Z a , T] a — > 7? a , 7T ab — > 7T a6 , 
T a ^ T a , Q a ► Qai ^a6 > Pab- 



(A.2) 



The above sign manifesto would help to translate equations from the signature 
of (— , +, +, +) to the other signature (+, — , — , — ). 



§A.5 Table of Symbols 

For all chapters using General Relativity 



Symbol 


Meaning of Symbol 


a(t) 


Scale Factor of the Universe with respect to cosmic time, t 


a(r) 


Scale Factor of the Universe with respect to conformal time, r 


H 


Hubble parameter with respect to cosmic time, t 


n 


Hubble parameter with respect to conformal time, r 


u a 


4- velocity/tangent vector 


v a 


Covariant derivative 
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n 

u a 


JrlUJcCLcCI I bpdlldl J LOVdlldlll QcllVdllvt; 


A 
1\ 


I AO T VT A 1 A (Tl / » o 1 f V\n OT O TIT" 

OOblilOlOglCdl OOIlbldlll 


A 


1 o n c l /"i n /~\T "f~ ri o r~»T*QT~io 

IcIlblOIl Ul Ulc Uldllt; 


D , , 

1 L abcd 


rLlclildllll IcIlbOl 


J^ab 


IxlCCl IcIlbOl 


D 
lh 


rUlXl bLdldl 


7? 
/V 


-< / - 1 iT*TroTii t*o o a n Iot l n /Toty ot* o 1 t*o 1 o fiintir 

O-CUI VdCUie bCdldl 111 geiieidl leidllVILy 




1 /~\ T Y O 1 O TY T T / - 1 1 Y*Y rO Til T*A T T Y fYO T Y AT' 1 O 1 T*0 1 fYTTT^TTtr 

OOIlStdllL OI CUI VdCUie Ill geiieidl leidtlVllTy 


^ab 


IH lTICITniTI TOTYOAT* 

rL/lIlbteill teilbOI 


J^ab 


1 H \~i OT* rfl r TYT / \ T VI f\v\ TIT WT T i ~\ Y~i O / \ T ' 

lilieigy-IIlOIIieilCUIIl teilbOI 


P 


1 1 OT~l nifir T O TT 111 /~1 1 /^lT* fYOTI OT* f\ 1 It r OT~l rWffT T /~1 O T Y nifir I 

ueiibiiy oi d nuiQ ^01 geneidiiy eneigy Qeiibiiyj 


A 


'■■i ~y~~\ n Ti A 1 rrv o rl l nTir i 1 1 t - n n t~\ r\ n o t t t r f av ati n-pm r nnncirir 1 at t ti o tH t t t a 

opdtldl gldQieilC ill lllH Qeilbliy ^OI eiieigy Geilbl iy J OI wie nuiQ 


p 
i 


1 'I'/ll 1 'lll'/^ AT O TT 111 n 

lr I ebbllie OI d I1U1Q 


Qa 


H/iit;igy nux oi d nuio 


T^ab 


IjULdl dlllbULIOpiL oLIcdd/ picboUIc 


Va 


OpdLldl gldUlcIlL OI lilt! O-CUI VdlUIc 


£>ab 


1 roinf n o opf I'm y \ .1 v| ^~\1 T"ri£i \A/£nrl fDncni 1 (hi "f"rio npo no 

VorldVllO-cltiLlIlL pdl I OI lllc Vvcyi LcIlbOI Oil II1K [Jldllc 


n ab 


1 t* n t n T a tyi rrn qt l a t~\ o t*t at Trio \A/oTrl Tono at* ati Trio hvoTici 

VjrldVlLO-IIldgllcllC pdl I Ul Lilts Vvcyi IcIIdOI Oil Lilt! Uldllc 


^ab 


\ / AT*T lAlTTr TOTIO AT* 

VOILlClLy LcIlbUI 


Cab 


Ollcdl 


V ) 


I / \ TYI / \ \ ' 1 Tl ( Y i \ \ r T "\ O T~l O 1 / \ Y~\ 

^yUIIlUVlllg, cXpclllolUll 




T "\ OTTO 1 (TT* 1 n 1 QTIT AT f lid A A TYI AT '111 fY 0"VT~\ QTICl ATI 

opdLldl ^IdCilclll Ul Lilt: CUIllUVlIl^ (^XpdllblUIl 


A 

^a 


-Ti-L/L/dt;! d-jlUll 


$ 


Grdvitdtioiidl Potentidl 


C s 


Adidbdtic bound speed 


n e 


Electron density 
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<7j< 


Thomson cross section 


Brane Worlds 


Symbol 


Meaning of Symbol 


S a b 


Symmetric tensor which contains the local quadratic energy correc- 
tions from the bulk to brane 




Symmetric projection of the bulk (5D)-Weyl tensor 


u 


Non-local energy density projected from bulk to brane 




Spatial gradient in non-local energy from bulk to brane 


Qa 


Non-local energy flux projected from bulk to brane 


v ab 


Non-local anisotropic stress projected from bulk to brane 



For the chapter using Gauge Theory of Gravity 



Symbol 


Meaning of Symbol 




Dirac vectors in Spacetime Algebra (STA) 


J 


Pseudoscalar 




Grade-r multivector 


V 


Vector derivative 


h(a) 


Position gauge field 


Q(a) 


Rotation gauge field 


V 


Covariant derivative with respect to a multivector 
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■?/) 


OpiilOl 111 vjr 1 \j 


T?(n A h \ 


OOVdllclill xvlclllcillll IcilbOl 111 Vjr 1 


/V-l tl i 


v^UVdlldlll rvi^L^l tclloUl 111 Vjt ± \j 


G{a) 


Covariant Einstein tensor in GTG 


T(a) 


Covariant energy momentum tensor in GTG 


S(a) 


Covariant spin torsion trivector 



* Gauge Theory of Gravity 



Appendix B 

Identities in Geometric Algebra 



In this appendix, we would like to summarize some of the important identities in 
geometric algebra. A more detailed exposition of the geometric algbera formalism 
is given in Hestenes et al[85, 86, 87, 88] and Lasenby, Doran and Gull in [121]. 

§B.l Mult i vectors 

Let A and B be multivectors of grade r and s, • be the interior product, A be 
the exterior product and x be the commutator product, where these operations are 
defined in Eqs. (2.2), (2.3) and (2.4) respectively in chapter 2, then the following 
identities hold: 




(B.l) 



r 



A r B s = (A r B s )\ r - s \ + (A r B s )\ r _ s \ + 2 + + (A r B s )\ r+s \ , 

A r ■ A = , 

(AB) = (BA) , 

A r ■ (B s ■ C t ) = (Ar A B s ) -C t (for r + s < t and r, s > 0) , 
A r ■ (B s ■ C t ) = (A r ■ B a ) ■ C t (forr + t<s), 
A x (B x C) + C x (A x B) + B x (C x A) = . 



(B.2) 
(B.3) 
(B.4) 
(B.5) 
(B.6) 
(B.7) 



With the 



pseudoscalar, X, we have the following identities, 



A r ■ (B S T) = A r A B S X (r + s < n) , 
A r A (B S T) = A r ■ B S X (r < s) . 



(B.8) 
(B.9) 



We now specialize to the case for the vector and bivector cases. Let {a, b, c} 
be a set of vectors and {A, B, C} be a set of bivectors, and the following identities 
will hold: 
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a ■ {b A c) = a ■ be — a ■ cb , (B.IO) 

(a Ab) • (c Ad) = a ■ db ■ c - a ■ cb ■ d , (B.ll) 

a ■ (b ■ B) = (a A b) ■ B , (B.12) 

(a A 6) x (c A d) = b ■ ca A d - a ■ cb A d + a • db A c - b ■ da A c , (B. 13) 

(aAb)xB = (a-B)Ab + aA(b-B), (B.14) 

BC = B-C + BxC + BAC , (B.15) 

CB = B-C-BxC + BAC, (B.16) 

B x A r = (B x A r ) r , (B.17) 

A x (.BC) = AxBC + BAxC. (B.18) 



§B.2 Geometric Calculus 

Let F be an arbitrary function of some multivector argument X, such that 
F = F(X). The derivative of F with respect to X in the A-direction is defined to 
be 

A tSx F(X) = ^ F< - X + xA) - F< - X) (B.19) 
where we define the multivector derivative dx to be 

d x = Yl ei A • • • A ei ( e ; A • • • A e *) * dx ( B - 20 ) 

i<-<j 

Most of the properties of the multivector derivative follow from the result 

dx(XA) =P X (A) (B.21) 

where Px{A) is the projection of the multivector A onto the graded element con- 
tained in X. The multivector derivative acts on the next object to its right unless 
brackets are present. Denoting the multivector which the derivative acts, we can 
write Leibniz rule in the form: 

dx(AB) = d x AB + d x AB (B.22) 
The derivative with respect to the vector a, d a is used to perform linear alegbra 



§B.3. LINEAR ALGEBRA 
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operations. The following identities are obtained as follows: 



d a b = b , 




(B.23) 


d a a = n , 




(B.24) 


d a A a = , 




(B.25) 


d a a = n , 




(B.26) 


d a a 2 = 2a , 




(B.27) 


(B ■ d a ) A a = 2B , 




(B.28) 


d a A d b (b A a) ■ B = 


2B, 


(B.29) 


d a A r a=(-l) r (n- 


2r)A r , 


(B.30) 


dad ' A r = v A T ; 




(B.31) 



§B.3 Linear Algebra 

Let a, b, c be vectors, f(A) be the linear function with the multivector A in its 
argument and f(A) be its adjoint, from the definition in Eq. (2.9), it follows that 

f(a) = d b (h(b)a) (B.32) 

Eq. (B.32) is frequently employed perform algebraic manipulations of linear func- 
tions. The advantage is that all manipulations are frame-frame. We note that the 
d a and a vectors can be replaced by the sum over a set of constant frame vectors 
and their reciprocals, if required. 

Linear functions extend to act on multivectors via 

f(a A b A ... A c) = 1(a) A /(6) A ... A /(c) (B.33) 

so that / is a grade-preserve linear function mapping multivectors to vectors. 
In particular, since the pseudoscalar X is unique to a scale factor, we define 



det(/) = /(J)J- 1 



(B.34) 
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From the Eq. (B.33) and (B.34), the following identities hold: 



/ (a A b A ... A c) = f(a) A f(b) A ... A f(c) , 


(B.35) 


f(I) = det(/)J , 


(B.36) 


f(b) = d a f(a) ■ b , 


(B.37) 


/(A) • B s = f[A r ■ f(B s )] (r>s), 


(B.38) 


A r ■ f(B s ) = f[f(A r ) ■ B s ] (r<s), 


(B.39) 


f-\A) = det(fr 1 lf(l- 1 A), 


(B.40) 


f-\A) = det(f)- 1 lf(I- l A), 


(B.41) 


d^deth^ 1 = —dethr 1 h~ 1 (a) , 


(B.42) 


d- K a)h-\b)-v = -h-\a)-vh-\b) . 


(B.43) 



Appendix C 

Identities in Cosmological 
Perturbation Theory 



C.l Differential Identities in the PSTF Approach 



We will summarize the following covariant linearized identities on a Friedmann 
background used in [31, 134, 135, 137, 138], which are useful for the derivation of 
the propagation and evolution equations in chapter 3. 



The following linearized identities with the operation of curl and div, which 
is valid for any arbitrary first order PSTF vector, V a , tensor, S a b, and arbitrary 
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scalar / with first order projected gradients are summarized as follows: 



curlAJ = -2fu a , (C.l 

2K 

D 2 (DJ) = D a (D 2 f) + — D a f + 2/curL; a , (C.2 

(X 

(Djy =DJ-HDJ + fA a , (C.3 

(aD a V b ...y =aD a V b ... , (C.4 

(D 2 fy = D 2 f - 2HD 2 f + fD a A a , (C.5 
K 

D [a D b] V c = —V [a h b]c , (C.6 

(X 

2K 

D [a D b] s cd = -rS [a ^h b f , (c.r 

D a cm\V a = , (C.8 

2K 

curlcurlK = D a (D b V b ) - D 2 V a + —V a , (C.9 

a 2 

(curlS^)' = curlS^ - -OcurlS^ , (CIO 

(D a S ab )' = D a S ab — ^QD a S ab , (C.ll 

D a (cm\S ab ) = ^cuT\(D a S ab ) , (C.12 

3 3K 

curl(curlS a5 ) = —D 2 S ab + -D {a D c S b)c + — S ab , (C.13 

Zi a 



where the vectors and tensors vanish in the background, i.e. S ab = and all 
the identities except (C.l) are linearized. 

The identities in [31] are also useful for the evaluation of the curl quantities for 
parity tensor harmonics are also summarized here: 



curlQ a6 
curlQ afe 



* 


(■♦ 


3K 






a 


T? 


* 




3K 






a 





1/2 



(k) 



1/2 



Q 



ab 



(C.14) 
(C.15) 
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§C.2 The identitites for the harmonics 

§C.2.1 Scalar Harmonics 

As in [22, 79, 108], we define the scalar harmonics as being eigenfunctions of 
the covariantly Laplace-Beltrami operator. 



D 2 Q= DjrQ [-) Q, (C16) 



where we have dropped the index k from the Qs. One can use the scalar harmonic 
to expand scalars, and with this we define the following vector which correspond 
to the definition in [108]. 

Q a = ~D a Q, (C.17) 
and a trace-free symmetric tensor, 

Q ah = D b D a Q+ l -h ah Q . (C.18) 

These harmonics are defined to have the following feature: 

Q = Qa = Qab = 0, (C.19) 

and by the above property, we can apply the commutations in [66] to get 

D [a D b] Q = -u ab Q = , (C.20) 

D [a D b] Q c = ~^(h ac Q b - h bc Q a ) , (C.21) 
K 

D[ a D b }Q cd = -^[(hacQbd — h bc Q ad ) + (h ad Q bc — huQac)] ■ (C.22) 



We derive the properties of the scalar harmonics, by making use of the projected 
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(spatial) covariant derivative, D a , 



D a Q a ( C - 23 ) 

D 2 Q a = - (^^) Qa , (C24) 

D b Q a = ~ (Qab - \h ab Q^j , (C.25) 

Q\ = , (C.26) 

D b Q ab = 1 (^) (k 2 - 3K)Q a , (C.27) 



D a D b Q ab = ) (k 2 - 3K)Q , (C.28) 



D b D c Q ac = \ (3K - k 2 ) (Q ab - \h ab Q ) , (C.29) 



3 Va 2 / V 3 

/ t\K — k 2 \ 

D 2 Q ab =1 q2 j Q ab . (C.30) 



§C.2.2 Vector Harmonics 

These harmonics are defined as eigenfunctions of the Helmholtz equation: 

D 2 Q a = -(^j Q a , (C.31) 

with Q a being the solenoidal vector, i.e. 

D a Q a = . (C.32) 
From this, we can construct a symmetric trace-free tensor, 

Qab = -^( D bQa + D a Q b ) , (C.33) 

and we note that both these harmonics are covariantly constant along the flow 
lines 

Qa = Qab = . (C.34) 

The following relations hold for the defined tensor Q ab (equation (C.33)) from 
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Q a in this section: 

Q\ = , (C.35) 

D b Q ab = l - (j^j (k 2 - 2K)Q a , (C.36) 
{ k 2 - 2K\ 

D b D c Q ac + D a D c Q bc = - f q2 j Q ab , (C.37) 

D 2 Q ab = - (^-=^) • (C38) 
Note that Q a and Q ab also satisfy relations (C.21) and (C.22) respectively. 

§C.2.3 Tensor Harmonics 

For tensor harmonics, we define again 

D 2 Q ab = -[^j Q ab , (C.39) 
and Qab = 0. The tensor is also trace- free and divergenceless, 



Q\ = , (C.40) 
D b Q ab = . (C.41) 
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